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EXECUTIVE  SUMMARY 


Optimal  Design 

Engineering  design  of  a  product  (such  as  an  aircraft)  invariably  seeks  to  achieve  several 
desiderata  -  several  features  and  performances  desired  by  the  customer  (actual  or 
perceived  customer).  Here,  a  “feature”  means  a  characteristic  that  the  product  either  will 
have  or  won’t  have  (such  as  the  capability  to  land  on  a  carrier),  while  a  “performance” 
indicates  an  amount  of  a  characteristic  that  the  product  must  have  (such  as  top  speed). 

Usually  the  desiderata  are  conflicting,  in  the  sense  that  improvement  of  one  requires 
some  sacrifice  of  another.  For  example,  one  typical  desideratum  of  a  vehicle  is  low  cost 
and  another  is  high  fuel  economy.  But  improvement  of  fuel  economy  entails  greater  cost. 

Therefore,  design  requires  that  a  means  be  employed  for  “balancing”  the  various 
desiderata,  in  order  to  determine  and  then  account  for  their  relative  importance. 
Traditionally,  this  has  been  accomplished,  implicitly  if  not  explicitly,  by  “jawboning”  - 
with  real  or  hypothetical  advocates  for  each  desideratum  arguing  for  theirs,  before  a 
“judge”  or  a  “court”  that  decides.  Such  a  procedure  has  been  effective  in  the  past,  even 
though  it  entails  much  subjectivity  and,  consequently,  takes  a  lot  of  time.  So,  a  rational 
means  is  needed  for  taking  the  various  desiderata  into  account. 

Optimal  design  seeks  to  provide  such  a  means.  The  key  is  to  establish  a  mathematical 
“objective  function”  with  a  measure  of  value  (to  the  customer)  as  the  dependent  variable 
and  the  desiderata  as  independent  variables.  A  common  measure  of  value  is  money 
(“dollar-value”).  However,  in  the  case  of  military  aircraft  -  where  the  total  budget  for  the 
“fleet”  to  be  produced  is  dictated,  not  by  the  customer  (say  the  Department  of  Defense), 
but  by  a  “third  party”  (Congress)  -  the  measure  of  value  might  well  be  taken  to  be  the 
size  of  the  fleet  (i.e.,  the  number  of  aircraft  to  be  produced,  the  “production”). 

More  difficult  than  the  selection  of  the  measure  for  value  is  the  establishment  of  the 
function  -  that  is,  the  relationship  that  determines  the  influence  upon  value  of  (a)  the 
existence  of  a  feature,  and  (b)  amounts  of  each  performance. 

This  report,  within  the  context  of  a  general  outline  for  thermal  design,  introduces  a 
procedure  for  establishing  an  objective  function,  rationally  accounting  for  each 
desideratum,  weighting  their  relative  importance  and,  in  the  case  of  a  performance, 
quantifying  the  value  of  varying  amounts  thereof  This  procedure  is  illustrated  by  an 
elementary  example,  applicable  to  certain  choices  made  in  the  design  of  a  light  aircraft. 

In  particular,  to  design  decisions  regarding  selections  -  choices  from  among  alternatives 
(which  differ  in  performance  or  features)  -  for  an  engine,  a  generator,  and  a  propeller. 

Subsystem  Decomposition 

When  designing  a  system  such  as  an  aircraft,  the  design  encompasses  several  subsystems, 
such  as  engine,  propeller,  generator, ....  For  a  simple  system,  like  a  light  aircraft,  the 
decisions  regarding  individual  subsystems  can  be  guided  by  direct  application  of  the 
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overall,  system  objective  function.  That  is  what  was  done  in  the  example  referred  to 
above. 

When  designing  a  complex  system  such  as  a  combat  aircraft,  in  practice  there  is  the  need 
for  “concurrent  engineering.”  The  individual  designs  of  different  subsystems  are  carried 
out  concurrently  by  distinct  engineers  (or  teams)  -  each  with  their  own  expertise. 
Unfortunately,  independent  optimization  of  the  subsystems,  separately,  in  no  way  insures 
optimal  design  of  the  overall,  complex  system. 

Again,  traditionally,  the  overall  “optimization”  has  been  carried  out  subjectively, 
depending  upon  the  experience  and  judgment  of  a  chief  “project”  engineer  to  make 
decisions  regarding  parameters  that  reflect  subsystem  interactions.  Here  too, 
“jawboning”  has  been  an  influence,  with  each  team  leader  advocating  parameter  values 
that  would  be  advantageous  to  the  team’s  subsystem. 

A  procedure  that  would  streamline  the  decision-making  and  make  it  objective,  is  called 
for.  Such  a  method  would  need  to  decompose  the  overall  objective  function  into 
individual  objective  fimctions  for  the  subsystems,  in  a  manner  such  that  somehow 
optimization  of  the  subsystem  objectives  would  assure  overall  optimization. 

This  report  also  presents  a  rational  methodology  for  “subsystem  integration.”  -  a 
procedure  for  establishing  subsystem  objective  functions  such  that,  when  each  is 
individually  optimized,  the  overall  objective  function  is  optimized.  For  each  subsystem 
objective,  Ae  procedure  deduces  weighting  and  cost  factors  from  the  basic,  overall 
objective  function.  This  methodology  is  also  illustrated  by  application  to  the 
aforementioned  light  aircraft.  The  application  of  the  methodology  requires  a  means  for 
assigning  costs  to  any  “energy”  flows  to  and  from  subsystems,  to  be  discussed  next. 

Thermoeconomic  Costin2 

In  the  overall  objective  for  an  aircraft,  fuel  consumption  is  invariably  a  factor;  one 
desideratum  is  always  to  reduce  fuel  consumption.  In  one  way  or  another,  evciy 
subsystem  consumes  “energy”  -  that  is,  every  subsystem  is  a  cause  of  fuel  consumption. 
Some  aircraft  subsystems  use  energy  only  passively,  energy  required  for  “lift”,  to  keep 
the  subsystem  aloft,  and  energy  required  for  acceleration.  Other  subsystems  actively 
consume  energy,  in  one  form  or  another,  and  many  of  these  also  deliver  energy,  in  a 
different  form.  The  design  objective  function  for  any  subsystem  needs  to  debit  the 
subsystem  for  every  imit  of  energy  consumed  and  give  a  credit  for  any  energy  delivered. 
The  value  of  a  unit  of  energy,  say  a  kilowatt-hour,  depends  upon  how  much  fuel  has  been 
consumed  to  generate  the  unit  and  upon  other,  associated  investments,  such  as  the  cost  of 
all  the  equipment  employed.  In  order  to  follow  all  the  costs  through  the  complex,  overall 
aircraft,  and  assign  imit  values  at  different  junctures,  the  proper  commodity  -  the  proper 
“energy”  to  be  tracked  -  is  not  textbook  energy,  but  the  concept  that,  today,  is  called 
exergy.  Exergy  (also  called  “availability”)  is  the  appropriate  measure  of  marketplace 
energy  -  that  is,  of  “potential  to  drive  processes.”  It  is  the  key  not  only  for  proper  cost 
accounting  of  “energy”  but  also,  even  beforehand,  to  analysis  -  i.e.,  to  pinpointing  and 
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quantification  of  inefficiencies  in  energy-conversion  systems.  (And,  basically,  an  aircraft 
is  one  example  of  an  energy-conversion  system.) 

In  this  report,  the  overall  objective  function  for  a  light  aircraft  is  decomposed  into 
subsystem  objective  functions,  and  exergy  costing  is  then  employed  to  account  for  debits 
and  credits  associated  with  “energy”  flows.  Each  of  the  subsystems  is  then  optimized 
independently.  And,  it  is  shown  that  the  subsystem  optima  that  are  achieved  lead  to  the 
same  decisions  as  those  made  without  decomposition  -  the  same  as  the  decisions  reached 
by  applying  the  overall  objective  function  alone.  This  comparison,  possible  because  of 
the  simple  subsystem  applications  within  the  context  of  a  relatively  simple  overall  system 
-  the  light  aircraft  -  provides  evidence  that  backs  up  (a)  the  precede  for  establishing  the 
overall  objective  function,  (b)  the  procedure  for  subsystem  composition,  (c)  the 
thermoeconomic  methods  of  exergy  costing,  and  (d)  the  means  developed  to  account  for 
“the  exergy  of  lift”. 


Exergy  of  Lift 

Clearly,  the  weight  of  any  subsystem  penalizes  the  performance  of  an  aircraft,  increasing 
the  fuel  consumption  and  decreasing  the  acceleration.  In  order  to  employ  exergy  costing 
to  achieve  subsystem  decomposition,  it  was  evident  that  a  means  was  needed  to  account 
for  weight.  This  led  to  a  new  theoretical  development,  the  deduction  of  a  new  component 
of  exergy,  needed  for  exergetic  analysis  and  costing  in  the  context  of  aircraft. 


Summary  of  Results  Obtained 

1.  A  procedure  for  the  development  of  an  overall  objective  function  for  rational 
optimal  design  of  an  aircraft,  taking  into  account  all  of  the  desiderata  and  their 
relative  importance. 

2.  A  procedure  for  “concurrent  engineering”  of  subsystems  -  “subsystem 
integration”  -  that  develops  an  objective  function  for  each  energy-conversion 
subsystem  (using  the  overall  objective  function  to  derive  weighting  factors 
employed  in  the  subsystem  objectives)  so  that  subsystem  optimization  leads  to 
overall  aircraft  optimization. 

3.  The  application  of  exergetic  costing  -  “thermoeconomics”  -  in  order  to  properly 
assign  the  impact  on  fuel  consumption  in  each  of  the  subsystem  objective 
functions  (i.e.,  the  costs  of  fuel,  of  energy-conversion  equipment  capital,  of 
energy-conversion  equipment  inefficiencies,  of  energy-conversion  equipment 
weight). 

4.  The  development  of  theoretical  relationships  for  evaluating  a  new  component  of 
exergy  -  “exergy  of  lift”  -  needed  to  properly  account  for  and  track  the  exergy 
costs  (i.e.,  the  costs  associated  with  energy-conversion  equipment  weight). 

Together,  the  foregoing  developments  establish  an  overall,  rational  procedure  for 

vehicle  design.  The  results  are  developed  in  detail  in  the  following  report  and,  more 

concisely,  in  the  technical  papers  appended  to  the  report. 
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Recommendations 

•  General.  It  is  suggested  that,  in  order  to  further  test  and  perfect  them, 

o  The  methods,  which  so  far  have  been  applied  only  to  the  selection  of 
subsystems,  be  applied  to  the  detailed  design  of  the  subsystems. 

o  The  methods,  which  so  far  have  been  applied  only  to  the  simple 
subsystems  of  a  simple  aircraft,  should  be  applied  to  the  design  of  a 
complete  aircraft  -  say  a  relatively  simple  one  such  as  an  Unmanned 
Reconnaissance  Vehicle  -  including  all  of  its  subsystems. 

•  Specific.  It  is  suggested  that  the  following  refinements  be  pursued. 

o  During  the  development  of  the  overall  objective  function,  more 

sophisticated  relations  might  well  be  employed  for  representation  of  the 
importance  of  the  various  desiderata. 

o  During  the  development  of  the  overall  objective  function,  the  “interaction” 
of  desiderata  could  be  taken  into  account.  That  is  the  “representations  of 
importance”  could  include  more  than  one  independent  variable,  more  than 
one  desideratum. 

o  During  the  application  of  exergetic  costing,  the  determination  of  marginal 
costs  in  lieu  of  average  costs  needs  to  be  pursued  further. 
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i.  NOMENCLATURE 


Variable 

Use 

a 

aspect  ratio 

A 

area 

C 

heat  capacity 

c 

cost 

Cd,Cl 

coefficient  of  drag,  lift 

E 

energy 

F 

fuel 

F 

force 

g 

acceleration  of  gravity 

J 

objective  function 

m 

mass 

p 

performance 

R 

range 

S,s 

entropy,  specific  entropy,  sensitivity  coefficient 

T 

time  factor,  temperature 

V 

value 

V 

speed 

V 

velocity  (scalar) 

w 

weighting  factor 

X,x 

exergy,  specific  exergy 

. .y _ _ 

decision  variable 

z 

capital 

Table  1:  Variables 


Greek  Variable 

Use 

77 

profit 

. . M  ....  . .  ^ 

viscosity 

_ ^ _ a _ 

density 

Table  2:  Greek  Variables 


2 


Subscript 

Use 

DF 

duct  firing 

HRSG 

heat  recovery  steam  generator 

m 

mass 

s 

destmction 

n 

production 

Table  3:  Subscripts 


il.  MOTIVATION 

This  work  was  performed  with  Air  Force  Office  of  Research  Support  funding  for  the 
optimization  of  aircraft  energy  conversion  subsystems.  The  optimization  of  these  subsystems 
proves  to  be  more  intricate  than  many  traditional  energy  system  design  problems,  as  figures  of 
merit  for  the  vehicle  as  a  whole  are  difficult  to  translate  into  values  that  will  be  meaningful  to  an 
on-board  energy  system  (or  subsystem  or  device).  These  figures  of  merit  are  manifold  for  a 
modem  fighter,  as  it  must  excel  in  several  competing  performance  areas.  One  objective  is  to  have 
a  low  signature,  i.e.,  be  stealthy.  Another  objective  is  rapid  acceleration,  typically  measured  from 
Mach  0.8  to  Mach  1 .2.  The  aircraft  should  have  a  long  range  and  a  great  payload. 

Here  these  problems  will  be  addressed  in  two  manners.  First,  an  overall  design 
methodology,  from  the  feasibility  study  through  detailed  design,  for  thermal  systems  will  be 
developed  in  a  general  fashion.  Secondly,  the  detailed  design  phase  will  be  covered  in  detail.  In 
order  to  support  these  two  areas,  a  means  of  developing  an  overall  objective  function  for  a  vehicle 
with  multiple  objectives  will  be  given.  Additionally,  it  wiU  be  shown  how  to  develop  subsidiary 
objective  functions  for  systems,  subsystems  and  devices  from  this  overall  function. 

Although  particular  attention  will  be  paid  to  vehicular  energy  systems,  namely  those  for 
aircraft,  the  generd  approach  for  thermal  system  design  outlined  here  will  be  as  valid  for  a 
stationary  plant  as  for  a  vehicle.  The  approach  will  be  boimded  in  the  following  manner.  It  is 
assumed  that  a  need  for  an  energy  system  has  been  given.  Only  thermal  design  will  be 
considered  explicitly.  Maintainability,  reliability,  architecture  and  mechanical  design  are  beyond 
the  scope  of  this  work.  Whenever  possible,  the  overall  design  processes  as  laid  out  in  Ostrofsky, 
1977  and  Woodson,  1966  are  followed.  This  is  to  allow  the  thermal  design  process  here  to  share 
the  same  vocabulary  and  jargon  with  general  design  as  a  standard.  In  many  instances,  however, 
the  steps  and  substeps  given  in  these  are  not  are  not  applicable,  or  out  of  place,  for  the  thermal 
design  process.  In  such  cases  they  are  not  used,  or  occur  at  different  times. 

For  the  design  process  to  occur  in  a  rational,  quantifiable  fashion,  a  mathematical 
objective  function  is  necessary.  Current  practice  in  vehicular  design  is  to  list  all  of  the  desired 
objectives  (e.g.,  speed,  payload,  climb  rate,  acceleration,  etc.)  and  then  to  express  some  of  them 
in  quantitative  terms  and  some  solely  in  qualitative  terms.  But  there  is  not  an  objective  procedure 
for  weighting  the  relative  importance  of  these  often-competing  desiderata.  The  “weighting”  is 
carried  out,  subjectively,  during  the  course  of  the  design  process,  primarily  as  a  consequence  of 
conversations  between  different  involved  personnel  —  such  as  the  customer,  the  end-users,  the 
chief  oigineer,  the  various  specialized  design  teams,  etc. 

While  such  conversations  will  always  be  important,  in  order  to  expedite  the  design 
I  process  and  make  it  more  efficient,  it  is  desirable  to  establish  a  rational  objective  function  for 
each  overall  vehicle  project,  including  a  quantitative  means  for  assessing  the  importance  of  the 
various  desired  objectives. 
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Extensive  work  has  been  done  in  the  area  of  optimized  design  of  energy  systems  that 
have  an  end  purpose  of  providing  a  mass,  heat  or  power  flow.  Examples  include  the 
“Thermoeconomics”  of  El-Sayed  and  co-workers  (e.g.,  El-Sayed  and  Evans,  1970,  Evans  and  von 
Spakovsky,  1984).  However,  many  energy  systems  are  used  in  applications  for  which  the  end 
product  is  not  energy,  such  as  vehicular  applications.  In  such  instances  the  energy  systems 
impact  not  only  lifetime  costs  of  the  apphcation,  but  also  its  performance.  It  is  necessary  to  take 
these  performance  impacts  into  account  in  order  properly  to  optimize  these  energy  systems. 
Indeed,  for  some  vehicles,  such  as  a  fighter  or  a  Formula-1  automobile,  performance  is  the 
overwhelming  objective.  Often  diese  performance  desiderata  are  competing,  as  per  the  example 
of  the  fighter  aircraft  above.  The  determination  of  the  optimum  balance  between  these  desiderata 
is  naturdly  a  greater  challenge  than  many  traditional,  stationary  energy  system  design  problems. 

Detailed  design  and  optimization  will  be  covered  in  depth.  It  is  in  this  phase  that  second 
law  techniques  enable  a  design  team  to  approach  an  optimal  design  while  concurrently 
performing  detailed  design  and  selection  of  components  to  an  energy  system.  The  application  of 
the  second  law  of  thermodynamics,  combined  with  economics,  allows  for  decomposition  of 
systems  into  devices  and  subsystems  that  may  be  independently  optimized.  This  theory  has  been 
developed  and  validated  by  El-Sayed  and  coworkers.  (See  for  example,  El-Sayed  and  Evans, 
1970;  El-Sayed,  1995.) 

Additional  second  law  and  thermoeconomic  considerations  will  be  developed  as 
necessary  to  support  the  design  of  vehicular  energy  systems.  These  new  developments  include 
die  exergy  of  lift,  the  flow  of  exergy  (and  its  cost)  in  an  aircraft,  the  reference  environment  and 
marginal  exergy  costing.  (Naturally  the  latter  two  areas  have  much  broader  applications  than 
solely  vehicle  design.) 

To  illustrate  the  concepts  of  this  report,  an  overall  objective  function  for  a  light  aircraft 
will  be  fonnulated.  The  exergy  flow  in  this  aircraft  will  be  diagrammed  for  several  realms  of 
flight.  Then  cost  balances  will  be  applied  to  find  the  marginal  costs  per  unit  exergy.  Finally,  an 
engine  and  an  alternator  will  be  selected  for  this  aircraft  using  this  information,  for  optimal  design 
of  the  aircraft. 

III.  DEFINITIONS 

(a)  System 

An  energy  conversion  system  is  a  device  or  set  of  devices  employed  to  produce  energy 
products  (i.e.,  exergy  in  some  form)  at  the  expense  of  the  consumption  of  fuel  (exergy  in  a 
different  form.  Here  the  word  system  is  used  for  the  largest  grouping  of  equipment  with  pure 
exergy  inputs  and  outputs.  One  example:  A  combined  cycle  power  plant.  For  aircraft,  a  system 
will  be  considered  a  collection  of  devices  used  to  provide  a  fundamental  requirement  for  flight. 

An  aircraft  example  of  a  system  would  be  the  propulsion  system. 

(b)  Subsystem 

An  energy  conversion  subsystem  is  a  set  of  devices,  contained  within  a  system, 
configured  togetiier  to  produce  products  used  within  a  system.  The  realms  of  the  system  and 
subsystem  are  not  clear-cut;  however,  the  subsystem’s  products  and  fuels  are  dependent  on  the 
requirements  of  the  overall  system.  Examples  in  a  combined  cycle  power  plant:  the  gas  turbine 
engine  and  the  heat  recovery  steam  generator  (HRSG).  For  aircraft,  this  report  will  define  an 
energy  conversion  subsystem  as  an  energy  conversion  system  not  directly  involved  in  the 
production  of  thrust.  One  aircraft  example  is  an  enviromnental  control  system. 
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(c)  Device 

An  energy  conversion  device  is  defined  here  as  hardware,  normally  contained  in  a  single 
housing,  designed  to  produce  a  single  product  with  only  one  or  two  fuels.  Examples  in  the 
combined  cycle  power  plant:  within  the  gas  turbine,  the  compressor,  combustion  chamber(s)  and 
the  turbine;  within  the  HRSG  the  various  heat  exchangers;  in  the  remainder,  pumps  and  steam 
turbines.  Examples  in  an  aircraft:  the  compressor,  a  hydraulic  pump,  an  alternator. . . 

(d)  System-level  Decision  Variable 

A  system  level  decision  variable  is  a  variable,  that,  when  changed,  will  greatly  affect  the 
(optimal)  design  and/or  selection  of  subsystems  and  devices.  In  the  combined  cycle  power  plant 
system  level  variables  may  include  gas  turbine  pressure  ratio,  firing  temperature  and  steam  side 
operating  pressines. 


(e)  Device-level  Decision  Variable 

A  device  level  variable,  when  perturbed,  has  little  or  no  effects  on  optimal  values  of 
system  level  or  other  device’s  decision  variables.  One  example  of  a  device  decision  variable  in  a 
combined  cycle  power  plant  would  be  the  peak  efficiency  of  a  HRSG  feed  pump. 

(f)  Fuel 

Any  exergy  input  to  a  system  or  subsystem  is  considered  a  fuel.  This  is  beyond  the 
traditional  thought  of  &el  being  something  “burned”  as  an  energy  source.  Of  coiuse,  the  exergy 
source  for  any  device  is  at  least  partially  “burned”  or  depleted  driving  whatever  process  occurs  in 
the  device. 


(g)  Product 

The  desired  exergy  outputs  of  an  energy  system  or  subsystem  are  products.  The  product 
of  a  vehicle  is  performance:  speed,  acceleration,  payload,  etc. 

(h)  Byproduct 

A  byproduct  is  an  exergy  output  of  a  system  that  is  not  necessarily  desirable  but 
unavoidable  due  to  the  process. 

(i)  Candidate  System 

A  candidate  system  is  a  system  that  satisfies  the  laws  of  physics  and  the  demands  of  the 
problem. 

(j)  Constraint 

A  constraint  limits  the  designer’s  fireedom  in  conceiving  a  design.  These  constraints  may 
be  on  the  design  itself,  e.g.  the  design  must  fit  in  a  certain  place,  or  on  the  design  process.  A 
typical  constraint  on  the  design  process  itself  is  time. 

iV.  A  GENERAL  OUTLINE  OF  THE  DESIGN  PROCESS  FOR  THERMAL 
SYSTEMS 

As  presented  here,  the  design  of  thermal  systems  has  three  phases.  In  the  first,  the 
feasibility  study,  a  set  of  possible  solutions,  consisting  of  candidate  systems,  is  developed.  In  the 
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next  phase,  preliminary  design,  the  best  design  of  the  candidate  systems  is  chosen  for 
development.  Finally,  during  detailed  design  all  hardware  is  selected  and  decision  variable 
concretized. 

(a)  Feasibility  Study 

According  to  Ostrofsky,  “The  purpose  of  the  feasibility  study  is  to  develop  a  set  of  useful 
solutions  to  meet  the  needs. . .  The  primary  importance  of  the  feasibility  study  is  that  it  is  the 
foundation  for  all  that  occurs  subsequently  in  the  design-plan  for  the  system.”' 

Applied  to  thermal  design,  the  feasibility  study  should  generate  a  set  of  solutions,  all  of 
which  meet  the  design  needs,  that  (a)  satisfy  the  laws  of  diermodynamics  and  (b)  are 
technologically  achievable  within  any  time  constraints  that  exist.  This  is  achieved  through  (i)  a 
product,  fuel  and  byproduct  analysis  which  yields  (ii)  a  problem  formulation.  After  these  are 
complete,  a  set  of  possible  solutions  are  proposed,  and  screened  so  that  they  meet  the  two 
conditions  stated  above. 

(i)  Product,  Fuel  and  Byproduct  Analysis 

Before  any  solutions  may  be  conceived,  the  energy  system  designer  must  know  (a)  what 
the  system  is  supposed  to  produce  (the  product(s)),  (b)  what  is  available  for  the  system  to 
consume  (the  fuel(s))  and  (c)  the  possible  byproducts  (both  undesirable  and  possibly  desirable). 

Typical  products  for  an  energy  system  include,  but  are  not  limited  to,  shaft  power, 
heating  or  cooling  loads  and  compressed  air  or  other  mass  flows. 

Fuels  may  be  fuels  in  the  conventional  sense,  such  as  fossil  or  nuclear  fuels,  or  mass,  heat 
or  work  streams  from  another  system  or  subsystem. 

At  this  stage  in  the  design  process,  undesirable  byproducts  should  be  listed.  These 
include,  among  others,  emissions,  noise,  weight,  volume  and  drag. 

(ii)  Problem  Formulation 

With  the  fuels,  products  and  byproducts  listed  the  problem  statement  may  be  crafted.  In 
addition  to  stating  the  information  generated  from  product,  fuel  and  byproduct  analysis,  any 
constraints  on  the  system  must  be  included.  Examples  of  such  constraints  are  output  (product) 
requirements  and  size,  emissions  and  noise  limitations.  Additionally,  time  constraints  on  the 
design  should  be  noted.  These  may  have  a  large  impact  on  the  screening  of  the  solutions 
generated  in  the  next  phase  of  die  feasibility  study. 

(Hi)  Solution  Generation 

Solution  generation  involves  the  “brainstorming”  of  possible  configurations  and 
processes  to  meet  the  problem.  At  this  stage,  no  consideration  should  be  given  to  the  viability  of 
the  brainstormed  solutions.  The  end  product  of  solution  generation  is  a  set  of  candidate  systems. 


(iv)  Solution  Screening 

Solution  screening  involves  the  asking  the  following  questions: 


*  Ostrofsky,  1977,  pg.  29 
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1 .  Is  the  system  physically  attainable?  That  is,  can  the  system  meet  the  needs  of  the 
problem  without  violating  the  laws  of  physics?  For  energy  systems,  the  candidate 
system  must  obey  the  first  and  second  laws  of  thermodynamics. 

2.  Is  the  system  capable  of  meeting  all  of  the  non-time  constraints? 

3 .  Can  the  system  be  realistically  developed  within  the  time  constraints  given?  If  new 
technology  is  involved  or  is  to  be  developed,  the  probability  of  completion  must  be 
assessed.  Is  the  risk  with  this  technology  acceptable,  or  does  a  backup  to  it  exist? 

In  order  for  a  candidate  system  to  advance  beyond  screening,  the  answers  to  all  of  these 
questions  must  be  yes. 

Question  1  should  be  answered  by  creating  a  basic  thermodynamic  model  of  the 
candidate  system.  Eqiupment  performance  may  be  assumed  to  have  some  reasonable  value,  or 
assumed  reversible.  (It  is  to  be  remembered  that  any  system  that  is  possible  should  not  be  failed 
here.)  Not  only  will  such  a  model  show  that  the  candidate  does  not  violate  the  laws  of  physics, 
but  may  also  serve  as  a  basis  for  optimization  during  preliminary  activities. 


(v)  Scenarios 

The  following  two  scenarios  will  illustrate  what  is  desired  as  a  result  of  the  feasibility 
study.  Furthermore,  both  of  these  scenarios  will  reappear  later  as  illustrations  of  points  developed 
later. 


1)  A  Light  Aircraft  Propulsion  System 

A  propulsion  system  is  needed  for  a  light  aircraft’s  airframe.  The  light  aircraft,  in  this 
case,  is  the  Glastar  kit-built  aircraft  (see  Figure  1  and  Appendix  A).  The  product,  fuel  and 
byproduct  analysis  is  relatively  straightforward:  the  product,  of  course,  is  thrust;  the  fuel  may  be 
automotive  fuel,  aviation  fuel  or  jet  fuel;  byproducts  include  noise  and  vibration,  as  well  as 
weight  and  drag. 

One  way  of  expressing  the  problem  statement  is  thus: 

Develop  a  propulsion  system  for  the  given  airframe  Cruise  should  be  in  the  140-knot  range;  good 
short-field  performance  should  be  achieved.  Endurance  should  be  in  the  3-6  hour  range.  As 
aviation  fuel,  automotive  fuel  and  jet  fuel  are  commonly  available  at  airports,  all  three  are 
acceptable  fuels. 

The  following  lists  of  solutions  is  generated: 

1.  Piston  engine  with  a  constant-speed  propeller 

2.  Piston  engine  with  a  fixed-pitch  propeller 

3.  Diesel  engine 

4.  Gas  turbine  engine  (turboprop) 

5.  Jet  engine  (turbojet  or  turbofan) 

Initial  screeiung  eliminates  number  5  immediately,  as  the  aiiframe  is  not  adaptable  to  a 
jet  engine.  Likewise,  numbers  3  and  4  will  be  eliminated  here,  as  there  currently  are  no  suitable 
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turboprop  installations  for  this  airframe.^  Thus,  options  1  and  2  will  advance  to  preliminary 
design. 


Figure  1:  Glastar®  Kit-built  Aircraft 

21  A  “Beer  and  Soup”  Plant^ 

What  follows  is  a  fictional  problem  statement  for  machinery  at  a  fictional  “beer  and 
soup”  plant: 

A  plant  is  producing  beer  and  soup.  In  the  final  processing,  it  is  necessary  to  hold  the  soup  at 
368K  and  the  beer  at  282K.  Typical  ambient  temperatures  are  between  these  two  values.  Natural 
gas  and  electricity  are  available  at  the  plant  sight  as  energy  sources.  166  kW  of  heat  energy  must 
be  supplied  to  the  soup  and  40  kW  removed  from  the  beer  in  the  process. 

One  can  easily  imagine  several  solutions.  The  beer  may  be  refrigerated,  while  the  soup  is 
heated  by  natural  gas-fired  burners.  The  soup  could  be  heated  by  electric  heating  coils.  The  soup 
could  be  heated  by  steam  generated  in  a  gas-fired  or  electric  boiler.  Or,  a  heat-pump/refrigeration 
unit,  as  shown  in  Figure  2,  could  be  employed.  Finally,  more  than  one  of  these  methods  could  be 
combined. 


^This  is  really  a  time  constraint,  as  a  turboprop  or  diesel  engine  could  certainly  be  developed  for  an  aircraft 
such  as  the  Glastar.  In  fact,  much  research  is  under  way  in  this  area. 

^  Although  here  described  as  a  “Beer  and  Soup”  plant,  such  a  system  could  also  be  employed  on  an  aircraft 
for  combined  cooling  and  heating  loads,  such  as  avionics  and  cabin  heat. 
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Figure  2:  Heat  Pump/Refrigeration  System 


(b)  Preliminary  Design 

The  end  result  of  the  preliminary  activities  should  be  the  selection  of  the  best 
configuration  from  all  of  the  candidate  systems.  In  order  to  accomplish  this,  an  objective 
function  must  be  formulated,  basic  optimization  must  be  performed  on  the  candidate  systems  and 
the  advantages  and  disadvantages  of  each  candidate  system  should  be  listed. 

(i)  The  Objective  Function 

It  is  not  possible  to  properly  evaluate  the  candidate  systems  without  a  function  that 
mathematically  quantifies  their  relative  merit.  For  many  energy  systems,  especially  stationary 
systems,  the  objective  function  is  the  lifecycle  cost,  and  it  is  to  be  minimized.  An  alternative  is  a 
profit  to  be  maximized. 

When  an  energy  system  is  used  in  a  vehicle,  either  as  the  prime  mover  or  for  other 
purposes,  such  as  climate  control,  the  objective  function  becomes  more  complicated.  This  is 
because  although  energy  drives  a  vehicle,  the  purpose  of  the  vehicle  is  not  simply  to  perform 
energy  conversion.  A  vehicle  exists  to  provide  performance,  and  this  cannot  be  expressed 
directly  in  energy  terms.  Tlierefore,  the  objective  function  for  a  vehicle  must  account  for  the 
benefits  of  performance.  If  an  energy  system  is  to  be  designed  for  use  on-board  a  vehicle,  the 
objectives  for  the  overall  vehicle  must  be  translated  into  terms  applicable  to  the  energy  system 
itself  Such  translation  involves  thermoeconomics  and  weighting  factors.  Additional  costs, 
outside  of  the  typical  fuel  and  capital  costs,  may  need  to  be  included.  One  example  of  an 
additional  cost,  highly  applicable  to  aircraft  design,  is  the  cost  of  weight. 
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Qualitative  objectives  for  Scenario  1,  presented  above,  are  listed  in  Table  4.  The 
formulation  of  objective  functions  is  discussed  in  detail  in  Section  4.  Additionally,  overall  and 
subsj'stem  objective  functions  are  developed  in  detail  during  Section  7,  in  a  continuation  of 
Scenario  1. 


Objective 

Extremization 

Cruise  speed  at  75%  power 

maximize 

Range  at  65%  power 

maximize 

Rate  0^^  climb  at  gross  weight 

maximize 

Take-off  roll 

minimize 

Useful  load 

maximize 

Cost 

minimize 

Table  4;  Qualitative  Objectives  for  Scenario  1. 


(ii)  Optimization  of  Candidate  Systems  and  the  Sensitivity  Analysis 

The  purpose  of  the  optimization  of  candidate  systems  is  to  compare  candidate  systems  on 
an  equal  basis.  Without  this  optimization,  one  may  be  comparing  a  poor  iteration  of  one 
candidate  to  a  good  iteration  of  another,  with  the  latter  appearing  to  be  tiie  better  candidate.  After 
optimization  of  both,  it  could  become  clear  that  the  initi^y  poorer  performing  system  is  indeed 
the  better  choice. 


The  purpose  of  the  sensitivity  analysis  is  to  evaluate  the  rate  of  change  of  the  objective 
function  to  each  decision  variable  of  each  candidate  system.  The  sensitivity  analysis  tells  the  , 
designer  how  critical  it  is  to  achieve  the  optimal  values  of  each  decision  variable.  If  a  decision 
variable  has  an  optimal  value  that,  for  whatever  reason,  may  be  difficult  to  achieve,  a  low 
sensitivity  to  this  value  is  naturally  desirable. 

During  tiie  course  of  optimization  and  the  sensitivity  analysis  it  will  be  necessary  to 
determine  which  design  variables  are  system-level  and  which  are  device-level.  This  information 
is  important  not  only  for  later  occurring  detail  activities,  but  also  here,  as  in  this  stage  of  the 
design  process  it  is  desirable  to  optimize  on  system-level  variables  only,  while,  in  this  stage, 
typical  values  should  be  used  for  device-level  variables.  Generic  models  should  be  used  for 
equipment  and  costing.  The  idea  is  to  provide  sufficient  information  for  candidate  screening; 
excess  effort  spent  on  the  candidates  not  chosen  is  wasted. 

Optimization  requires  an  optimization  algorithm,  a  simulation  algorithm  and  cost  models. 
These  topics  are  well  developed  io  the  literature.^ 

Mathematically  the  sensitivity  analysis  is  performed  by  taking  the  partial  derivative  of  the 
objective  function  with  respect  to  each  system-level  decision  variable. 


where  S',  is  the  sensitivity  coefficient  for  objective  function  J  and  decision  variable  _y,. 

The  sensitivity  analysis  and  the  optimizations  are  intricately  intertwined.  For  the 
sensitivity  analysis  to  be  accurate,  the  values  of  the  decision  variables  should  be  near  their 
optimiun  values.  If  it  is  not  clear  which  variables  are  system-level  and  which  are  device-level, 
the  sensitivity  analysis  may  need  to  be  extended  so  as  to  investigate  the  effects  of  decision 


See,  for  example  El-Sayed  and  Evans,  1970;  El-Sayed,  1996 
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variables  on  one  another.  Namely,  partial  derivatives  of  the  optimal  value  of  a  decision  variable 
with  respect  to  the  variable  in  question  may  be  taken. 


(2) 


These  derivatives  should  be  of  relatively  low  magnitude  for  the  variable  in  question  to  be 
considered  device-level. 


After  optimization  the  value  of  the  objective  function,  the  value  of  all  relevant  decision 
variables  and  the  sensitivity  of  liie  objective  function  to  these  should  be  listed  for  each  candidate 
system. 


(Hi)  Listing  of  Advantages  and  Disadvantages 

Ideally,  all  advantages  and  disadvantages  of  the  candidate  systems  will  be  mathematically 
included  in  the  objective  function.  However,  in  reality  this  will  not  always  be  the  case.  Many 
desirable  traits,  such  as  maintainability  and  ease  of  replacement  part  procurement  may  be  difficult 
to  quantify.  Any  of  these  traits  should  be  listed  along  with  the  values  of  the  objective  fimetions. 


To  continue  with  the  example  presented  by  scenario  1,  Table  3.5  shows  some  possible 
piston  engine  combinations  along  with  their  advantages  and  disadvantages. 


Engine-Propeller 

Advantages 

Disadvantages 

Lycoming  0-320  or  0-360 
engine^,  fixed-pitch  propeller 

Simplicity,  less  maintenance, 
engine  parts  are  readily 
available,  aircraft  mechanics 
are  very  familiar  with  engine. 
Lightweight  propeller. 

Propeller  will  always 
compromise  speed,  climb  or 
both. 

Lycoming  0-320  or  0-360 
engine,  constant-speed 
propeller 

Engine  parts  are  readily 
available,  aircraft  mechanics 
are  very  familiar  with  engine. 
Very  good  all-round 
performance.  Lower  cruise 
fuel  consumption. 

Constant  speed  propeller  will 
add  complexity,  weight  and 
maintenance  costs. 

Walther-LOM  160 
horsepower  engine®,  constant 
speed  propeller 

i 

Overhead-camshaft  design  is 
more  efficient.  Fuel  injection 
eliminates  manual  leaning  of 
mixture.  Supercharger 
provides  better  power  at 
altitude  than  the  Lycoming  O- 
320  does.  Lightweight  engine. 

Constant  speed  propeller  will 
add  complexity  and 
maintenance  costs.  LOM 
engine  not  nearly  as  common 
as  Lycoming.  The  engine’s 
configuration  requires  a  longer 
cowling.  Engine  is  only  rated 
at  140  horsepower  continuous. 

Table  5:  Engine-propeller  combination  advantages  and  disadvantages 


*  The  O  here  stands  for  opposed,  and  the  number  for  the  engine  displacement.  These  engines  are 
carbureted,  produce  from  150-180  horsepower,  are  built  in  the  United  States  and  are  very  standard  in  light 
aircraft. 

®  The  LOM  engine  is  a  fuel  injected  and  supercharged  inline  engine. 
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(iv)  Choosing  the  Best  Candidate 

The  best  candidate  will  nonnally  have  the  highest  value  of  the  objective  function. 
However,  before  simply  choosing  the  candidate  with  the  highest  value,  the  following  questions 
should  be  answered.  (This  list  should  not  be  construed  as  all-inclusive.) 

1 .  Do  any  of  the  decision  variables  have  values,  after  optimization,  that  may  be  difficult  to 
achieve  due  to  technological  difficulties  or  constraints.  Until  detail  design  is  complete,  there 
is  always  some  doubt  as  to  whether  the  system  can  be  built  as  envisioned  during  this  phase. 
The  effect  on  the  value  of  the  objective  toction,  whi^h  occurs  due  to  movement  from  the 
optimal  value,  may  be  estimated  from  the  results  of  the  sensitivity  analysis  with 

AJ^S^Ay,  (3) 

With  this,  the  effects  of  non-achievement  of  the  results  of  the  optimization  may  be  estimated, 
and  the  risk  associated  with  developing  a  candidate  may  be  evaluated. 

2.  Does  the  candidate  with  the  highest  value  of  the  objective  function  have  any  major  (non- 
quantified)  disadvantages  to  the  other  candidates? 

3.  Is  there  uncertainty  in  meeting  time  constraints? 

Yes  answers  to  any  questions  like  these  may  be  valid  reasons  for  choosing  a  candidate 
other  than  that  with  the  highest  value  of  the  objective  function.  Nonetheless,  with  a  good 
objective  function,  a  cost  may  be  put  on  a  choice.  For  example,  one  might  say  choosing  design  X 
over  Y  will  cost  so  much,  but  will  yield  a  greater  level  of  maintainability.  Likewise,  risk-benefit 
analysis  may  be  accomplished. 

If  multiple  candidates  have  similar  niunerical  values  for  the  overall  objective,  additional 
effort  is  necessary  to  make  the  selection.  It  may  be  necessary  to  perform  a  more  detailed 
optimization  to  confidently  choose  one  candidate  system.  Also,  for  similar  values  in  the  objective 
fhnctions,  generally  the  system  with  smaller  sensitivities  will  be  better. 

At  the  end  of  this  stage,  the  designer  will  have  picked  a  single  candidate  for  advancement 
to  detailed  design. 

(c)  Detailed  Design 

Detailed  design  begins  with  (or  in  a  sense  is  preceded  by)  “generic  optimization”.  The 
purpose  of  generic  optimization  is  to  yield  a  design  that  serves  as  a  starting  point  for  refinements 
(to  the  decision  variables)  to  be  carried  out  during  the  detailed  design.  One  reason  the  generic 
initial  optimization  is  needed  is  because  the  detailed  design  and  selection  of  different  systems, 
subsystems  and  devices  will  be  carried  out  by  “independent”  design  teams. 

(i)  Generic  Optimization 

This  step  is  performed  by  the  “overall-system  design  team”’. 

The  selected  candidate  design  is  “optimized”  generically.  Here  the  optimization  should 
include  device-level  variables.  Ideally,  this  would  be  accomplished  using  software,  which 
includes  cost  and  performance  models. 


’  The  different  teams,  as  referred  to  here,  may  in  fact  consist  of  only  a  single  individual  or  there  may  be 
overlap  in  personnel.  A  team  may  even  be  from  a  different  firm. 
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In  the  absence  of  such  software,  the  initial,  “optimized”  design  could  be  a  design  penned 
by  an  engineer  with  sufficient  “fingertip-feel”  to  produce  a  good  design.  This  “optimization” 
could  simply  be  the  existing  design  of  a  current,  state-of-the-art  vehicle,  a  perturbation  of  such  an 
existing  design  or  a  composite  of  several  existing  designs. 

In  this  step,  exergy  costs,  expressed  in  the  units  of  the  objective  fimction,  associated  with 
mass  or  power  streams  need  to  be  determined.  Also,  each  device  in  the  design  needs  to  have  its 
“optimized”  parameters  listed.  Examples  of  these  parameters  are  capital  costs,  efficiencies  and 
weights. 

(ii)  Subsystem  and/or  Device  Design  and/or  Selection 

The  optimized  parameters  from  generic  optimization  (i)  are  given  to  the  “device  design 
teams”  as  initial  design  values,  along  wife  the  relevant  exergy  costs  and  an  objective  function  for 
their  device,  which  includes  only  local  variables.  System-level  decision  variables  are  to  be 
considered  as  constraints.  Device  design  teams  are  responsible  for  the  detailed  design  of  a  device 
or  the  selection  of  a  device.  The  designers  are  encouraged  to  better  the  results  (e.g.  cost,  weight, 
efficiency  of  their  device)  from  the  initial  system  design  and  optimization,  striving  for  the 
optimum  of  die  local  objective  fimction.  (The  method  for  developing  the  local  objective  fimction 
will  be  presented  shortly.)  Thus,  the  generic  design’s  optimized  device-level  decision  variables 
serve  as  a  point  of  departure  for  the  detailed  design. 

Several  different  scenarios  are  foreseeable  for  a  design.  For  examples: 

•  The  device  will  be  designed  in-house. 

•  The  device,  for  which  a  design  team  is  responsible,  is  to  be  purchased  from  an  outside  source 
as  an  “off-the  shelf’  component.  One  example;  An  architect-engineer  firm  is  designing  a 
combined-cycle  power  plant.  The  team  responsible  for  the  gas  turbine  engine  will  be 
selecting  an  existing  model  from  a  firm,  which  specialized  in  building  such  engines. 

•  The  device  is  to  designed  by  an  outside  firm.  Here  the  outside  firm  becomes  the  device 
design  team,  and  should  be  supplied  with  the  same  information  that  an  in-house  design  team 
would  receive.  Heat  exchangers  are  an  example  of  a  device,  which  often  fit  in  this  category. 

(Hi)  Perturbations 

There  may  be  instances  where  a  device  design  team  may  need  or  propose  to  perturb 
system-level  variables.  Also,  in  the  search  for  improvements,  there  may  be  times  when  it  is 
desirable  to  do  so.  For  example,  it  may  propose  to  drop  an  operating  pressure  somewhat,  because 
a  drastic  reduction  in  weight  (or  material  cost)  could  be  achieved,  due  to  materials  being  often 
available  only  in  discrete  sizes.  When  this  happens,  the  device  team  must  consult  with  the  system 
team.  Decisions  can  be  facilitated  by  reruiming  the  simulation  or  using  the  results  of  the 
sensitivity  analysis. 


(iv)  Simulation  and  Proofing 

After  all  initial  device  designs  are  completed  the  simulation  is  rerun  with  models  of  the 
proposed  equipment  to  ensure  proper  on  and  off-design  operation.  At  this  point,  it  is  desirable  to 
globally  optimize  once  again  any  operating  variables.  If  these  do  not  vary  from  the  values 
resulting  from  the  iiutial,  generic  optimization,  the  thermal  design  is,  for  all  practical  purposes, 
finished.  It  is  even  conceivable  that  fiulher  iterations  of  the  detailed  design,  with  updated  exergy 
costs,  would  diverge  further  from  the  tiieoretical,  optimal  design.  This  re-optimization  puts  such 
variables  in  line  with  the  actual  devices  used,  as  opposed  to  being  optimized  for  generic,  non¬ 
existent  equipment. 
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If  in  the  re-optimization  of  operational  variables,  noticeable  excursions  from  the  initial 
design  values  occur,  it  may  be  desirable  to  reconsider  the  global  variables  and  repeat  steps  2 
through  4. 

fvj  Updating  of  Generic  Models 

Device  design  will  be  documented  so  as  to  improve  and  update  weight  and  cost  models. 
Now,  methods  for  developing  the  objective  fimctions  (overall,  and  device)  referred  to  above  will 
be  proposed. 

V.  DEVELOPING  AN  OVERALL  OBJECTIVE  FUNCTION 

(a)  A  Proposed  Overall  Objective  Function 

Before  beginning  any  design,  an  engineer  (or  team)  must  determine  the  constraints  and 
objectives  of  the  project.  For  a  meaiungful  optimization  the  objectives  need  to  be  expressed  in  a 
single  function. 

A  traditional  objective  function  for  an  overall  energy  system  with  a  single  feed  and  single 
product  is  (El-Sayed,  1995) 

J  =  CfF  +  Y,Z-v^P  (4) 

where  J  is  to  be  minimized  (-J  is  the  profit),  F  symbolizes  the  amount  of  feed,  Z  symbolizes 
capital  equipment  and  P  is  the  product.  The  lower  case  c’s  and  v’s  represent  value,  per  unit,  of 
the  feed  and  product,  respectively. 

Although  a  vehicle  has  energy  feeds,  its  product  is  performance  rather  than  energy. 
Nevertheless,  equation  (4)  remains  appropriate  for  the  overall  vehicle,  when 

J  =  (5) 

Each  P  is  now  used  to  signify  the  various  performance  benefits.  For  the  case  of  an 
aircraft.  P’s  might  accoimt  for  climb  rate  and  cruising  speed,  for  examples. 

Taking  equation  (5)  and  negating  it  yields  a  “profif’  to  be  maximized,  which  will  be 
referred  to  as  n. 

n-2;v/-[c,F+Xz]  «) 

This  objective  function  is  useful  in  conceptual  and  preliminary  design,  but  cannot  directly 
be  applied  to  systems  or  subsystems,  without  simulation  of  the  entire  vehicle.  The  system  or 
subsystem  will  have  inputs  and  outputs  that  are  mass,  momentum  or  entropy  streams.  In  order  to 
evaluate  such  a  system  or  subsystem,  (6)  must  be  translated  into  a  function  that  contains  terms  for 
the  exergy  flows  are  being  produced/consumed.  This  process  will  be  demonstrated  later. 

(b)  Determining  the  Value  of  Performance 

Clearly,  before  one  can  ^ply  (6)  to  a  conceptual  design,  one  must  find  the  numerical 
values  of  the  various  Vp’s.  These  values  are,  in  tiie  most  general  case,  functions  of  the  levels  of 
performance  themselves.  One  would  expect  a  very  low  level  of  performance  for  a  given 
desideratum  would  be  zero,  that  is,  it  is  unacceptable. 
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As  an  example  of  this,  consider  an  automobile.  A  person  residing  in  the  United  States 
(with  typical  speed  limits  in  the  100-125  km/hr  range)  might  be  unwilling  to  consider  purchasing 
a  vehicle  incapable  of  at  least  125  km/hr. 

However,  this  person  would  likely  be  willing  to  pay  more  for  an  automobile  that  could 
achieve  a  higher  speed  than  this  minimal  amoimt.  Speed  limits  (or  more  correctly,  the 
enforcement  thereof)  cause  a  reduction  in  the  marginal  value  of  a  1  km/hr  increase  in  speed  as  the 
top  speed  of  the  car  becomes  higher  and  higher.  This  hypothetical  person,  therefore,  might  not  be 
willins  to  pay  any  more  money  for  an  automobile  capable  of  reaching  240  km/hr  than  an 
automobile  capable  of  reaching  200  km/hr.  (Acceleration  would  be  a  separate  desideratum.) 

Returning  to  the  example  of  the  fighter  plane,  one  can  imagine  that  a  certain  level  of 
overall  performance  could  be  reached  that  would  insine  victory  over  any  opponent  and 
penetration  of  any  airspace.  It  is  not  rational  to  invest  further  resources  in  performance  beyond 
this  point.  So  we  expect  that  the  fimctions  Vp(P)  become  flat  after  a  certain  value  for  a 
desideratum  is  reached.  In  between  these  points  (the  minimum  acceptable  and  the  maximum 
useful)  lies  some  continuous  function. 

A  method  is  laid  out  below  to  estimate  these  functions.  It  consists  of  five  steps,  and  may 
involve  iteration.  It  should  be  noted  that  in  order,  ultimately,  to  optimize  a  subsystem,  it  is 
absolutely  necessary  to  develop  this  information,  that  is,  the  functions  Vp(PJ,  in  some  form,  even 
if  it  is  not  done  with  the  method  used  in  this  dissertation. 

(i)  Step  1 :  Determining  Median  Performance 

The  foregoing  algebraic  “tradeoff  functions”  s>p(P)  for  representing  the  value  of 
performance  may  have  an  arbitrary  shape  between  the  points  of  minimum  acceptable  and 
maximum  desirable  performance.  However,  if  a  limited  range  of  performance  is  considered,  the 
assumption  of  a  linear  relationship  is  reasonable.  One  way  in  which  this  linear  function  may  be 
constructed  is  around  a  median  point.  If  a  linear  function  is  undesirable,  this  information  will  still 
be  of  use  in  the  construction  of  the  function. 

For  a  military  combat  aircraft,  one  way  to  find  a  median  point  is  by  considering  the 
performance  of  the  aircraft’s  adversaries,  both  current  and  projected*.  For  each  desideratum,  at 
least  one  of  tiie  aircraft  in  the  adversarial  group  has  a  best  value.  That  value  could  be  selected  as 
the  median  of  acceptable  performance  levels.  The  set  of  medians  would  form  a  “standard”  of 
comparison  for  further  investigations. 

The  same  basic  idea  is  applicable  to  civilian  vehicles.  However,  the  median  values  of  a 
market  segment  might  instead  be  chosen  to  set  the  “standard  values.” 

(ii)  Step  2:  Projected  Units  Costs  and  Projected  Production 

A  realistic  estimate  of  both  the  cost  per  vehicle  and  the  total  production  quantity  of  the 
vehicle  should  be  made.  As  will  be  seen  below,  this  step  is  an  aid  to  the  fourth  step. 


*  Here,  the  aircraft’s  adversaries  are  considered  to  be  other  aircraft.  This  view  could  of  course  be  expanded 
to  consider  enemy  air  defense  capabilities. 
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(Hi)  Step  3:  Projected  Research,  Design  and  Development  costs 

The  projected  research,  design  and  development  costs  should  be  listed.  This  is  also  an 
aid  to  the  following  step.  When  employed  with  the  information  attained  above,  the  total  project 
cost  may  be  estimated. 


(iv)  Step  4:  Algebraic  Tradeoff  Functions 

The  algebraic  tradeoff  functions  for  each  desideratum  must  be  determined  through 
questioning  of  the  customer  o’'  '^nd-user.  As  a  bare  minimum,  the  following  three  questions 
should  be  asked  (in  some  form) 

1 .  What  is  the  minimum  acceptable  value  for  each  of  the  desiderata? 

2.  Is  there  a  point  beyond  which  further  improvement  is  not  necessary? 

3.  How  much  would  a  given  improvement,  over  and  above  the  median  value  determined  in 
Step  1,  be  valued? 

With  the  answers  to  these  three  questions  determined,  the  simplest  tradeoff  function, 
linear,  may  be  determined.  Once  again  visiting  the  automobile  example,  let  us  imagine  that  the 
marketing  department  has  asked  the  above  three  questions  to  potential  customers  regarding  the 
top  speed  of  an  automobile  of  a  certain  class.  The  average  (or  weighted  average)  answers  were, 
5=125  km/hr  as  a  minimiim  acceptable  top  speed,  5=200  km/hr  as  a  ceiling  beyond  which 
improvement  has  little  or  no  value,  and  a  willingness  to  pay  800  dollars  for  an  improvement 
A5=10  km/hr  over  a  median  165  km/hr  top  speed.  One  can  imagine  a  function,  v^S),  which 
would  have  a  value  of  zero  dollars  up  to  125  km/hr  and  rise  with  a  slope  of  80  dollars/(km/hr)  to  a 
maximum  of 6000  dollars. 

However,  for  many  cases  the  questions  may  not  be  best  asked  in  terms  of  dollars.  For  a 
military  vehicle,  say  an  air  superiority  fighter,  dollars  would  be  a  poor  choice  of  units.  In  this 
example,  the  end-user  (the  Department  of  the  Air  Force  or  Navy)  is  not  the  same  as  the  purchaser 
(Congress).  The  pinchase  costs  have  reached  such  high  amounts  that  it  is  difficult  for  the  average 
person  to  comprehend  the  sums  in  rational  terms.  Furthermore,  neither  body  (Congress  or  the 
end-user)  is  spending  their  own  money. 

In  such  a  case  the  questions  may  be  rephrased  in  terms  of  production  sacrifices;  Question 
Number  3  could  be  changed  to:  “What  reduction,  in  number  of  aircraft  delivered  to  you,  would 
you  accept  in  order  to  obtaiu  a  specified  improvement  firom  the  median  value  determined  in  Step 
1?”  A  miUtary  leader  should  have  a  good  grasp  of  tradeoffs  between  quantity  and  quality.  The 
information  firom  Steps  2  and  3  allows  production  tradeoffs  to  be  converted  to  a  dollar  amount  (or 
a  cost  to  a  production  adjustment). 

If  a  computer  simulation  were  to  be  available  that  would  predict  aircraft  survivability  as  a 
function  of  measured  desiderata,  it  could  be  used  to  develop,  or  help  develop,  these  trade-off 
functions. 


(v)  Step  5:  Relative  Weighting 

At  the  overall  vehicle  level,  the  algebraic  tradeoff  functions  are  sufficient  for 
optimization.  As  shown  below,  in  order  to  decompose  the  vehicle  into  systems  and  subsystems, 
weighting  factors  will  be  needed.  Then  a  distinct  objective  function  can  be  defined  for  each 
system,  subsystem  or  device,  each  to  be  designed  by  a  distinct  “team”. 

The  relative  weighting  factor  for  the  desideratum  ,  as  employed  in  this  dissertation,  is 
deduced  fi'om  the  foregoing  information  and  is  given  by 
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Where  $  represents  total  money,  P  a  performance  desideratum,  and  AS  the  increase  in 
price  an  end  user  would  pay  for  a  AP  increase  in  a  desideratum.  The  numerator,  then,  represents 
the  percent  of  expenditure  increase  the  customer  is  willing  to  make,  per  percent  increase  in  a 
given  performance  desideratum  Pi.  The  denominator  is  the  sum  over  all  desiderata;  so  Wt 
represents  the  relative  importance  of  Pi.  $  cancels,  leaving 


(vi)  Iteration 

It  may  be  necessary  to  repeat  Steps  1  through  5  as  a  preliminary  design  is  completed. 

This  is  due  to  an  implicit  assumption  as  to  the  independence  of  the  individual  desiderata.  An  end- 
user  may  sacrifice  far  more  total  resources  than  intended,  as  this  person  was  looking  at  only  a 
single  desideratum  at  a  time  when  proceeding  through  the  five  steps.  Or,  the  standard 
performance  may  be  so  great  that  the  user  erred  in  the  opposite  direction.  Therefore,  there  must 
be  good  commimication  between  designer  and  end-user  at  all  times.  As  a  conceptual  design 
evolves,  the  questions  fi'om  Step  3  may  need  to  be  repeated.  A  design  may  have  become,  in  the  . 
end-users  opinion,  too  expensive;  or  its  performance  may  simply  be  inadequate.  This  is 
especially  true  if  the  median  performance  was  far  fi-om  the  optimum. 


(vii)  Notes  on  the  Development  of  an  Overall  Objective  Function 

In  the  methods  presented  here,  the  units  of  each  term  in  the  objective  functions  are 
monetary.  For  commercial  vehicles  and  transport  vehicles,  monetary  units  are  the  obvious 
choice.  For  combat  vehicles,  other  units  may  be  bettev.-,  such  as  “production  quantity”. 

Vi.  THERMOECONOMICS  AND  DETAILED  DESIGN 

(a)  Concurrent  Engineering  and  the  Need  for  Decomposition 

Once  a  conceptual  design  has  been  accepted,  the  detailed  design  should  proceed  in  an 
efficient  fashion.  It  is  not  possible  to  optimize  and  design  an  aircraft,  or  even  an  energy 
subsystem,  as  an  entirety.  Vehicular  energy  systems  are  too  complex  to  be  designed  by  a  single 
individual  or  team,  as  the  number  of  decision  variables  becomes  unmanageable.  Therefore  the 
design  or  selection  of  individual  devices  and/or  subsystems  is  delegated  to  subordinate  teams  or 
individuals.  It  is  desirable  for  these  designers  to  have  a  methodology  and  information  that 
allows  them  to  make  decisions  in  accordance  with  the  overall  goals  of  the  vehicle.  Furthermore, 
it  must  be  determined  which  decision  variables  are  local  and  which  are  global.  That  is,  which 
variables  a  device  design  team  is  free  to  vary  in  optimizing  its  individual  device,  and  which  are 
constraints  to  the  design  team  so  as  not  to  affect  the  designs  being  carried  out,  simultaneously,  by 
other  design  teams.  These  determinations  may  at  times  be  made  through  “common  sense”,  but  at 
other  times  may  require  a  sensitivity  analysis. 
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A  vehicle  typically  relies  on  one  fuel  to  achieve  its  goals.  Because  the  subsystem  or 
device  design  team  is  not  optimizing  a  whole  vehicle,  but  something  which  may  produce  or 
consume  commodities  not  considered  when  looking  at  the  whole  vehicle,  its  objective  function 
will  vary  from  that  of  the  overall  vehicle.  For  example,  ttie  alternator  on  a  light  general  aviation 
aircraft  does  not  directly  consume  fuel,  nor  does  it  directly  influence  performance.  Nonetheless, 
an  aircraft  with  an  alternator  that  is  both  lighter  and  more  efficient,  will,  with  all  else  remaining 
equal,  perform  better  (and  cost  more).  The  person  selecting  an  alternator  should  not  be  burdened 
with,  nor  at  this  stage  of  the  design  be  necessarily  capable  of,  directly  calculating  the  impact  upon 
aircraft  performance 

The  alternator  mentioned  above  still  has  only  one  product  (electrical  energy).  An  onboard 
energy  device  or  subsystem  may  use  several  “fuels”  and/or  supply  several  “products”.  These 
fuels  and  products,  besides  electricity  and  shaft  power,  include,  but  are  not  limited  to,  compressed 
air,  hydraulic  power  and  heat  (or  cooling).  Additionally,  one  must  accoimt  for  the  lift  required  to 
hold  the  device  in  the  air  (or  to  make  it  go  up).  The  costs  and  benefits  of  these  mass  and  energy 
flows  must  be  included  in  the  device  or  subsystem  design  team’s  objective  function. 

Before  proceeding  in  developing  an  objective  function,  for  a  device  or  subsystem,  which 
is  to  be  designed  by  a  team  subordinate  to  the  overall-vehicle  design  team,  methods  are  needed 
for  determining  the  costs  and  benefits.  This  is  done  through  accounting  for  mass  and  energy 
streams  with  exergy,  and  using  the  methods  of  thermoeconomics  to  place  unit  costs  on  them.^ 

(b)  The  Role  for  Exergy  and  Thermoeconomics 

The  role  for  exergy  and  thermoeconomics  discussed  within  this  paper  is  for  employment 
during  detailed  design.  In  this  paper  detailed  design  is  the  “filling  in”  of  black  boxes  such  as  a 
propulsion  system,  environmental  control  system  or  a  hydraulic  pump.  The  methodology 
presented  here  assumes  that  the  overall-vehicle  designers  have  penned  a  “reasonably  optimized” 
preliminary  design.  Additionally,  it  is  assumed  that  allowance  has  been  made  for  the  weight  and 
exergy  consumption  of  systems  and  subsystems.  Before  employing  this  thermoeconomic 
methodology,  it  is  absolutely  necessary  to  have  the  previously  defined  overall  objective  function, 
weighting  factors  and  values  of  performance.  Thus  all  design  teams  will  strive  for  the  same 
common  goals. 

(i)  Exergy 

For  comparing  energy  converters  that  have  single  feed  and  product  streams,  and  perform 
the  same  function,  performance  measures  such  as  “thermal  efficiency”  and  “COP”  often  suffice. 
Such  energy-based  performance  measurements  are  referred  to  as  first  law  performance  measiues. 
As  devices  or  systems  add  additional  feeds  and/or  products,  these  definitions  become  inadequate. 
Moreover,  even  if  a  device  has  a  single  feed  and  a  single  product,  but  is  one  component  in  a 
complex  system,  such  first  law  measures  are  inadequate  for  comparing  the  relative  impact  of  any 
one  device  upon  the  overall  system  performance.  That  is,  first  law  measures  of  feeds  and 
products  do  not  represent  the  true  values  of  these  commodities.  What  is  needed  in  analyzing  and 
optimizing  complex  systems  is  a  common,  consistent  measure  for  representing  the  potential  of 
mass  and  energy  flows. 

This  proper  method  of  quantifying  fuels  and  products  is  with  exergy  (sometimes  called 
“availability”  or  “available  energy”).  Exergy  measures  the  ability  of  a  commodity  (matter, 
momentum,  entropy  or  “heat”,  charge...)  to  cause  change. 


®  As  a  general  reference  to  exergy  and  thermoeconomics  see,  for  example,  Bejan  et  al,,  1996. 
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After  preliminary  design,  all  exergy  flows  and  destructions  within  the  proposed  vehicle 
should  be  diagrammed.  As  a  vehicle  may  be  evaluated  under  several  different  operating  modes 
(typically  related  to  the  desiderata),  it  may  be  necessary  to  create  several  exergy  flow  diagrams. 

If  one  were  developing  a  fighter  aircraft,  diagrams  may  be  necessary  for  cruise,  maneuver,  climb, 
etc. 

Exergy  is  dependent  on  the  environment  in  which  the  vehicle  operates.  If  the  range  of 
operating  conditions  is  not  too  large,  an  average  (or  weighted  average)  of  the  operating  conditions 
may  suffice  to  serve  as  the  reference  environment.  For  vehicles  operating  in  a  wide  variety  of 
conditions  (i.e.,  aircraft)  a  number  of  diagrams,  sufficient  to  cover  the  variations  in  operating 
environment,  will  be  necessary.  It  may  even  be  necessary  to  express  values  as  functions  of  time 
for  such  realms  of  operation  as  a  maximum  performance  climb,  where  ambient  pressure  and 
temperature  are  changing  rapidly. 

The  purpose  of  these  exergy  diagrams  is  twofold.*®  The  first  benefit  is  that  this  diagram 
allows  the  designers  to  truly  comprehend  the  interactions  and  inefficiencies  of  an  energy  system. 
The  relative  magnitude  of  Ae  destructions  (and  wastes)  will  point  to  the  areas  of  a  system  where 
improvements  will  be  most  beneficial.  This  topic  will  not  be  addressed  directly  here,  but  is 
addressed  extensively  in  the  literature.  The  second  benefit  of  creating  the  exergy  flow  and 
destruction  diagram  is  that  costs  may  be  assigned  to  each  stream’s  exergy  flow.  That  is,  the  cost 
per  imit  of  exergy  can  be  evaluated  at  each  station  on  the  exergy  flow  diagram.  Thus,  not  only 
exergy,  but  also  costs  can  be  tracked  through  the  overall  vehicle.  Moreover,  the  unit  costs  of 
exergy  entering  a  subsystem  (or  device)  are  of  special  relevance  to  that  subsystem’s  design  team. 
Before  proceeding  to  die  construction  of  an  exergy  flow  diagram  for  an  aircraft,  the  methods  for 
establishing  the  unit  costs,  thermoeconomics,  will  now  be  discussed. 


(ii)  Thermoeconomics 

Thermoeconomics  allows  a  cost  to  be  associated  with  each  exergy  flow.  These  unit  costs 
are  absolutely  necessary  in  the  development  of  the  objective  functions  for  an  on-board  device  or 
subsystem. 

After  the  overall  exergy  flow  diagrams  are  finished,  complete  exergy  costing  and  cost 
flow  diagrams  should  be  created  for  each  exergy  flow  diagram.  Capital  costs  should  not  be 
included  in  the  costing;  as  the  local  optimization  should  be  changing  design  values  only  by 
relatively  small  amounts,  the  proper  monetary  charges  are  for  marginal  costs.  (This  is  explained 
in  more  detail  in  Section  3.06.) 

For  a  subsystem,  the  costing  will  be  performed  one  of  two  ways,  depending  on  whether  it 
is  at  a  maximum  performance  condition  or  a  “cruise”  condition.  These  correspond  to  costing  for 
fixed  resources  or  fixed  output,  respectively  (El-Sayed,  1995).  In  a  maximum  performance 
condition,  the  vehicle’s  primary  energy  conversion  system  (i.e.,  its  propulsion  system)  is 
producing  as  much  power  as  possible.  Any  exergy  drawn  to  fuel  the  subsystem  will  be  detracting 
from  the  vehicle’s  performance.  Therefore  all  exergy  flows  must  be  assigned  a  per  unit  cost 
based  on  the  cost  to  vehicle  performance.  For  example,  a  light  aircraft  in  full  throttle  cruise  at  a 
specified  altitude,  where  the  desired  performance  is  a  cruise  speed,  can  have  a  marginal  value  of 
timist  exergy  assigned  with 


dv  an 

dXj.  dV 


(9) 


*®  See,  for  example  Moran,  1989  and  Szargut  et  al.,  1988. 
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Here,  — r—  is  the  partial  derivative  of  velocity  with  respect  to  thrust  exergy  and  — 

dXj-  dV 

of  the  cruise  speed  “tradeoff’  curve.  Multiplying  these  two  partial  derivatives  yields 


which  clearly  is  the  marginal  profit  associated  with  thrust  exergy. 

This  information  allows  the  complete  exergy  costing  to  be  completed. 


is  the  slope 

an 

dXr  ’ 


On  the  other  hand,  for  the  case  of  “part-throttle”  cruise,  costing  is  based  on  the  cost  of 
fuel  Engine  power  is  no  longer  the  limiting  commodity;  any  drop  in  performance  may  be 
accommodated  by  increasing  the  power  output  of  the  engine.  One  resulting  cost  is  the  direct  cost 
for  additional  fuel.  Furthermore,  if  range  is  a  desideratum,  the  cost  of  an  increase  of  fuel  bum  to 
range  may  be  added  to  the  fuel’s  base  cost.  This  will  in  general  be 


range 


dR  an 

drhj^,  dR 


(10) 


If  a  device  is  a  member  of  the  vehicle’s  propulsion  system  (e.g.,  a  piston  engine  in  a 
propeller  driven  aircraft),  it  will  be  shown  that  exergy  costing  will  have  to  be  done  from  both  of 
the  preceding  viewpoints. 

The  application  of  the  second  law  of  thermodynamics,  combined  with  economics,  allows 
for  decomposition  of  systems  into  devices  and  subsystems  that  may  be  independently  optimized. 
This  theory  has  been  developed  and  validated  by  El-Sayed  and  coworkers.  (See  for  example,  El- 
Sayed  and  Evans,  1970;  El-Sayed,  1995.) 


(c)  The  Objective  Function  for  an  on-board  Subsystem 

With  the  exergy  flow  and  cost  diagrams  complete,  there  is  enough  information  to 
optimize  a  subsystem.  Here  a  subsystem  is  defined  as  an  on-board  energy  conversion  system  not 
directly  involved  in  producing  the  necessary  basic  forces  of flight  -  thrast  and  lift.  The 
subsystem  is,  in  general,  required  to  deliver  a  certain  output  as  a  design  constraint.  Then,  the  local 
objective  function  for  a  subsystem,  to  be  minimized,  is  a  total  cost. 


Uj 


(11) 


Here,  i  represents  each  desideratum,  7  each  feed  stream,  Cy  are  unit  costs  of  exergy,  c„  is  the  unit 
cost  of  mass,  and  m  is  the  mass  of  the  component.  For  an  aircraft,  the  cost  of  mass  is  related  to 
the  exergy  of  lift.  As  was  stated  previously,  this  exergy  expenditure  used  to  hold  mass  in  die  air 
must  be  taken  into  account.  Suitable  relations  will  be  given  in  the  example  that  follows.  W 
represents  the  various  weighting  factors,  determined  according  to  the  procedure  described  in 
Section  3.04. 

In  general,  costs  will  be  functions  of  the  decision  variables  themselves.  For  example,  the 
marginal  cost  of  weight  on-board  an  aircraft  is  a  function  of  the  aircraft’s  gross  weight. 

Therefore,  as  a  subsystem  designer  strays  from  die  estimates  made  about  the  subsystem  during 
preliminary  design,  iterations  may  be  necessary,  so  as  to  update  exergy  and  weight  costs. 


(d)  The  Objective  Function  for  a  Propulsion  System  Component 

While  a  subsystem,  in  general,  will  have  its  required  output  a  design  constraint,  the 
propulsion  system  itself  is  a  balance  between  performance  supplied  to  the  aircraft  (as  thrust). 
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weight,  capital  costs,  range  and  operating  costs.  Its  objective  fimction  differs  from  that  of  a 
subsystem  in  that  there  is  now  a  profit  to  be  maximized  -  the  difference  between  the  performance 
delivered  and  the  cost  of  providing  it.  This  takes  the  form  of  a  “profit”  to  be  maximized: 


n=2;»', 


m 


\>.J 


i,k 


i  a 


(12) 


It  is  necessary  to  cost  exergy  in  both  manners.  That  is  to  say,  costing  will  be  performed 
once  as  if  the  input  (fuel)  were  fixed,  and  once  as  if  the  output  (performance)  were  fixed.  Exergy 
outputs,  with  subscript  j,  are  determined  by  placing  a  value  on  thrust.  Inputs,  with  subscript  k,  are 
determined  using  the  price  of  fuel  (including,  if  necessary,  cost  adjustment  for  range).  The 
variable  T  is  a  time  factor,  the  percentage  of  time  a  vehicle  spends  in  each  operating  mode.  This 
follows  from  the  overall  vehicle’s  objective  function;  we  charge  a  device  in  the  propulsion  system 
for  only  the  fuel  consumption  for  which  it  is  responsible.  In  contrast,  performance  is  available  all 
of  the  time;  therefore  the  output  values  are  multiplied  by  the  weighting  factors  of  the  various 
desideratum. 


VII.  NOTES  ON  EXERGY  FOR  THE  DESIGN  PROCESS 


(a)  Exergy  of  lift  and  exergy  flow  in  aircraft 


(i)  Exergy  of  Lift 

Although  the  exergy  produced  by  the  wing  in  level  cruise  (at  the  expense  of  drag,  i.e., 
through  an  input  of  exergy  with  momentum)  is  ultimately  destroyed,  it  is  necessary  to  calculate 
the  value  of  the  lift  exergy  so  as  to  properly  charge  individual  components  for  their  lift 
requirements. 

In  general,  for  any  force,  the  exergy  transfer  associated  with  it  may  be  defined  as 
^/«.=F(V-V„)  (13) 


During  steady  state  flight,  the  force  of  lift  is  equal  to  the  weight  of  the  airplane  (or  component 
being  lifted).  The  problem  remaining  is  to  properly  define  V# .  An  ideal  wing,  producing  no 

entropy,  travelling  forward  at  speed  Vx  will  fall  towards  the  earth  at  a  certain  rate,  Vyo.  This  is  a 
thermostatic  state  (Gaggioli  et.  al.,  1999),  and  will  be  used  as  a  dead  state  here. 

McCormick  (1995)  states  that  the  minimum  induced  drag  coefficient  for  a  wing  is 


C 


(14) 


where  a  is  the  aspect  ratio.  Therefore  the  minimum  thrust  required  to  generate  lift  is  given  by 

rit 

^  (15) 


pA\^C, 
2na 


The  power  to  maintain  level  flight  with  this  ideal  wing  is  then,  in  accordance  with  equation  (15) 


pavXl  .  p 


(16) 
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The  rate  of  climb  for  a  typical  light  aircraft  is  well  approximated  with 

P-P 


'■  req 


mg 


The  use  of  (13),  (16)  and  (17)  yields  the  following  expression  for  the  exergy  of  lift. 

\2 


+ 


2(mg) 


. pV,Ana 

The  marginal  cost  of  mass  is  then 

dm  ^ 


(17) 


(18) 


(19) 


Note  that  ci^  is  not  the  coefficient  of  lift,  but  the  unit  exergy  cost  of  lift,  determined  via  exergy 
costing  (thermoeconomics).  Performing  the  differentiation  yields 


= 

m 


4g^/w 

pSAna 


'lift 
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(ii)  The  Creation  of  an  Exergy  Flow  Diagram 

The  expressions  for  the  exergy  of  lift  allow  exergy  flow  diagrams  to  be  developed.  In 

creating  these  Vagrams,  the  following  points  should  be  remembered: 

•  The  general  source  of  exergy  is  the  fuel.  In  general,  this  will  be  converted  in  the  propulsion 
system  to  thrust,  shaft  power,  compressed  air,  etc. 

•  Thrust  is  used  to  produce  lift  (overcome  induced  drag),  overcome  parasitic  drag,  overcome 
trim  drag*^  and  provide  ram  air.  In  the  exergy  flow  diagram,  thrust  exergy  will  flow  out  of 
the  aircraft  and  back  into  various  portions  of  the  aircraft. 

•  Lift  flows  from  the  wing  to  the  remainder  of  the  a  j  craft.  The  amount  of  lift  an  item  receives 
is  assumed  proportional  to  its  weight.  As  this  is  definitely  the  case  for  climb,  it  seems 
reasonable  to  assiune  this  for  all  conditions.  If  the  aircraft  is  at  a  steady  altitude,  the  exergy 
of  lift  is  consumed  by  the  “aircraft  pieces”  in  the  process  of  holding  them  aloft. 

•  These  exergy  diagrams  will  vary  with  the  loading  of  the  aircraft.  Therefore,  the  exergy  flow 
diagram  will  continually  change  during  flight.  An  exergy  flow  diagram  for  the  average  flight 
condition,  although  a  good  start,  may  not  always  be  adequate. 

Simple  exergy  flow  diagrams  for  the  Glastar  aircraft  are  shown  in  the  Case  Studies. 


''  An  exergy  balance  on  a  mass  in  flight  takes  the  form  ~  ■  Substitution  of 

(13)  yields  Xs  =  -mgVyQ. 

In  order  to  control  the  airplane  and  overcome  the  moment  produced  by  the  wing,  fixed-wing  aircraft 
generally  have  a  stabilizer.  Drag  produced  in  producing  the  lift  of  this  stabilizer  (which  actually  normally 
is  directed  downward)  is  called  trim  drag. 
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(b)  Choosing  the  ‘‘Reference  Environment’’  for  Design  and  Costing 

For  exergy  calculations,  it  is  necessary  to  choose  a  reference  environment.  For  a 
stationary  plant,  if  the  conditions  of  its  surroundings  change  significantly  during  its  operation,  it 
may  be  necessary  to  integrate  over  a  variety  of  environmental  conditions,  or  to  weight,  by 
percentage  of  time  in  a  given  environment,  finite  blocks  of  conditions  using  average  conditions 
for  that  block  (say,  with  the  temperature  fi'om  0-5°C,  5-10°C  and  so  on).  An  aircraft,  as  stated 
previously,  operates  in  a  wide  variety  of  conditions.  Not  only  will  different  reference  conditions 
be  used  for  different  modes  of  flight  (cruise,  climb,  etc.),  but  also  it  may  be  necessaiy'  to  integrate 
over  changing  conditions  in  a  single  mode  of  flight.  An  example  of  this  is  a  climu  to  cruise 
altitude. 

However,  it  will  be  argued  here  that  the  outdoor  environmental  conditions  are  not  always 
correct  as  a  reference  “dead  state”  for  evaluating  a  system  or  accounting.  This  argument  will  be 
made  through  the  example  of  a  candidate  system  for  Scenario  2  in  3.03,  the  “Beer  and  Soup” 
plant.  Figure  3  shows  this  system  operating  at  its  design  condition’^,  along  with  associated 
exergy  flows  and  cost  flows,  with  a  reference  temperature  of  288K.  The  management  of  the 
“Beer  and  Soup”  plant  evaluated  the  different  systems  available  to  them,  and  decided  that  this 
particular  system  was  in  the  plant’s  best  interest. 

The  “Beer  and  Soup”  company’s  cost  accountants  require  that  both  the  beer  department 
and  the  soup  department  pay  for  their  exergy.  However,  it  is  observed  that  as  the  outdoor 
temperature  changes  (which  is  being  used  as  a  reference  temperature)  exergy  charges  \vill  vary  to 
the  departments,  as  shown  in  Figure  4*'‘.  The  department  heads  argue  that  the  charges  should  be 
the  same  regardless  of  the  weather.  They  say  that  their  beer  and  soup  are  always  essentially  the 
same  temperature,  and  the  system  operating  cost  is  not  varying.  Furthermore,  the  management  of 
the  “Beer  and  Soup”  company  evaluated  the  different  systems  available  to  them,  and  decided  that 
this  particular  system  was  in  the  plant’s  best  interest.  The  system  had  been  chosen  as  the  best 
overall  for  the  company  as  a  whole. 


The  thermodynamic  and  thermoeconomic  models  are  given  in  Appendix  B.  Note  that  the  turbine 
efficiency  is  zero.  The  author  argues  that  the  proper  method  for  costing  when  a  valve  is  used  in  lieu  of  a 
turbine  is  as  a  turbine  with  zero  efficiency.  This  method  allows  capital  to  be  played  off  against  efficiency 
during  the  generic  optimization.  -  ' 

The  thermodynamic  model  and  the  costing  equations  are  given  in  Appendix  B. 
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Figure  3:  "Beer  and  Soup"  Design-Point  Exergy  Flows 

The  available  energy  theory  of  Gibbs  (Gaggioli  et.  al.,  1999),  along  with  an  examination 
of  the  system’s  constraints,  resolves  the  apparent  paradox  of  varying  charges  to  the  departments 
without  a  variance  in  the  overall  operating  cost  of  the  system. 
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Figure  4:  Variation  in  Exergy  and  Cost  Flows  with  Dead  State  Temperature 


(i)  System  Information 

Before  moving  on,  it  will  be  necessary  to  have  more  information  about  the  beer  and  soup 
plant.  At  base  load,  it  is  designed  to  cool  3 1.9  kg/s  of  beer  from  288K  to  276K  and  heat  33.3  k^s 
of  soup  from  300K  to  34 7K.  (Calculations  of  the  log  mean  temperature  of  each  will  yield  the 
temperatures  of  the  beer  and  soup  from  previous  exaniples  of  this  plant.)  The  plant’s  circulating 
fluid,  R1 14,  has  a  mass  flow  of  88.5  kg/s  and  operate.:  between  862  and  40.3  kPa.  The 
compressor  requires  5000  kW  of  shaft  power  and  has  an  isentropic  efficiency  of  70%. 

(ii)  Employing  Gibb 's  Available  Energy  Theory 

At  any  instant,  available  energy  for  an  object,  or  collection  of  objects,  is  defined  by  Gibbs 
as 


A  =  E  —  E^^^  (Constraints^  (21) 

Given  the  instantaneous  values  of  the  constraints,  the  available  energy  is  obtained  by  reducing  the 
energy  to  the  minimum  possible  values  consistent  with  the  constraints.  The  “dead  state”  is  the 
minimum-energy  state,  the  constrained  equilibrium  state,  at  these  values  of  the  constraints. 
Therefore,  in  order  to  find  the  dead  state,  the  constraints  must  be  defined. 

11  Constraints 

The  beer  and  soup  plant  never  interacts  thermally  with  the  ambient  environment.  It 
interacts  with  only  the  beer  and  soup  department.  Its  performance  (i.e.,  power  requirement  and 
exergy  delivery)  depends  only  on  the  beer  and  soup  temperatures.  The  composite  of  the  beer  and 
soup  entering  the  system  has  available  energy  because  the  beer  and  soup  are  not  in  equilibrium 
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with  each  other.  The  purpose  of  the  plant  is  to  increase  the  temperature  of  the  soup  and  decrease 
the  temperature  of  the  beer,  driving  them  further  out  of  equilibrium  and  hence  increasing  the 
available  energy  of  the  composite. 

A  typical  rej&igeration  system’s  working  fluid  will  be  in  a  two-phase  state  at  normal 
operating  temperatures  when  the  system  is  not  running.  Therefore,  the  dead  state  of  the 
refrigerant  should  be  saturated  at  the  reference  temperature.  (Wepfer  and  Gaggioli,  1980) 

Finally,  chemical  exergy  should  not  be  considered  for  this  system.  As  long  as  the 
refrigerant  is  contained  inside  of  the  system,  it  is  constrained  from  chemical  reactions. 


21  Calculation  of  the  Equilibrium  State 

If  we  assume  the  two  fluids  are  incompressible  and  have  constant  specific  heats  and 
assuming  a  reference  temperature  of  absolute  zero,  the  energy  flow  may  be  defined  as 

E  =  C^er^beer  +  C^oup'^^oup  (22) 

where  C  =  mc^ .  Then  is  foimd  by  minimizing  (22)  with  the  constraint  of  the  entropy  flow 
remaining  constant.  The  constraint  yields 


In^^  +  G  ln-^=0 


soup 


^ beer^eq 


(23) 


Solving  equation  (23)  for  Tsoup.eq  and  substituting  into  (22), 


E  =  C^  T.  +C  T 

beer  beer,eq  soup  soup 


T 

V  beer, eg  J 


(24) 


The  minimum  energy  is  obtained  by  taking  the  derivative  of  this  with  respect  to  Tbeer.eqmi  setting 
it  equal  to  zero. 


dE 

0  = - -T 

Solving  for  yields 


T 

\  beer. eg  J 


^  beer, eg 


beer.eq  soup 


J 


(25) 


(26) 


As  expected,  when  this  result  is  substituted  into  (23)  to  get  Tsoup.eq^^  same  result  is  obtained. 
For  the  conditions  of  the  beer  and  soup  entering  the  plant,  this  temperature  is  294K.  For  the  beer 
and  soup  exiting  the  plant  it  is  3  lOK.  In  practice,  either  of  these  would  suffice  as  an  appropriate 
dead  state  temperature;  in  theory  a  proper  average  is  required. 


31  Calculation  of  the  Theoretical  Dead  State  Temperature 
The  reason  that  the  equilibriiun  temperature  changes,  as  the  beer  is  cooled  and  the  soup  is 
heated,  is  because  the  increase  in  available  energy  of  the  beer-soup  composite  is  less  than  the  net 

energy  input  (the  compressor  power).  That  is,  because  is  increasing  as  a  result  of  entropy 
production.  The  change  in  is  given  by  dE^  =  T^^dS  where 

Sou,  -  =  (Ssoup.ou,  -  S,^p.,n)  +  {Sbeurpu,  ”  ^quals  the  ovetall  rate  of  entropy  production. 

Thus, 
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AE  =E  -E  =  dS 

^^min  ^mirijout  ^minjn  ^ 


(27) 


Hence  the  appropriate  average  equilibrium  temperature  is  equal  to 

7;=^""”/.  (28) 

This  temperature  is  that  of  a  hypothetical  entropy  reservoir  from  which  the  entropy 

could  be  obtained  (i)  while  transferring  a  net  amount  of  energy  to  the  system  (iij 

with  no  net  transfer  of  available  energy  used  in  order  to  change  .  This  is  illustrated  in  Figure 


5;  the  area  imder  the  curve  is  .  When  To  equals 


AE. 


•  ,  the  two  shaded  areas  are  equal. 


Available  energy  obtained  while  raising  Teg  from  Tgg.m  to  To  is  returned  upon  raising  Teg  from  To 

to  Teg,out‘ 


Finally,  inasmuch  as 


and 


p  _  p  zz  C  ( T 

^mm,out  ^ min, in  ^ total  eg, out 


^min  put  ^  min  Jn 


c. 


total 


(29) 

(30) 


This  is  the  logarithmic  mean  of  the  equilibrium  temperatures,  as  calculated  by  (26),  before  and 
after  the  process  occurs.  For  this  example,  that  temperature  is  302K.  The  corresponding  exergy 
flow  and  costing  is  shown  in  Figure  6. 


Figure  5:  Theoretical  T-S  Curve 

As  a  closure  to  the  “beer  and  soup”  story  the  following  should  be  mentioned.  After  the 
previous  calculations  had  been  performed,  the  beer  department  manager,  well-trained  in  exergy, 
calculated  that  his  heat  exchanger  efficiency  was  only  5.5%  while  his  colleague’s,  in  the  soup 
department,  had  an  efficiency  of  26%.  He  asked,  “Why  should  I  be  penalized  with  higher  cost 
per  unit  exergy  just  because  management  stuck  me  with  a  low  efficiency  heat  exchanger?’’ 
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Indeed,  he  was  correct.  Costing  was  revised  to  the  equivalence  method,  with  the  unit  price  of 
exergy  the  same  for  both  the  beer  and  soup  departments. 


Figure  6:  Exergy  Flows  and  Unit  Costs,  Calculated  with  Equilibrium  Temperature 


(c)  Mai^inal  Exergy  Costs 

The  most  common  exergy  cost  is  an  average  cost,  which  includes  capital  costs.  This  cost 
is  at  times  misused,  such  as  when  operating  decisions  are  being  made  for  an  existing  plan  (El- 
Sayed  and  GaggioU,  1989).  At  such  time  the  proper  cost  is  the  marginal  exergy  cost, 


(32) 


There  is  another  instance,  essentially  unnoticed  before,  when  marginal  costs  should  be 
used.  The  marginal  cost  is  also  the  proper  costing  for  detailed  design,  as  during  detailed  design 
the  device  or  subsystem  design  team  is  operating  under  the  assumption  that  the  remainder 
equipment  for  the  system  is  fixed  The  team 's  design  has  no  effect  on  the  capital  costs  of  the 
remainder  of  the  plant.  These  capital  costs  must  therefore  be  considered  sunk. 

Consider  once  design  scenario  one,  the  aircraft  example.  Including  capital  costs  in  the 
cost  equations  (which  will  be  described  in  detail  in  the  case  studies)  yields  a  cost  of  weight  of 
himdreds  of  dollars.  This  high  of  a  cost  is  sensible  diuing  preliminary  design,  as  reduced  weight 
can  reduce  engine  power  requirements  (as  well  as  stmcturd  requirements).  However,  during 
detailed  design,  an  engine,  or  at  least  range  of  engines  has  been  selected.  Weight  reduction  affect 
only  performance.  In  the  case  studies  it  will  be  seen  that  not  including  capital  cost  reduces  the 
cost  of  weigh  by  an  order  of  magnitude.  One  would  be  hard  pressed  to  find  an  owner  of  a  light 
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plane  willing  to  spend  hundreds  of  dollars  to  reduce  the  empty  weight  of  his  plane  by  a  couple  of 
pounds,  but  many  would  be  willing  to  spend  tens  of  dollars. 

The  most  accurate  analysis  would  integrate  marginal  cost  from  the  initial  design 
parameters  (passed  to  a  team  at  the  beginning  of  detailed  design)  to  the  actual  operating  point,  as 
the  marginal  cost  is  a  function  of  the  design  parameters.  If  deviations  are  small,  the  marginal 
costs  may  be  assumed  constant. 

Furthermore,  if  the  deviations  from  the  initial  design  parameters  remain  small, 
incremental  costs  may  ofien  be  used  in  place  of  the  marginal  costs.  Here,  the  incremental  cost  is 
defined  as  the  average  exe^'gy  cost  ignoring  capital  costs.  This  is  simply  the  marginal  cost  if 
equipment  performance  is  considered  constant. 

These  points  are  illustrated  with  the  following  example.  Consider  an  existing  combustion 
gas  turbine  power  plant.  The  owner  requires  an  additional  power  output  from  the  plant;  it  is 
decided  to  add  a  Rankine  bottoming  cycle  as  illustrated  in  Figure  7.  The  genericaUy  “optimized” 
initial  design  calls  for  “duct  firing”,  which  is  combustion  in  the  turbine  e^diaust  before  it  enters 
the  heat  recovery  steam  generator  (HRSG). 


Figure  7:  Rankine  Bottoming  Cycle 

A  component  design  team  is  responsible  for  selecting  a  suitable  steam  turbine-generator. 
They  are  given,  from  the  system  design  team,  the  design  values  from  the  generically  optimized 
system,  such  as  mass  flow,  condenser  vacuiun,  superheat,  expected  capital  and  exergy  costs,  etc. 
It  is  highly  unlikely  that  a  steam  turbine  exists  that  will  perfectly  match  these  values,  nor  does  it 
take  a  detailed  economic  analysis  to  realize  that  designing  and  custom  building  a  steam  turbine 
makes  no  economic  sense.  The  required  output  of  the  turbine-generator  is  fixed.  The  initial  total 
cost,  from  the  “generic  optimization”,  based  on  the  thermodynamic  model  and  generic  cost 
models,  is 
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'^initial  ^sl,milial  ^^l,miUaI 


(33) 


The  problem  of  selecting  a  steam  turbine-generator  becomes  finding  the  match  that 
satisfies  system-level  constraints  while,  hopefiilly,  decreasing  the  total  cost  over  the  initial  design 
value,  or,  if  that  proves  unfeasible,  minimizing  the  increase  in  cost.  With  the  approximation  of 
constant  margind  costs,  this  is  written  as  the  following  fimction  to  be  minimized. 

~  i^^sl  ~  ^sl,MtiaI )  ■*■  ^6  ("^6  ~  (^7  ~  ^IJmtiaJ )  (3“^) 

The  initial  design  terms  may  be  eliminated  from  the  problem,  as  they  have  been  determined  in  the 
generic  optimization  process. 

To  calculate  the  marginal  cost  of  steam  supplied,  the  cost  of  the  steam  is  first  calculated 
with  a  money  balance  around  the  HRSG  and  duct  firing. 

=  c^Xi  +  C2X2  +  +  Zdf  +  ^hesg  (35) 


It  is  assumed  here  that  the  feedwater  exergy  is  negligible  compared  to  that  of  the  combustion 
gases. 


The  marginal  cost  of  the  steam  supplied  is 


The  three  terms  cancel  because  the  flow  of  gas  turbine  exhaust  is  fixed  and,  at  this  stage,  the 
selection  of  the  steam  turbine  is  being  done  independently  from  that  of  the  HRSG  and  duct  firing 
design.*^  (Changes  in  the  amoimt  of  duct  firing  and  flow  of  cooling  water  through  the  condenser 
will  be  used  to  compensate  for  deviation  from  the  initial  design.) 

For  this  problem,  the  second-law  efficiencies  of  the  duct  firing  and  HRSG  are  defined  as 


‘Hhusg  ~ 


^3 


(37) 


Vdf  ~ 


X. 


(38) 


From  these  two  efficiencies,  and  once  again  with  the  assumption  that  the  feedwater 
exergy  is  negligible. 


^6  ~  ^HRSG^DF  (^1  "^2) 


(39) 


Rearranging  and  taking  the  partial  derivative,  assuming  constant  efficiencies,  yields 
^^6  ^HFSgVdF 


(40) 


As  discussed  in  Section  3.03,  before  detailed  design  is  complete,  simulation  and  proofing  are  done  to 
ensure  that  independently  designed  or  selected  devices  function  properly  together.  Moreover,  an  additional 
optimization  on  operating  parameters  may  be  performed. 
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and 


^hrsgVdf 

It  will  now  be  demonstrated  that  this  is  a  case  that  cannot  be  well  approximated  by  the 
incremental  cost.  As  gas  turbine  exhaust  is  currently  being  “thrown  away”,  the  cost  associated 
with  it  is  zero.  The  cost  balance  (35)  without  capital  cost  is  then 

CgXg  =  yfXy  +  C2X2 + y^X^  (42) 

Substituting  (41)  from  above,  and  solving  for  Cg  yields 


^6  = 


VhrsgHdf  ^2 


(43) 


The  reason  for  this  inaccmacy  is  that  in  this  case,  all  additional  exergy  as  steam  can  only 
be  had  at  the  expense  of  additional  duct  firing.  Equation  (43)  would  assume  Aat  it  would  be  split 
from  turbine  exhaust  and  duct  firing.  However,  there  is  no  physical  way  to  increase  the  gas 
turbine  exhaust‘d. 

Similar  methodology  would  be  used  to  find  c'-j .  After  selection  of  the  individual  devices 

in  this  system  is  complete,  the  design  team  moves  on  to  the  simulation  and  proofing  phase.  At 
this  stage,  it  would  be  desirable  to  re-optimize  on  operating  variables.  For  this  system,  these 
could  include  superheat,  operating  pressure  and  condenser  vacuum. 

Marginal  exergy  costs  are  the  proper  exergy  costs  to  use  in  detailed  design.  Additionally, 
if  the  last  stage  of  a  global-local  optimization  scheme  were  a  local  optimization,  the  use  of 
marginal  costs  here  could  likely  improve  accuracy.  In  many  cases,  but  not  all,  incremental  costs 
may  be  used  in  place  of  marginal  costs.  However,  if  the  engineer  does  not  have  a  good  physical 
sense  of  how  the  system  performs,  it  is  wise  to  perform  some  mathematical  investigation  as  was 
performed  above. 


VIII.  CASE  STUDIES  FOR  A  LIGHT  AIRCRAFT 

The  case  studies  that  follow  are  all  based  on  scenario  1,  as  discussed  in  Section  3.03.  The 
aircraft  is  the  Glastar  aircraft  (see  Appendix  A  for  details).  Here,  an  overall  objective  fimction 
will  be  formulated.  Exergy  flow  and  cost  diagrams,  as  required  to  begin  the  detailed  design 
process,  are  developed.’’  Then,  the  applicability  of  the  methods  previously  developed  here  to  the 
detailed  design  process  wiU  be  demonstrated  through  alternator  and  engine  selection. 


(a)  A  User  Questionnaire 

In  order  to  construct  objective  functions  it  was  necessary  to  question  the  user  about  the 
value  of  different  desiderata. 


An  analagous  scenario  would  occur  when  an  aircraft  engine  must  use  its  afterburner  in  order  to  increase 
thrust. 

”  When  detailed  design  is  discussed,  it  is  assumed  that  “generic”  optimization  has  selected  the  160 
horsepower  Lycoming  0-320  engine  and  a  constant-speed  propeller. 
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It  was  decided,  based  on  the  light  aircraft’s  typical  “mission”  profile,  to  consider  full 
power  cruise  at  8000  ft.,  economy  cruise  at  8000  ft.  and  climb  rate  as  desiderata.  Climb  rate  was 
considered  a  measure  of  short  field  performance.  Although  it  would  have  been  advantageous  to 
add  a  more  direct  measure  of  this  part  of  the  mission,  it  would  have  involved  using  a  far  more 
complicated  aircraft  simxilation.  A  summary  of  conditions  for  the  measiuement  of  the  various 
desiderata  and  their  preliminary  design  values  is  given  in  Table  6**. 


Cargo 
Mass  (kg) 

r  I  ? 

Speed 

(km/hr) 

Climb 

(m/min) 

Fuel  Bum 
(kg/hr) 

Full 

power 

cruise 

204 

54 

259 

0 

22.9 

Economy 

cruise 

221 

54 

244 

0 

18.8 

Climb 

221 

109 

147 

372 

38.1 

Table  6:  Preliminary  Design  Values 

The  following  questions  were  posed  to  the  end-user,  and  the  following  answers  received. 

1.  When  you  were  selecting  an  aircraft  kit,  what  was  the  range  of  advertised  cruise  speeds  of  the 
aircraft  (at  75%  power,  8000  ft.)  in  which  you  were  interested?  Answer:  120  to  200  mph 

2.  Approximately  how  much  additional  money  would  you  be  willing  to  spend  for  the  aircraft  in 
order  to  increase  the  cmise  speed  of  your  airplane  by  5  mph  above  the  average  of  the  figures 
in  question  1?  Answer:  $350 

3.  When  you  were  selecting  an  aircraft  kit,  what  was  the  range  of  advertised  ranges  of  the 
aircraft  in  which  you  were  interested  (at  65%  power,  8000  ft.)?  Answer:  500  to  1000  miles 

4.  Approximately  how  much  money  would  you  be  willing  to  spend  in  order  to  increase  the 
range  of  your  airplane  by  50  miles  over  the  average  of  the  figures  in  question  3?  Answer: 
$300 

5.  When  you  were  selecting  an  aircraft  kit,  what  was  the  range  of  advertised  climb  rates  of  the 
aircraft  in  which  you  were  interested?  AnsM’er:  700  to  1800 fpm 

6.  Approximately  how  much  money  would  you  be  willing  to  spend  in  order  to  increase  the 
climb  rate  of  your  airplane  by  100  fpm  above  the  average  of  the  figures  in  question  1? 
Answer:  $250 

For  all  calculations,  it  was  assumed  that  the  aircraft  would  have  a  lifespan  of  20  years,  be 

flown  200  hours  per  year,  with  50%  of  those  hours  at  economy  cruise  and  the  remaining  50%  at 

high-speed  cruise.  Hours  spent  in  climb  are  considered  negligible.  An  interest  rate  of  8%  was 

used. 


**  These  values  are  from  a  simple  aircraft  performance  model.  The  program  itself  is  in  Appendix  B.  The 
values  computed  for  speed  and  range  are  very  close  to  manufacturer’s  test  data.  However,  climb  rates 
calculated  by  the  program  appear  approximately  300  fpm  too  low.  For  the  purpose  of  the  example,  values 
from  the  model  will  be  used  as  is. 

Although  exergy  and  performance  calculations  for  economy  cruise  were  performed  with  half  fuel  so  as  to 
approximate  an  average,  range  calculations  were  performed  with  the  flight  beginning  with  fijll  fuel  (109 
kg) 
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(b)  Cost-benefit  Curves 

With  this  information  linear  functions  for  the  three  desiderata  were  created,  and  are 
shown  in  Figures  8-10. 


Figure  8:  Value  versus  Cruise  Speed 


Figure  9:  Value  versus  Rate  of  Climb 


Figure  10:  Value  versus  Range 
(c)  The  Overall  Objective  Function 

The  combination  of  the  three  cost-benefit  curves  allows  an  overall  objective  fimction  to 
be  formed.  The  value  of  each  of  the  desiderata  may  be  written  as: 

=  {Kruise  -  ^  120mph  <  <  200mph  (44) 

=  5600dol  €  200mph  < 

^Rc  =  ^  700fpm>i?C 

Z^c  =  i^C  -  700mph)2.5-^  e  700fpm  <RC<1 800fpm  (45) 

fpm 

Zjfc  =  2750dol  G 1  SOOfpm  <RC 
=  Odol  e  500mile>i? 

Z„  =  (7?  -  500mph)6-^^  e  SOOmile  <R<  lOOOmile  (46) 

mile 

Zj^  =  SOOOdol  €  lOOOmile  <  R 

Here  Fc  represents  cruise  speed  at  75%  power  and  8000  ft.,  RC  the  rate  of  climb  at  seal  level  and 
R  the  range  at  65%  power  and  8000  ft. 

Now  the  objective  fimction,  to  be  maximized,  may  be  written  as 

~  ^RC  “  ^Capital  ■*"  ^mem  (47) 

It  is  desired  to  maximize  the  value  of  the  aircraft’s  performance  minus  the  cost  of  the  aircraft. 
Zmean  is  the  mean  value  of  a  light  aircraft  to  the  user. 

Although  the  questioning  of  the  user  did  not  determine  Zmem,  (47)  is  adequate  for  making 
decisions  between  candidate  systems  and  for  preliminary  design,  as  Zmem  is  a  constant  which  will 
cancel  in  compressions  or  optimization. 
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(d)  Exergy  Flow  Diagrams 

Figures  1 1, 12  and  13  show  the  exergy  flows  and  per  unit  exergy  costs  in  a  light  aircraft 
for  full  power  cruise,  economy  cruise  and  climb.  Exergy  destructions  are  also  given. 

These  figmes  were  developed  by  (i)  writing  an  exergy  balance  on  each  device,  (ii) 
writing  cost  balances  on  each  device  and  adding  auxiliary  equations.^®  Additionally,  it  is 
necessary  to  fix  a  unit  cost  of  exergy  for  either  thrust  or  fuel;  this  was  done  either  equation  (9)  or 
with  fuel  cost  and  equation  (10),  depending  on  the  flight  conditions. 

Additionally,  engine  and  propeller  performance  relations  were  used.  These  were  in  the 
form  of  a  user’s  manual  for  the  Lycoming  engine  and  a  computerized  propeller  map  supplied  by 
Hartzell.  Aircraft  performance  was  estimated  from  airfoil  and  manufacturer’s  test  data. 


Figure  11;  Exergy  Flows  and  Unit  Costs,  75%  Power,  8000  ft. 


In  the  interest  of  conserving  space,  the  systems  of  equations  are  listed  in  Appendix  B.  A  good 
explanation  of  the  methods  employed  to  develop  the  equations,  along  with  simple  examples,  is  found  on 
pages  3-50  of  Efficiency  and  Costing,  1983,  edited  by  R.  A.  Gaggioli.  An  additional  valuable  reference  is 
Bejan  et  al.,  1996. 
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Figure  12:  Exergy  Flows  and  Unit  Costs,  65%  Power,  8000  ft. 


Figure  13:  Exergy  Flows  and  Unit  Costs,  Maximum  Climb  at  Sea  Level 

(e)  Application  to  Subsystem  Design/Selection:  Alternator  Selection 

Here  it  will  be  decided  whether  it  is  better  to  purchase  a  “standard”  alternator  or  a 
lightweight  model.  As  no  efficiency  data  is  available,  both  will  be  assumed  equally  efficient. 

The  standard  alternator  has  an  initial  cost  of  $294  and  a  mass  of  5.9-kg  (13-lbm).  The  alternative 
costs  $450,  but  has  a  mass  of  only  2.7-kg  (6-lbm).  The  selection  will  first  be  made  using  the 
foregoing  decomposition  process.  The  results  will  be  checked  by  application  of  the  overall 
objective  function.  (The  relative  simplicity  of  modeling  a  light  aircraft  makes  this  application 
straightforward.) 

To  make  this  decision  via  decomposition,  equation  (1 1)  is  used.  Terms  for  shaft  power 
are  unnecessary,  as  the  power  inputs  to  both  units  are  the  same.  The  cost  of  switching  fi’om  the 
standard  to  the  lightweight  alternator  will  be  given  by 
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^lightweight  ^standard 

^i^mi  {^^lighlweighl  ^standard  ) 

/ 

The  lifespan  of  20  years,  with  an  interest  rate  of  8%  and  200  operating  hours  per  year 
yields  0.229  and  0.150  dollars/hr  for  values  of  Z  for  the  lightweight  and  standard  alternators, 
respectively. 

Using  the  performance  values  from  Table  1,  along  with  the  information  from  Figures  8- 
10^',  yields  weighting  factors  of  0.565  for  high-speed  cruise,  0.283  for  range  and  0.153  for  climb 
rate. 

Equation  (20)  is  used  to  find  the  costs  of  mass,  with  the  cost  of  lift  taken  from  the 
preceding  exergy  flow  diagrams.  This  yields  marginal  mass  costs  of  4.17-10'^,  3.33*10’^  and 
5.32- 10"^  dollars/kg  hr  for  full  power  cruise,  economy  cruise  and  rate  of  climb,  respectively. 

Finally,  evaluating  equation  (48)  with  these  values  yields  an  increase  in  cost  for  the 
aircraft  of  0.0608  dollars/hr,  or  a  present  value  of  1 18  dollars,  associated  with  the  lightweight 
alternator.  Therefore  it  is  not  desirable  to  use  the  lightweight  alternator  over  the  standard  model. 

The  foregoing  conclusion  was  checked  by  direct  application  of  the  overall  aircraft 
objective  function.  Evaluating  the  aircraft  performance  by  simulating  the  aircraft  with  the  model 
used  for  design  point  calculations,  along  with  the  overall  objective  function  in  equation  (47) 
yields  an  increase  of  cost  of  12 1  dollars  for  the  lightweight  alternator.  In  terms  of  overall  aircraft 
cost,  or  even  alternator  cost,  the  difference  between  the  values  found  through  decomposition  and 
through  overall  simulation  is  insignificant. 

(f)  Application  to  a  Propulsion  System:  Engine  Selection 

The  160-hp  engine  is  the  “standard”  engine,  chosen  most  widely  by  builders.  The  180-hp  engine 
will  produce  gains  in  full  power  cruise  and  climb  rate,  but  at  the  expense  of  increased  weight, 
decreased  range,  increased  fuel  consumption  and  increased  capital  cost.  (Table  4) 


Engine  (hp) 

Cost 

(k  dollars) 

Mass  (kg) 

Full  Power 
Cruise 
(km/hr)/ 
Fuel  Bum 
(kg/hr) 

Fuel  Bum  @ 
244  km/hr 
(kg/hr) 

Rate  of 
Climb 
(m/min) 

160 

22.3 

116 

259/22.9 

18.8 

372 

180 

26.2 

122 

271/29.9 

21.8 

442 

Table  7:  Performance  with  different  engines 


This  information  may  be  applied  directly  to  the  overall  objective  function,  which  results  in 
decrease  of  profit  (present  value)  of  10,800  dollars. 

Alternatively,  equation  (12)  may  be  employed  in  the  form 


The  slopes  of  these  figures  yield  the  values  of  $350/5-mph,  $3 00/50-miles  and  $250/100-fpm. 
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^ power  J  ('^pcitw, 1,180  ^  pan'erj,]6Q^ 

tj 


(49) 


The  unit  costs  of  power  (i.e.,  exergy)  may  be  taken  from  the  exergy  flow  diagrams  and  the 
weighting  factors  remain  the  same  as  the  previous  examples.  The  time  factors  are  0.5  for  full 
power  cruise,  0.5  for  economy  cmise  and  0  for  climb.  Fuel  costs  are  calculated  on  a  per  unit 
mass  basis  here,  as  fuel  is  the  raw  energy  source.  The  value  used  is  0.735  dollars/kg.  (Exergy 
costing  with  the  fuel  exergy  costs  assigned  would  yield  tiie  same  values.)  The  cost  of  fuel  bum  to 
range  must  be  added  to  the  cost  of  fuel  for  evaluation  of  range  at  economy  cruise,  as  per  equation 
(10),  increasing  the  per  kilogram  cost  of  fuel  to  0.912  dollars. 

Evaluation  of  equation  (49),  shows  a  decrease  of  present  value  profit  of  1 1,400  dollars  for 
the  180-hp  engine.  Although  this  differs  from  the  10,800-dollar  non-decomposed  value  by  5.5%, 
the  same  decision  would  be  reached.  The  error  likely  results  from  the  large  shift  from  the 
preliminary  design  values,  shifts  in  engine  power  of  12%,  (engine)  weight  of  6%  and  fuel 
consumption  of  14%.  In  such  cases  iteration  in  the  exergy  costing  step  would  improve  results  if 
necessary. 


IX.  CONCLUSIONS 

A  general  outline  of  the  design  process  has  been  laid  out  in  a  manner  that  is  applicable  to 
the  design  of  thermal  systems  in  general,  and  vehicular-based  energy  systems  in  particular.  In 
support  of  this  general  process,  means  of  constmcting  an  overall  objective  function  have  been 
discussed.  One  area  in  the  thermal  design  process  that  has  been  particularly  lacking  is  the 
detailed  design  process.  Much  work  has  been  accomplished  in  the  area  of  generic  optimization; 
this  leaves  off  before  an  actual  working  system  is  penned.  A  thermoeconomic  approach  to 
detailed  design  is  presented  here.  Finadly,  notes  are  made  (i)  on  the  selection  of  a  “dead  state” 
for  exergy  analysis  (ii)  as  to  when  marginal  exergy  costs  should  be  used. 

(a)  The  Overall  Objective  Function 

One  method  was  outlined  for  finding  an  overall  objective  function.  In  any  case,  the 
subjectivity  of  the  customers’  desires  for  a  product  must  be  objectified  for  optimization  to  occur. 
Without  an  overall  objective  function  for,  for  example,  an  aircraft,  it  is  impossible  to  develop 
objective  functions  to  guide  engineers  in  the  design  of  energy  systems,  subsystems  and 
conversion  devices. 


(b)  Detafled  Design 

Second-law  based  decomposition  strategies  do  not  only  aid  in  mathematical,  “generic” 
optimization,  but  allow  detailed  design  to  proceed  in  a  logical,  efficient  manner.  It  allows 
objective  functions  to  be  written  for  subsystems  and  devices.  Without  these  objective  functions, 
rational  decisions  may  not  be  made  during  the  design  of  these  subsystems  and  devices  without 
recourse  to  the  simulation  of  the  whole  system  or  vehicle.  Such  simulation  is  time  and  resource 
consuming.  Within  this  dissertation,  objective  functions  for  aircraft  energy  systems  were 
developed  in  detail. 

So  as  to  apply  exergy  concepts  to  aircraft,  exergy  flow  diagrams  were  developed  for  a 
single-engine  general  aviation  aircratft.  To  do  this  it  was  necessary  to  derive  an  expression  for 
the  exergy  of  lift.  These  exergy  flow  diagrams,  along  with  thermoeconomics  were  used  to  place 
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costs  on  weight,  and  to  add  costs  to  fuel  to  account  for  range  as  a  desideratum.  It  was  also  shown 
how  to  develop  weighting  factors  for  different  realms  of  performance. 

Although  the  objective  functions  for  detailed  design  as  presented  here  were  developed 
specifically  for  aircraft  design,  they  are  readily  extensible  to  other  vehicles  as  well  as  stationary 
systems. 

(c)  The  “Dead  State” 

Before  choosing  a  “dead  state”  for  exergy  calculations,  the  constraints  on  the  system 
must  be  considered.  For  oxie  such  set  of  constraints  an  example  of  the  calculation  of  the  dead 
state,  using  the  available  energy  theory  of  Gibbs,  was  made. 

The  calculation  was  made  here  for  a  system  with  fixed  constraints  independent  of  the 
environment;  in  general,  these  constraints  will  vary  with  time  and  may  be  dependent  on  the 
environment  This  is  often  the  result  of  changing  environmental  conditions,  such  as  varying 
weather,  or  the  system  itself  moving  (such  as  an  aircraft  in  climb).  If  this  is  the  case,  integration 
of  exergy  calculations  over  time  or  weighting  of  calculations  with  different  dead  states  may  be 
necessary. 

A  system  vwth  constraints  independent  of  the  environment  need  not  be  a  fixed  plant, 
although  the  example  in  this  dissertation  was.  The  equilibrium  temperature  as  defined  through 
available  energy  theory  has  applicability  to  vehicle  design  as  well.  An  aircraft  subsystem  may  be 
operating  between  two  thermal  “sinks”,  for  example,  the  avionics  and  the  fuel,  and  never  “see” 
the  external  environment. 

The  fundamental  conclusion  is  this.  The  available  energy  concept  introduced  by  Gibbs  is 
more  fundamental  than  the  exergy  concept,  which  is  a  special  case  of  available  energy  of  a  “body 
and  a  (large)  medium”  in  Gibbs  terminology.  As  shown  by  Gaggioli  (1995, 1998),  in  any  case  (i) 
the  equations  representing  contributions  of  subsystems  to  the  overall  available  energy  of  a  system 
can  be  derived  without  reference  to  any  environment;  (ii)  for  each  subsystem,  that  equation 
includes  properties  of  that  subsystem’s  dead  state;  (iii)  the  subsystem’s  dead  state  properties  are 
found  by  minimization  of  the  overall  system  energy  subject  to  the  constraints  on  each  of  the 
subsystems  while  maintaining  overall  system  entropy.  The  example  presented  here  represents  an 
application  of  that  theory,  and  shows  an  example  when  that  theory  is  needed,  practically. 

(d)  Marginal  Exergy  Costs 

When  the  design  process  moves  to  detailed  design  marginal  exergy  costs  must  be  used  by 
the  device  design  teams.  These  may  often  be  approximated  by  incremental  costs.  However,  one 
example  was  given  when  the  incremental  costs  would  not  be  proper,  a  combined-cycle  power 
plant  with  duct  firing.  Many  times  it  will  be  obvious  if  approximations  may  be  made.  If  it  is  not 
obvious,  mathematical  checks  should  be  made  before  the  assumptions! 


(e)  Case  Studies 

The  tools  for  detailed  design,  namely  the  creation  of  exergy  flow  diagrams  for  an  aircraft 
and  the  development  of  overall  and  device  objective  fimctions  have  been  demonstrated  in  a  real- 
world  case  study  of  a  homebuilt  aircraft. 

(f)  Future  Work  and  Recommendations 
Future  work  should  be  centered  in  two  areas: 
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1.  The  development  of  the  overall  objective  function  for  combat  aircraft.  In  this  dissertation, 
metiiods  by  which  this  may  be  approached  have  been  proposed,  but  not  tested.  These 
proposals  should  be  tested  and  refined. 

2.  The  device  and  subsystem  objective  functions  should  be  applied  to  past  decisions  in  industry'. 
Ideally,  the  decisions  should  have  involved  computing  the  effects  on  the  aircraft’s 
performance.  This  will  help  validate  the  methodology.  The  aircraft  would  ideally  be  a 
supersonic  aircraft,  so  as  to  ensure  that  the  exergy  of  lift,  as  developed  here,  is  applicable  to 
supersonic  flight. 
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APPENDIX  A:  GLASTAR  AIRCRAFT  INFORMATION^ 


Psrforma-nce 

125  h.p. 

160  h.p. 

180  h.p. 

Top  Speed  (TAS  at  sea 
level) 

156  m.p.h.  /  136 

167  m.p.h.  /  145 

171  m.p.h.  /  149 

kts. 

kts. 

kts. 

Cruise  Speed  (TAS) 

75%  power  at  8,000  ft. 

151  m.p.h.  /  131 
kts. 

161  m.p.h.  /  140 
kts. 

167  m.p.h.  /  145 
kts. 

65%  power  at  8,000  ft 

140  m.p.h.  /  122 

153  m.p.h.  /  133 

159  m.p.h.  /  138 

kts. 

kts. 

kts. 

Staii  Speed  (at  max. 
gross) 

No  flaps  (Vs) 

56  m.p.h.  /  49  kts. 

Full  flaps  (Vso) 

Rate  of  Climb 

49  m.p.h.  /  43  kts. 

Solo 

1,300  f.p.m. 

2,075  f.p.m. 

2,150  f.p.m. 

Max.  gross 

Range  (at  65%  power) 

1,000  f.p.m. 

1,390  f.p.m. 

1,500  f.p.m. 

Standard  tanks 

598  mi.  /  520  n.m. 

553  mi.  /  481  n.m. 

436  mi.  /  379  n.m. 

With  auxiliary  tanks 

1,021  mi.  /  888 
n.m. 

953  mi,  /  829  n.m. 

762  mi.  /■663  n.m. 

Fuei  Consumption  (at 
65%  power) 

5.8  g.p.h. 

6.7  g.p.h. 

8.5  g.p.h. 

Service  Ceiling 
(estimated) 

17,000  ft. 

20,000  ft. 

21,000  ft. 

Specifications 

Fuselage  Length 

Continental  10-240  (rigged  for  flight) 

Continental  10-240  (wings  folded) 

Lycoming  engine  (rigged  for  flight) 

Continental  10-240  (wings  folded) 

Wing  Span 

Rigged  for  flight 

Wings  folded  &tail  removed 

Other  Wing  Data 

Airfoil 

Area 

Aspect  ratio 

Wing  loading  (at  max.  gross) 
Structural  limit  loads  (at  max.  gross) 

Maximum  Height 


22.3  ft. 
24.5  ft. 
22.8  ft. 
25.0  ft. 


35.0  ft. 

8.0  ft. 

NASA  GA(W)-2 
128.0  sq.  ft. 

9.6 

15.3  lbs.  per  sq.  ft. 
+3.8  /  -1.5  Gs 


'  Information  is  from  the  former  Stoddard-Hamilton  company’s  website.  Performance  data  is  from 
company-built  aircraft. 


Tricycle  (on  gear) 

9.1  ft. 

Tricycle  (wings  folded) 

7.1  ft. 

Taildragger 

6.9  ft. 

Cabin  Dimensions 

Width  (at  hips) 

44.0  in. 

Width  (at  shoulders) 

46  in. 

Door  width 

37.0  in. 

Door  height 

31.5  in. 

Baggage  space 

32.0  cu.  ft. 

Weights 

Maximum  gross  weight  (on  wheels) 

1,960  lbs. 

Maximum  gross  weight  (on  floats) 

2,100  lbs. 

Empty  weight  (typical) 

1,200  lbs. 

Useful  load  (typical) 

760  lbs. 

Full-fuel  payload  (standard  tanks) 

520  lbs. 

Maximum  baggage  capacity 

250  lbs. 

Fuel  Capacity  (usable) 

Fuel  cells 

40.0  gals. 

APPENDIX  B.  EES  MODELS 


B.Ol  “Beer  and  Soup”  Models 


(a)  With  To  Varied 
"Refrigeration  Cycle" 

"Compressor" 

m_dot* (h_VI-h_VII) +W_dot_V=0  " [kW] " 
m_dot* (h_VI-h_VII_s) +W_dot_V_s=0  " [kW] " 
m_dot* (s_VI-s_VII) +S_PI_V=0  " [kW/K] " 
m_dot* (s_VI-s_VII_s)=0  " [kW/K] " 
X_dot_VI-X_dot_VII+W_dot_V-X_dot_d_V=0  " [kW] " 
h_VI=enthalpy(rll4,  p=p_VI,  s=s_VI)  "[kJ/kg]" 
h_VII=enthalpy (rll4,  p=p_VII,  s=s_VII)  "[kJ/kg]" 
h_VII_s=enthalpy (rll4 ,  p=p_VII,  s=s_VII_s)  "[kJ/kg] 
X_dot_VI=m_dot* (h_VI-h_0-T_0* (s_VI-S_0) )  " [kW] " 
X_dot_VII=m_dot* (h_VII-h_0-T_0* (s_VII-S_0) )  " [kW] " 
eta_V=W_dot_V_s /W_dot_V 

"Condenser" 

m_dot* (h_VII-h_TI) -Q_dot_C=0  " [kW] " 
m_dot* (s_VII-s_TI) -S_tau_c+S_PI_C=0  " [kW/K] " 
X_dot_VlI-X_dot_C-X_dot_TI-X_dot_d_C=0  " [kW] " 
X_dot_TI=m_dot* (h_TI-h_0-T_0* (s_TI-S_0) )  " [kW] " 
h_TI=enthalpy (rll4,  p=p_TI,  s=s_TI)  "[kJ/kg]" 
S_tau_C=Q_dot_C/T_soup  " [kW/K] " 

X_dot_C= {T_soup-T_0) *S_tau_C  " [kW] " 

"Turbine" 

m_dot* (h_TI-h_TII ) -W_dot_T=0  " [kW] " 
m_dot* (h_TI-h_TII_s) -W_dot_T_s=0  " [kW] " 
m_dot* (s_TI-s_TII) +S_PI_T=0  " [kW/K] " 
m_dot* (s_TI-s_TII_s)=0  " [kW/K] " 
X_dot_TI-X_dot_TII-W_dot_T-X_dot_d_T=0  " [kW] " 
h_TII=enthalpy (rll4 ,  p=p_TII,  s=s_TII)  "[kJ/kg]" 
h_TII_s=enthalpy (rll4 ,  p=p_TII,  s=s_TII_s)  "[kJ/kg] 
X_dot_TII=m_dot* (h_TII-h_0-T_0* (s_TII-S_0) )  " [kW] " 
eta_T=W_dot_T/W_dot_T_s 

"Evaporator" 

in_dot *  ( h_TI  I -h_VI )  +Q_dot_E=0  "  [  kW]  " 
m_dot* (s_TII-s_VI) +S_tau_E+S_PI_E=0  " [kW/K] " 
X_dot_TII-X_dot_VI— X_dot_E-X_dot_d_E=0  " [kW] " 
S_tau_E=Q_dot_E/T_beer  " [kW/K] " 

X_dot_E=- (T_beer-T_0 ) *S_tau_E  " [kW] " 

"Equilibrium" 

s_0=entropy (rll4 ,  T=T_0,  x=0)  "[kJ/kg-K]" 
h_0=enthalpy (rll4,  T=T_0,  x=0)  "[kJ/kg]" 


c_p=4.2  "[kJ/kg-K]" 

m_dot_beer*c_p* (T_beer_I-T_beer_II) -Q_dot_E=0 
m_dot_soup*c_p* (T_soup_I-T_soup_II ) +Q_dot_C=0 
T_soup_I=300  "[K]" 

T_beer_I=288  "[K]" 

T_b  e  e  r= ( T_be  e  r_ I - T_be e  r_I I ) / 1 n ( T_be e  r_I / T_be  e  r_I I ) 
T_soup= (T_soup_I-T_soup_II ) /In (T_soup_I/T_soup_II ) 


"Boundary  Conditions" 

T_soup=323  "[K]" 

T_beer=282  "IK]" 

p_VII=pressure (rll4,  T=T_soup+27, 
p_TII=pressure (rll4,  T=T_beer-27, 
p_TI=p_VII  "[kPa]" 
p_VI=p_TII  "[kPa]" 
T_sat_C=temperature (rll4,  p=p_VII, 
T_sat_E=temperature (rll4,  p=p_TII, 
s_VI=entropy (rll4 ,  p=p_VI,  T=T_VI) 
s_TI=entropy (rll4 ,  p=p_TI,  x=q_TI) 
q_TI=0. 3 

T_VI=T_sat_E+5  " [K] " 
eta_T=0 . 0 
eta  V=0.7 


x=0)  "[kPa]" 
x=0)  "[kPa]" 


x=0)  "[K]" 
x=0)  "[K]" 

" [kJ/kg-K] " 
" [kJ/kg-K] " 


{X_dot_E=l. 544*22  "[kW]"} 
m_dot=88.46  "[kg/s]" 


"Calculated" 

eta_II= (X_dot_E+X_dot_C) / (W_dot_V-W_dot_T) 
T_Vll=temperature (rll4,  p=p_VII,  s=s_VII)  "[K]" 
T_TII=temperature (rll4,  p=p_TII,  s=s_TII)  "[K]" 
q_TII=quality (rll4 ,  p=p_TII,  s=s_TII) 
q_VII=quality (rll4,  p=p_VII,  s=s_VII) 
W_dot_net=W_dot_V-W_dot_T  " [kW] " 
X_dot_VI_out=-X_dot_VI  "[kW]" 

X  dot_TII_in=-X_dot  Til  "[kW]" 


"Costing" 

Z_V+c_power*W_dot_V=c_VI *X_dot_VI_out+c_VI I *X_dot_VI I 
Z_C+c_VI I *X_dot_VI I=c_TI *X_dot_TI +c_X_C*X_dot_C 
Z_T+c_TI*X_dot_TI+c_TII*X_dot_TII_in=c_power*W_dot_T 
Z_E+c_VI *X_dot_VI_out=c_T I I *X_dot_TI I_in+c_X_E*X_dot_E 
{Z_V_A=25000* . 15* (m_dot* convert (kg/s, 

lb_m/s) ) * (p_VII/p_VI) ^0.45* (eta_V/ {l-eta_V) ) ^0. 65  " [$] " 
Z_T_A=25000*  .  32*  (in_dot* convert  (kg/s, 

lb_in/s)  )  *  (p_TI/p_TII)  ''0.5*  (eta_T/  (l-eta_T)  )  "0.85  "  [$]  " 
Z_C_A=15*15*Q_dot_C*15" (-1) *6" (- . 16) * . 005" (- . 04 ) 
Z_E_A=15*15*Q_dot_E*15" (-1) *6" (-.16) *.005" (-.04)  } 
Z_V_A=0  "[$]" 

Z_T_A=0  "[$]" 

Z_C_A=0 

Z_E_A=0 

Z  V=Z  V  A*. 1175/(365. 25*24*  3600)  "[$/s]" 


Z_T=Z_T_A*. 1175/ (365.25*24*  3600)  "[$/s]" 
Z_C=Z_C_A*. 1175/ (365. 25*24*  3600)  "[$/s]" 
Z_E=Z_E_A*. 1175/ (365.25*24*  3600)  "[$/s]" 
C_dot_beer=c_X_E*X_dot_E*convert (1/s, 1/hr)  " [$/hr] " 
C_dot_soup=c_X_C*X_dot_C*convert ( 1/s, 1/hr)  " [$/hr] " 

c_power=0. 07/3600  "[$/kJ]" 
c_VII=c_Vl 
c  TII=C  TI 


fij  With 

"Refrigeration  Cycle" 

"Compressor" 

m_dot* (h_VI-h_VII ) +W_dot_V=0  " [kW] " 
m_dot* (h_VI-h_VII_s) +W_dot_V_s=0  " [kW] " 
m_dot* (s_VI-s_VII) +S_PI_V=0  " [kW/K] " 
m_dot* (s_VI-s_VII_s)=0  " [kW/K] " 
X_dot_VI-X_dot_VII+W_dot_V-X_dot_d_V=0  " [kW] " 
h_VI=enthalpy(rll4,  p=p_VI,  s=s_VI)  "[kJ/kg]" 
h_VII=enthalpy (rll4,  p=p_VII,  s=s_VII)  "[kJ/kg]" 
h_VII_s=enthalpy (rll4 ,  p=p_VII,  s=s_VII_s)  " [kJ/kg] 
X_dot_VI=m_dot* (h_VI-h_0-T_0* (s_VI-S_0) )  " [kW] " 
X_dot_VII=m_dot* (h_VII-h_0-T_0* (s_VII-S_0) )  " [kW] " 
e t a_V=W_dot_V_s /W_dot_V 

"Condenser" 

m_dot* (h_VII-h_TI) -Q_dot_C=0  " [kW] " 
m_dot* (s_VII-s_TI) -S_tau_c+S_PI_C=0  " [kW/K] " 
X_dot_VII-X_dot_C-X_dot_TI-X_dot_d_C=0  " [kW] " 
X_dot_TI=m_dot*  (h_TI-h_0-T_0*  (s_^TI-S_0)  )  "  [kW]  " 
h_TI=enthalpy (rll4 ,  p=p_TI,  s=s_TI)  "[kJ/kg]" 
S_tau_C=Q_dot_C/T_soup  " [kW/K] " 

X_dot_C=  (T_soup-T_0)  *S_tau_C  '■  [kW]  " 

"Turbine" 

m_dot* (h_TI-h_TII) -W_dot_T=0  " [kW] " 
m_do t  * ( h_T I -h_T I I_s ) -W_do  t_T_s= 0  " [ kW ] " 
m_dot* (s_TI-s_TII) +S_PI_T=0  " [kW/K] " 
m_dot* (s_TI-s_TII_s)=0  " [kW/K] " 
X_dot_Tl-X_dot_TII-W_dot_T-X_dot_d_T=0  " [kW] " 
h_TII=enthalpy (rll4,  p=p_TII,  s=s_TII)  "[kJ/kg]" 
h_TII_s=enthalpy (rll4,  p=p_TlI,  s=s_TII_s)  "[kJ/kg] 
X_dot_TII=m_dot* (h_TII-h_0-T_0* (s_TII-S_0) )  " [kW] " 
e  t  a_T=w_do  t_T / W_do  t_T_s 

"Evaporator" 

m_dot* (h_TII-h_VI) +Q_dot_E=0  " [kW] " 

m_dot* (s_TII-s_VI) +S_tau_E+S_PI_E=0  " [kW/K] " 

X_dot_TII-X_dot_VI— X_dot_E-X_dot_d_E=0  " [kW] " 

S  tau  E=Q  dot  E/T  beer  "[kW/K]" 
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X_dot_E=- (T_beer-T_0 ) *S_tau_E  " [kW] " 

"Equilibrium" 

s_0=entropy (rll4 ,  T=T_0,  x=0)  "[kJ/kg-K]" 
h_0=enthalpy (rll4,  T=T_0,  x=0)  "[kJ/kg]" 
c_p=4.2  "[kJ/kg-K]" 

m_dot_beer*c_p* (T_beer_I-T_beer_II ) -Q_dot_E=0 
m_dot_soup*c_p* (T_soup_I-T_soup_II) +Q_dot_C=0 
T_soup_I=300  "[K]" 

T_beer_I=288  " [K] " 

T_beer= (T_beer_I-T_beer_II) /In (T_beer_I/T_beer_II ) 
T_soup= (T_saup_I-T_soup_II) /In {T_soup_I/T_soup_II ) 
T_0_I=T_soup_I* (T_soup_I/T_beer_I) ^ ( {- 
m_dot_beer*c_p) / (m_dot_beer*c_p+m_dot_soup*c_p) )  " [K] " 
T_0_II=T_soup_II*  (T_soup_II/T_beer_II ) ''  (  (- 
m_dot_beer*c_p) / (m_dot_beer*c_p+m_dot_soup*c_p) )  " [K] " 
T_0==  ( T_0_I -T_0_I I )  / In  ( T_0_I /T_0_I I ) 

"Boundary  Conditions" 

T_soup=323  "[K]" 

T_beer=282  " [K] " 

p_VII=pressure (rll4,  T=T_soup+27,  x=0)  "[kPa]" 
p_TII=pressure (rll4,  T=T_beer-27,  x=0)  "[kPa]" 
p_TI=p_VII  "[kPa]" 
p_VI=p_TIl  "[kPa]" 

T_sat_C=temperature (rll4 ,  p=p_VII,  x=0)  " [K] " 
T_sat_E=temperature (rll4,  p=p_TII,  x=0)  " [K] " 
s_VI=entropy (rll4,  p=p_VI,  T=T_VI)  "[kJ/kg-K]" 
s_TI=entropy (rll4,  p=p_TI,  x=q_TI)  "[kJ/kg-K]" 
q_TI=0.3 

T_VI=T_sat_E+5  " [K] " 
eta_T=0 . 0 
eta_V=0.7 

m_dot=88.46  "[kg/s]" 

"Calculated" 

eta_II= (X_dot_E+X_dot_C) / (W_dot_V-W_dot_T) 
T_VII=temperature (rll4,  p=p_VII,  s=s_VII)  " [K] " 
T_TII=temperature (rll4,  p=p_TII,  s=s_TII)  " [K] " 
q_TII=quality (rll4,  p=p_TII,  s=s_TII) 
q_VII=quality (rll4,  p=p_VII,  s=s_VII) 
W_dot_net=W_dot_V-W_dot_T  " [kW] " 

X_dot_VI_out=-X_dot_VI  "[kW]" 

X_dot_TII_in=-X_dot_TII  "[kW]" 

"Costing" 

Z_V+c_power*W_dot_V=c_Vl*X_dot_VI_out+c_VII*X_dot_VlI 
Z_C+c_VI I *  X_do  t_VI I = c_T I *  X_do  t_T I +c_X_C  *  X_do t_C 
Z_T + c_T I *  X_do t_T I + c_T 1 1 *  X_do t_T 1 1 _i n= c_po we  r * W_do  t _T 
Z_E+c_VI*X_dot_VI_out=c_TII*X_dot_Tll_in+c_X_E*X_dot_E 

Z  V  A=0  " [$] " 
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Z_T_A=0  "[$]" 

Z_C_A=0 

Z_E_A=0 

Z_V=Z_V_A*.1175/ (365.25*24*  3600) 
Z_T=Z_T_A*. 1175/ (365,25*24*  3600) 
Z_C=Z_C_A*. 1175/ (365.25*24*  3600) 
Z  E=Z  E  A*. 1175/(365. 25*24*  3600) 


"[$/s]" 

"[$/s]" 

"[$/s]" 

"[$/s]" 


C  dot  beer=c  X  E*X  dot  E*convert (1/s, 1/hr) 


C_dot_soup=c_X_C*X_dot_C*convert (1/s, 1/hr) 


c_power=0. 07/3600  "[$/kJ]" 


c  VII=c  VI 


c  TII=c  TI 


" [$/hr] " 
"[$/hr]" 


B.02  Aircraft  Models 


(a)  Cruise,  75%  Power,  8000ft. 

"Light  A/C  Exergy  Balance  and  Costing,  Cruise" 

"Exergy  Balances" 

"Fuel:"  X_dot_L_F-X_dot_d_F=0  "[kW]" 

"Engine"  X_dot_Ch_F-X_dot_P_P-X_dot_P_A-X_dot_Exh+X_dot_L_E- 
X_dot_d_E=0  "[kW]" 

"Propeller"  X_dot_P_P+X_dot_L_P-X_dot_T-X_dot_d_P=0  "[kW]" 

"  A1 1  e  mat  or "  X_dot_P_A+X_do t_L_A-X_dot_e 1 e  ct -X_dot_d_A=  0 
"  [kW] " 

"Wing"  X_dot_T_W-X_dot_L_net-X_dot_d_W=0  "[kW]" 

"Fuselage  etc.:"  X_dot_T_M+X_dot_L_M-X_dot_d_M=0  "[kW]" 
"Cargo:"  X_dot_L_C-X_dot_d_C=0  "[kW]" 

"Auxiliary  Relations" 

X_dot_L_net=X_dot_L-X_dot_L_W 

X_dot_L_F=m_F/m_G*X_dot_L 

X_do t_L_E=m_E /m_G  *  X_dot_L 

X_do  t  _L_P=in_P  /m_G  *  X_do  t_L 

X_do  t_L_A=m_A/m_G  *  X_do t_L 

X_dot_L_W=m_W/m_G*X_dot_L 

X_do  t_L_M=m_M/m_G  *  X_do t_L 

X_dot_L_C=m_C  /in_G  *X_dot_L 

X_dot_T=X_dot_T_W+X_dot_T_M 

X_dot_T=eta_P*X_dot_P_P 

x_dot_elect=eta_A*X_dot_P_A 

X_dot_elect=volt*i*convert (W,  kW) 

X_dot_L=m_G*g* (RC-RC_0) *convert (W,  kW) 
X_dot_T_W=V*rho_air*V''2/2*C_D*A_wing*convert  (W,  kW) 

X_do  t_Ch_F=m_d  o  t_F  *  X_Ch_F 

X_dot_T_M=V*rho_air*V''2/2*A_equiv_M*convert  (W,  kW) 
P_req=X_dot_T_W+X_dot_T_M  " [kW] " 

P_ex=X_dot_T-P_req  "[kW]" 

RC_0=-rho_air*V^3*C_L^2*A_wing/ (2*pi*a) / (m_G*g) 


"Boundary  Conditions" 
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"Masses" 

m_F=120*convert {lb_m,  kg)  "[kg]" 
m_E=2 60* convert {lb_m,  kg)  "[kg]" 

m_P=4 0* convert (lb_m,  kg)  "[kg]"  "Includes  propeller  & 
governor" 

m_A=6.  l*convert  (lb_in,  kg)  "[kg]" 
m_W=  4  0  0  *  con ve  r t ( lb_m,  kg )  " [ kg ] " 
ni_C=450*convert  {lb_m,  kg)  "[kg]" 
m_M=520*convert (lb_m,  kg)  "[kg]" 
m  G-m  F-m  E-m  P-m  A-m  W-m  C-m_M=0  "[kg]" 


"A/C  Geometry" 

A_wing=128*convert  (ft''2,  m^2)  "[m''2]" 

"Aspect  Ratio:  "  a=9.6 

"Equipment  Performance" 

eta_P=0.814 

eta_A=0 . 8 

"A/C  Performance" 

"Speed:  "  {V=161*convert (mph,  m/s)}  "[m/s]" 

"Engine  Power:"  X_dot_P_A+X_dot_P_P=120* convert (hp,  kW)  "[kW]" 
"Fuel  Flow:"  m_dot_F=8 . 4*6*convert (lb_m/hr,  kg/s)  "[kg/s]"  . 
"Current  Draw:  "  i=30  " [amp] " 

"Voltage:"  volt=14  "[volt]" 

{C_D=0.0075} 
k=0.0038 
C_D_min=0.007 
e=l/  (H-delta+k*pi*a) 

C_D=C_D_min+C_L^2/ (pi*a*e) 
delta=0.17 

C_L=C_l_alpha* (a/ {a+ (2* {a+4 ) / (a+2 ) ) ) ) 

A_equi v_M=0 .2998 

V_mph=V*convert (m/s,  mph)  "[mph]" 

RC=0 


"Mi sc" 

X_Ch_F=46000  "[kJ/kg]" 

X_dot_Exh=0  "[kW]"  "Lump  in  with  engine  destruction" 
rho_air=0 . 75*density (air,  p=101,  T=288) 
rho_air*V''2/2*C_L*A_wing=m_G*g 
g=9.8  "[m/s"2]" 


"Costing" 

"Fuel : "  c_F_total*X_dot_Ch_F=c_F*m_dot_F+c_Li*X_dot_L_F 
Z_dot_F=c_F*m_dot_F*convert (1/s, 1/hr) 
Z_dot_F_total=c_F_total*X_dot_Ch_F* convert (1/s, 1/hr) 
"Engine"  c_F_total*X_dot_Ch_F-c_P*X_dot_P_P- 
c_P*X_dot_P_A+c_Li*X_dot_L_E+Z_E=0 
Z_dot_P_P=c_P*X_dot_P_P* convert (1/s, 1/hr) 

Z  dot  P  A=c  P*X  dot  P  A*convert (1/s, 1/hr) 


Z_dot_L_E=c_Li*X_dot_L_E* convert (1/s, 1/hr) 

"Propeller"  c_P*X_dot_P_P+c_Li*X_dot_L_P-c_T*X_dot_T+Z_P=0 
Z_dot_T=c_T*X_dot_T* convert (1/s, 1/hr) 
Z_dot_L_P=c_Li*X_dot_L_P* convert (l/s, 1/hr) 

"Alternator"  c_P*  X_do  t_P_A+c_Li *  X_do t_L_A- 
c_elect*X_dot_elect+Z_A=0 
Z_dot_L_A=c_Li*X_dot_L_A* convert (1/s, 1/hr) 
Z_dot_elect=c_elect*X_dot_elect*convert (1/s, 1/hr) 

"Wing"  c_T*X_dot_T_W-c_Li*X_dot_L_net+Z_W=0 
Z_dot_T_W=c_T*X_dot_T_W*convert (1/s, 1/hr) 

"Fuselage  etc.:"  c_T*X_dot_T_M+c_Li*X_dot_L_M+Z_M=C_dot_M 
Z_dot_T_M=c_T*X_dot_T_M*convert (1/s, 1/hr) 

Z_dot_T_M=Z_S 

Z_dot_L_M=c_Li*X_dot_L_M*convert (1/s, 1/hr) 

"Cargo : "  c_Li *X_dot_L_C=C_dot_C 
Z_dot_L_C=c_Li*X_dot_L_C*convert (l/s, 1/hr) 

"Capital  Costs" 
c_T=8 . le-6 
Z_E=0 
Z_P=0 

Z_A=4.159e-5 

Z_W=0 

Z_M=0 

(b)  Cruise,  65%  Power,  8000 ft. 

"Light  A/C  Exergy  Balance  and  Costing,  Cruise" 

"Exergy  Balances" 

"Fuel:"  X_dot_L_F-X_dot_d_F=0  "[kW]" 

"Engine"  X_dot_Ch_F-X_dot_P_P-X_dot_P_A-X_dot_Exh+X_dot_L_E 
X_dot_d_E=0  "[kW]" 

"Propeller"  X_aot_P_P+X_dot_L_P-X_dot_T-X_dot__d_P=0  "  [kW.]  " 
"Alternator"  X_do t_P_A+X_do t_L_A-X_do  t_e 1 e  c t -X_do  t_d_A=0 
"  [kW]  " 

"Wing"  X_dot_T_W-X_dot_L_net-X_dot_d_W=0  "[kW]" 

"Fuselage  etc.:"  X_dot_T_M+X_dot_L_M-X_dot_d_M=0  "[kW]" 
"Cargo: "  X_dot_L_C-X_dot_d_C=0  " [kW] " 

"Auxiliary  Relations" 

X_dot_L_net=X_dot_L-X_dot_L_w 

X_dot_L_F=m_F/m_G*X_dot_L 

X_do  t_L_E=m_E  /iia_G  *  X_do  t_L 

X_do  t_L_P=m_P /m_G  *  X_do t_L 

X_dot_L_A=m_A/m_G*X_dot_L 

X_dot_L_W=m_W/m_G*X_dot_L 

X_dot_L_M=m_M/m_G*X_dot_L 

X_do  t_L_C=m_C  /in_G  *  X_do  t_L 

X_dot_T=X_dot_T_W+X_dot_T_M 

X_dot_T=eta_P*X_dot_P_P 

X_dot_elect=eta_A*X_dot_P_A 

X  dot  elect=volt*i*convert (W,  kW) 
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X_dot_L=m_G*g* (RC-RC_0) ^convert (W,  kW) 
X_dot_T_W=V*rho_air*V"2/2*C_D*A_wing*convert (W, kW) 
X_dot_Ch_F=in_dot_F*X_Ch_F 

X_dot_T_M=V*rho_air*V^2/2*A_equiv_M*convert (W, kW) 
P_req=X_dot_T_W+X_dot_T_M  " [kW] " 
P_ex=X_dot_T-P_req  "[kW]" 

RC_0=-rho_air*V"3*C_L"2*A_wing/ (2*pi*a) / (m_G*g) 

"Boundary  Conditions" 

"Masses" 


m_F=12 0* convert (lb_m,  kg) 

"[kg]" 

m_E=2 60* convert (lb_m,  kg) 

" [kg] " 

m  P=4 0* convert (lb_m,  kg) 

"[kg]"  "Includes  propeller  & 

governor" 

m_A=l 3* convert (lb_m,  kg) 

" [kg] " 

m  W=400*convert (lb  m,  kg) 

" [kg] " 

m_C=220.9  "[kg]" 

m_M=520*convert (lb_m,  kg) 

" [kg] " 

m  G-m_F-in_E-m_P-m_A-m_W-m_C-in_M=0  "[kg]" 


"A/C  Geometry" 

A_wing=128*convert  (ft^2,  m'^2)  "[m''2]" 

"Aspect  Ratio:  "  a=9.6 

"Equipment  Performance" 

eta_P=0.814 

eta_A=0 . 8 

"A/C  Performance" 

"Speed:  "{  V=161*convert (mph,  m/s)}  "[m/s]" 

"Engine  Power:"  X_dot_P_A+X_dot_P_P=104*convert (hp,  kW)  "[kW]" 
"Fuel  Flow:"  m_dot_F=6 . 9*6*convert (lb_m/hr,  kg/s)  "[kg/s]" 
"Current  Draw:  "  i=30  " [amp] " 

"Voltage: "  volt=14  " [volt] " 

k=0.0038 

C_D_min=0 ,007 

e=l/ {l+delta+k*pi*a) 

C_D=C_D_min+C_L''2/  (pi*a*e) 
delta=0.17 

C_L=C_l_alpha* (a/ (a+ (2* (a+4 ) /  (a+2) ) ) ) 

A_equi v_M=0 .2988 

V_mph=V*convert (m/s,  mph)  "[mph]" 

RC=0 


"Mi sc" 

X_Ch_F=46000  "[kJ/kg]" 

X_dot_Exh=0  "[kW]"  "Lump  in  with  engine  destruction" 
rho_air=0 . 75*density (air,  p=101,  T=288) 
rho_a  ir*V'^2/2*  C_L*  A_wing=m_G*  g 
g=9.8  "[m/s^2]" 

Ratio_DL= (X_dot_T*convert (kW,  W) /V) / (m_G*g) 
BSFC=m_dot_F/X_dot_P_  P  " [kg/kW-s] " 


"Costing" 

"Fuel : "  c_F_total*X_dot_Ch_F=c_F*m_dot_F+c_Li*X_dot_L_F 
Z_dot_F=c_F*m_dot_F* convert (1/s, 1/hr) 
Z_dot_F_total=c_F_total*X_dot_Ch_F*convert (1/s, 1/hr) 
"Engine"  c_F_total*X_dot_Ch_F-c_P*X_dot_P_P- 
C_P*X_dot_P_A+c_Li*X_dot_L_E+Z_E=0 
Z_dot_P_P=c_P'*X_dot_P_P*convert  (1/s,  1/hr) 
Z_dot_P_A=c_P*X_dot_P_A* convert (1/s, 1/hr) 
Z_dot_L_E=c_Li*X_dot_L_E*convert (1/s, 1/hr) 

"Propeller"  c_P*X_dot_P_P+c_Li*X_dot_L_P-c_T*X_dot_T+Z_P=0 
Z_dot_T=c_T*X_dot_T* convert (1/s, 1/hr) 
Z_dot_L_P=c_Li*X_dot_L_P*convert (1/s,  1/hr) 

"Alternator"  c_P  *  X_do  t_P_A+  c_Li *  X_do t_L_A- 
c_elect*X_dot_elect+Z_A=0 
Z_dot_L_A=c_Li*X_dot_L_A*convert (1/s, 1/hr) 
Z_dot_elect=c_elect*X_dot_elect*convert (1/s, 1/hr) 

"Wing"  c_T  *  X_do  t_T_W - c_Li *  X_dot_L_ne t + Z_W=  0 
Z_dot_T_W=c_T*X_dot_T_W*convert (1/s, 1/hr) 

"Fuselage  etc.:"  c_T*X_dot_T_M+c_Li*X_dot_L_M+Z_M=C_dot_M 
Z_dot_T_M=c_T*X_dot_T_M* convert (1/s, 1/hr) 

Z_dot_T_M=Z_S 

Z_dot_L_M=c_Li*X_dot_L_M* convert (1/s, 1/hr) 

"Cargo: "  c_Li*X_dot_L_C=C_dot_C 
Z_dot_L_C=c_Li*X_dot_L_C* convert (1/s, 1/hr) 

"Capital  Costs" 
c_F=(2+. 4 14) *1/6*2.2046 
Z_E=0 
Z_P=0 

Z_A=4.159e-5 

Z_W=0 

Z_M=0 

(c)  Climb,  Sea  Level  (Maximize  RC  as  a  function  of  V) 

"Light  A/C  Exergy  Balance  and  Costing,  Climb" 

"Exergy  Balances" 

"Fuel:"  X_dot_L_F-X_dot_d_F=m_F*g*RC*convert (W,  kW)  "[kW]" 
"Engine "  X_dot_Ch_F-X_dot_P_P-X_dot_P_A-X_dot_Exh+X_dot_L_E 
X_do t_d_E=m_E  *  g  *  RC  *  con ve  r t ( W ,  kW )  " [ kW ] " 

"Propeller"  X_dot_P_P+X_dot_L_P-X_dot_T- 
X_dot_d_P=m_P* g*RC* convert (W,  kW)  "[kW]" 

"Alternator"  X_do  t_P_A+X_do  t_L_A-X_do  t_e 1 e  c t - 
X_do t_d_A=m_A*  g  *  RC  *  conve  r t ( W ,  kW )  " [ kW ] " 

"Wing"  X_dot_T_W-X_dot_L_net- 

X_dot_d_W+P_ex=m_W*g*RC*convert (W,  kW)  "[kW]" 

"Fuselage  etc.:"  X_dot_T_M+X_dot_L_M- 

X_dot_d_M=m_M*g*RC*convert (W,  kW)  "[kW]" 

"Cargo:"  X_dot_L_C-X_dot_d_C=m_C*g*RC*convert (W,  kW)  "[kW] 


"Auxiliary  Relations" 
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X_clot_L_net=X_dot_L-X_dot_L_W 
X_dot_L_F=m_F/in_G*X_dot_L 
X_do  t_L_E=n:i_E  /xn_G  *  X_do  t_L 
X_do  t_L_P=m_P /m_G  *  X_do  t_L 
X_do  t_L_A=m_A/m_G  *  X_do t_L 
X_dot_L_W=m_W/m_G*X_dot_L 
X_do t_L_M=m_M/m_G  *  X_do  t_L 
X_do  t_L_C=m_C /m_G  *  X_do t_L 

X_dot_T=X_dot_T_W+X_dot_T_M+m_G*g*RC*Gonvert (W,  kW) 

X_do  t_T=e  t  a_P  *  X_do  t_P_P 
X_do  t_e 1 e  ct =e  t  a_A*  X_do  t_P_A 
X_dot_elect=volt*i*convert (W,  kW) 

X_dot_L=m_G*g* (RC-RC_0) *convert (W,  kW) 
X_dot_T_W=V*rho_air*V^2/2*C_D*A_wing*convert (W, kW) 

X_do t_Ch_F=m_do  t_F  *  X_Ch_F 

X_dot_T_M=V*rho_air*V''2/2*A_equiv_M*convert  (W,  kW) 
P_req=X_dot_T_W+X_dot_T_M  " [kW] " 

P_ex=X_dot_T-P_req  "[kW]” 

RC_0=-rho_air*V''3*C_L''2*A_wing/  {2*pi*a)  /  (m_G*g) 
P_req=X_dot_T_W+X_dot_T_M  " [kW] " 

P_ex=X_dot_T-P_req  "[kW]" 

"Boundary  Conditions" 

"Masses" 

m_G=1960*convert (lb_m,  kg)  "[kg]"  "Gross  weight" 
m_F=240*convert (lb_m,  kg)  "[kg]" 
m_E=275*convert {lb_m,  kg)  " [kg] " 

m_P=40*convert (lb_m,  kg)  "[kg]"  "Includes  propeller  & 
governor" 

m_A=13*convert (lb_m,  kg)  "[kg]" 
m_W=4 00* convert (lb_m,  kg)  "[kg]" 
m_G-m_F-m_E-m_P-m_A-m_W-m_C=m_M  " [kg] " 
m_M=235.9  "[kg]" 

"A/C  Geometry" 

A_wing=128*convert (ft^2,  m^2)  "[m^2]" 

"Aspect  Ratio:  "  a=9.6 

"Equipment  Performance" 

eta_P=0. 197+8. 25e-3* (V_mph/1 . 15) -2 . 723e-5* {V_mph/1.15) ^2 
eta_A=0 . 8 

"A/C  Performance" 

"Speed:  "  V_mph=V* convert (m/s,  mph)  " [mph] " 

"Engine  Power:"  X_dot_P_A+X_dot_P_P=l 60* convert (hp,  kW)  "[kW]" 
"Fuel  Flow:"  m_dot_F=16*6*convert (lb_m/hr,  kg/s)  "[kg/s]" 
"Current  Draw:  "  i=30  " [amp] " 

"Voltage:"  volt=14  "[volt]" 

"Climb  Rate"  RC_fpm=RC*convert (m/s,  ft/min)  "[ft/min]" 

k=0.0038 

C_D_min=0.007 

e=l/ (l+delta+k*pi*a) 


C_D=C_D_min+C_L"'2/  (pi*a*e) 
delta=0 . 17 

C_L=C_l_alpha* (a/ (a+(2* (a+4) / (a+2) ) ) ) 

A_equiv_M=0 . 2998  ''[m"'2]" 

Theta_C=arctan (RC/V) 

"Mi sc” 

X_Ch_F=46000  "[kJ/kg]" 

X_dot_Exh=0  "[kW]"  "Lump  in  with  engine  destruction" 
rho_air=density (air,  p=101,  T=288) 
r  ho_a  ir*V''2/2*  C_L  *  A_wi  ng=m_G  *  g 
g=9.8  "[m/s "2]" 


"Costing" 

"Fuel : "  c_F_total*X_dot_Ch_F=c_F*m_dot_F+c_Li*X_dot_L_F 
Z_dot_F=c_F*m_dot_F* convert (1/s, 1/hr) 
Z_dot_F_total=c_F_total*X_dot_Ch_F*convert (1/s, 1/hr) 
"Engine"  c_F_total*X_dot_Ch_F-c_P*X_dot_P_P- 
c_P*X_dot_P_A+c_Li*X_dot_L_E+Z_E=0 
Z_dot_P_P=c_P*X_dot_P_P* convert (1/s, 1/hr) 
Z_dot_P_A=c_P*X_dot_P_A* convert (1/s, 1/hr) 
Z_dot_L_E=c_Li*X_dot_L_E* convert (1/s, 1/hr) 

"Propeller"  c_P*X_dot_P_P+c_Li*X_dot_L_P-c_T*X_dot_T+Z_P=0 
Z_dot_T=c_T*X_dot_T* convert (1/s, 1/hr) 
Z_dot_L_P=c_Li*X_dot_L_P*convert (1/s, 1/hr) 

"Alternator"  c_P*X_dot_P_A+c_Li*X_dot_L_A- 
c_elect*X_dot_elect+Z_A=0 
Z_dot_L_A=c_Li*X_dot_L_A* convert (1/s, 1/hr) 
Z_dot_elect=c_elect*X_dot_elect*convert (1/s, 1/hr) 

"Wing"  c_T  *  X_d o  t _T_W - c_L i *  X_do  t_L_n e  t + Z_W=  0 
Z_dot_T_W=c_T*X_dot_T_W* convert (1/s, 1/hr) 

"Fuselage  etc.:"  c_T*X_dot_T_M+c_Li*X_dot_L_M+Z_M=C_dot_M 
Z_dot_T_M=c_T*X_dot_T_M*convert (1/s, 1/hr) 

Z_dot_L_M=c_Li *X_dot_L_M* convert (1/s, 1/hr) 

"Cargo: "  c_Li*X_dot_L_C=C_dot_C 
Z_dot_L_C=c_Li*X_dot_L_C*convert (1/s, 1/hr) 

"Capital  Costs" 

Z_E=0 

Z_P=0 

Z_A=0 

Z_W=0 

Z_M=0 

c  T=6.944e-6  " [dol/kJ] " 


53 


APPENDIX  C.  YEAR  2001  PUBLICATIONS 

Paulus,  D.M.  Jr,  and  R.A.  Gaggioli,  2001,  “The  Exergy  of  Lift,  and  Aircraft  Exergy  Flow  Diagrams”, 
ASME  AES-Vol  41,  ASME,  New  York  (to  be  submitted  to  ASME  Journal  of  Energy  Resource 
Technology) 

Paulus,  D.M.  Jr.  and  R.A.  Gaggioli,  2001,  “Rational  Design  of  Vehicles:  I.  Multi-Objective 
Optimization”,  ASME  AES-Vol  41,  ASME,  New  York 

Paulus,  D.M,  Jr.  and  R.A.  Gaggioli,  2001,  “Rational  Design  of  Vehicles:  II.  Subsystem 
Decomposition”,  ASME  AES-Vol  41,  ASME,  New  York 


54 


THE  EXERGY  OF  LIFT,  AND  AIRCRAFT  EXERGY  FLOW  DIAGRAMS 

David  M.  Paul  us,  Jr.  and  Richard  A.  Gaggioli 

Marquette  University 
Milwaukee,  W1 


ABSTRACT 

Aside  from  incidental,  auxiliary'  loads,  in  level 
flight  the  principal  load  on  the  aircraft  propulsion 
engine  is  the  power  required  to  provide  the  continuous 
lift.  To  construct  an  exergy  flow  diagram  for  an 
aircraft  -  for  example,  for  the  purpose  of  exergy’ 
costing  -  an  expression  is  needed  for  the  exergy 
delivered  to  {and  by)  the  wings.  That  is,  an  expression 
is  needed  for  the  exergy  of  lift.  In  this  paper  an 
expression  is  developed  for  this  quantity,  applicable 
not  only  in  level  flight  but  in  other  modes  of  flight  as 
well.  Then,  three  exergy  flow  diagrams  are  presented 
for  a  light  aircraft,  operating  under  three  conditions: 
climb,  level  economy  cruise,  and  level  maximum 
power.  The  expression  for  lift  exergy  and  these 
exergy  flow  diagrams  are  relevant  to  the  optimal 
design  of  aircraft  subsystems  [1], 


NOMENCLATURE 


Symbol 

Description 

a 

Aspect  ratio 

A 

Area 

c. 

Coefficient  of  drag  for  an  infinite  span 

Cd 

Coefficient  of  drag  for  a  finite  span 

c, 

Coefficient  of  lift  for  an  infinite  span 

c. 

Coefficient  of  lift  for  a  finite  span 

e 

Oswald’s  efficiency  factor 

F 

Force 

g 

Acceleration  of  gravity 

P 

Power 

V 

Velocity 

X 

Exergy 

X 

Specific  exergy 

Y 

Altitude 

Table  1:  Symbols 


Greek 

Symbol 

Description 

5 

Constant  for  determining  Oswald’s 
efficiency  factor 

V 

Efficiency 

^ _ 

Air  density 

Table  2:  Greek  Symbols 


Subscript 

Description 

0 

Dead  state 

d 

Destruction 

X 

x-component 

y  : 

y-component 

Table  3:  Subscripts 


INTRODUCTION 

The  concept  of  exerg>%  a  special  case  of  Gibbs' 
available  energy  for  a  body  and  large  medium,  is  a 
valuable  tool  in  both  the  analysis,  and  the  optimization 
of  energy  conversion  systems.  Its  value  in  the 
optimization  of  energy’  systems  lies  in  their 
decomposition,  breaking  a  large  system  into 
subsystems  and  devices  to  be  optimized  separately  [2], 
typically  with  the  objective  of  minimizing  cost. 

Aircraft  energy  systems,  however,  have 
traditionally  been  optimized  so  as  to  minimize  gross 
takeoff  weight.  Although  this  objective  does  reflect 
aircraft  lifecycle  cost  as  well  as  performance,  it  is  not 
necessarily  the  objective  variable  of  choice  for 
optimization,  nor  does  it  lend  itself  to  decomposition\ 
Recently  there  has  been  interest  in  applying  second- 
law  methodologies,  including  exergy’  analysis, 
thermoeconomics  and  entropy'  generation 
minimization  [3],  to  aircraft  energy  systems^. 

Aircraft  energy  systems  are  unique,  in  that  exergy^ 
is  required  not  only  to  operate  them  but  also  to  lift 
them  and  to  hold  them  aloft.  In  order  to  complete  an 
exergy  analysis,  or  apply  therraoeconomics,  it  is 
necessary^  to  create  exergy  flow  diagrams.  Required 
for  this  task  is  an  expression  for  the  exergy  of  lift. 

Traditionally,  the  exergy  associated  with  the 
application  of  a  force  has  been  given  as 

X;.=F.V 

This  is  in  fact  a  special  case;  the  general  expression  is 
actually 

x^=r-(v-v,) 

It  is  only  because  the  ‘‘dead  state”  velocity  is  typically 
taken  to  be  zero  that  (1)  is  then  correct.  However,  for 
aircraft,  equation  (1)  incorrectly  leads  to  a  conclusion 
that  the  exergy  transport  associated  with  lift  is  zero  in 
level  flight;  such  a  conclusion  does  not  charge  an  on¬ 
board  energy  system  for  the  exergy  it  consumes  in 
“staying  aloft”. 


*  For  a  discussion  of  objective  functions  for  vehicles,  see 
Paulus  and  Gaggioli  [1] 

For  a  general  overview  of  second-law  methodologies,  see 
Bejan  et  al.  [4] 
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Finding  a  proper  Fyo,  or  dead  state  velocity  in  the 
vertical  direction,  allows  a  proper  calculation  of  the 
exergy  of  lift. 

FINDING  THE  “DEAD  STATE”  VELOCITY 

McCormick  [5]  states  that  the  minimum  induced 
drag  coefficient  of  a  wing  is 


c 


=  ^1/ 


Dimn  (3) 

where  a  is  the  aspect  ratio.  For  drag  to  be  the 
minimiim  while  producing  a  given  amount  of  lift,  no 
parasitic  drag  may  be  present,  and  the  induced  drag 
coefficient  must  be  given  by  (3).  This  minimum  drag 
force  is  given  by 

pAVlCl 

Ina  (4) 

The  minimum  exergy  input  to  the  wing,  assuming 
Q  =  0  and  steady  flight,  is  therefore 

pAVlCl 

Ina  (5) 

An  exergy  balance  on  this  wing  in  level  flight  results 


X, 


thrust. min 


dX, 


^thrustjn  ^ d 


0: 

(6) 

The  exergy  destruction  for  the  minimum  thrust  exergy 
input  must  be  zero.  Using  this,  and  with  the 

substitution  of^  (2)  for  X^jj  and  (5)  for 

yields  the  following  expression  for  level  flight  (where 
and  the  lift  force  equals  Wag,  v^here  nta  is  the 
mass  of  the  aircraft): 

yS)  ^ 

2tCCI  (7) 

The  dead  state  velocity  is  then  found  to  equal 

Y  ^  ^  '^^aS 

iTtamg  TtapAV^ 

With  (8),  the  expression  for  the  steady  state  exergy^  of 
lift  for  any  component  of  mass  ntc  is 


( 


=  ^cg 


^  napAV^ 


\ 


xj 


(9) 


A  NOTE  ON  THE  DEAD  STATE  VELOCITY 

The  rate  of  climb  in  a  light  aircraft  is  well 
approximated  [5]  with 

P-P 


K  =  - 


req 


mg 


(10) 


^  in  the  y-direction 


where  P  is  the  power  supplied,  and  Preq  is  the  power 
required  for  level  flight. 

If  the  power  supplied  is  zero,  and  the  aircraft  is 
ideal"^, 

y  P^K’Cl 

Inamg 

This  is  noted  to  be  the  same  as  Vyo  given  b\'  equation 
(8).  Thus,  at  a  given  forward  velocit>’  the  dead 
state  velocity  corresponds  to  the  speed  of  descent  of 
an  ideal  wing,  at  zero  power  input,  that  is  to  sa>’,  in  a 
steady-state  glide. 

CREATION  OF  EXERGY  FLOW 
DIAGRAMS 

The  relationship  for  the  exergy  of  lift  will 
presently  be  applied  to  a  light  aircraft,  the  Glastar™ 
homebuilt  aircraft,  developed  b\^  Arlington  Aircraft 
Development,  Inc.  A  Lycoming  0-320  or  0-360 
engine  producing  between  112  and  149  kW  (150  to 
200  hp),  turning  either  a  constant-speed  or  fixed  pitch 
propeller,  typically  powers  this  aircraft.  The  aircraft 
has  a  gross  weight  of  889  kg  (1960  pounds).  The 
Glastar  has  a  NASA  GAW-2  airfoil,  a  wingspan  of 
(10.7  m)  35  feet,  and  an  aspect  ration  of  9.6.  For  the 
purpose  of  this  article,  the  aircraft  uall  be  assumed  to 
be  powered  by  a  119  kW  (160  hp)  0-320  engine  with 
a  constant-speed  Hartzell  propeller  and  to  have  an 
empty  weight  of  544  kg  (1200  pounds).^ 

Three  exergy  flow  diagrams  will  be  presented, 
one  for  full-throttle  cruise  (75%  power)  at  an  altitude 
of  2440  m  (8000  feet),  one  for  65%  power  cruise  at 
2440  m,  and  one  for  maximum  rate  of  climb  at  sea 
level.. 

The  aircraft  was  divided  into  the  following 
components:  Engine,  propeller,  alternator,  wing, 
fuselage  and  empennage  (the  horizontal  and  vertical 
stabilizers),  and  cargo.  Tnm  drag  and  non-wing 
parasitic  drag  were  attrib^ui  to  the  fuselage  and 
empennage  component. 


The  Aircraft  Flight  Model 

The  aircraft  performance  model  employed  is 
vastly  simplified.  The  wing  performance  was 
assumed  to  follow  the  equation  [5] 

r  =  r  +  .^..L 

me  (12) 

where  e,  Oswald’s  efficiency  factor  is  given  by 


1 


1  +  (5  +  km 


(13) 


^  Here,  imagine  an  aircraft  producing  no  parasitic  or  trim 
drag. 

^  As  these  aircraft  are  amateur-built  from  a  basic  airframe  kit, 
each  individual  aircraft  is  unique.  Variations  in  power  plant 
and  propeller  selection,  avionic  selection,  construction 
technic^ue,  painting  and  interior  finish  lead  to  widely  varying 
empty  weights. 
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The  value  ^  was  set  to  0.17  based  on 
McCormick  Figure  4.21.  The  value  k  was  found  from 
a  best  fit  to  published  Cd  versus  C/  data  [6]. 

The  parasitic  drag  of  the  remainder  of  the  aircraft 
and  the  trim  drag  were  lumped  together  and  calculated 
via 


F. 


drag,  fuse  ^t'^^equiv' 


(14) 


The  equivalent  area  was  estimated  to  be  0.3  m^  (3.225 
ft^)  by  tuning  the  model  to  best  approximate  published 
aircraft  performance. 

Although  the  performance  model  is  perfectly 
adequate  for  the  purpose  of  this  paper,  it  tends  to 
underestimate  the  aircraft’s  climb  performance. 


Exergy  Balances 

The  exergy  flow  diagrams  were  created  by 
applying  the  following  exergy  balances. 

Fuel: 

^/dt  ^  ^Jvet^fuehch  = 

^  ^  fuel  ^  fuel, ch  ~  ^S,fuel 


Engine; 

dX, 


engine  j 


^  =m  = 

^fuel^fueLch  ^ lift, engine  ^ shaft, prop 

^  Y 

^  shaft  ,alt  ^  S  ,engine 


(16) 


Propeller: 

dX 


prop/  _ 

^dt~ 


m 


propi 


g 


dY/ 


'dt ' 


X  lift  ^  prop ^  shaft,  prop  X 


thrust 


-X. 


S,prop 


(17) 


Alternator: 


Fuselage  and  Miscellaneous: 

„dY /  - 

/dt  ~  /dt 

^lifl,Juse  ^  thrust,  fitse  ~  ^  S,  fuse 


Cargo: 

^caigo/_^  gdY/  - 
/dt  "'cargo*  /(li 


^  lift. cargo  ^  • 


<5. cargo 


(21) 


Other  Relations 

As  stated  previously,  the  exergy'  of  lift  provided 
to  each  component  is  proportionate  to  the  weight  of 
each  component.  For  example,  the  exerg\'  of  lift 
provided  to  the  engine  was  found  with 

Y  _  engine  &  ^ 

^  lift, engine  ^  lift 

m  2 

gross  b  Q2) 

Likewise,  the  net  exerg>^  of  lift  delivered  by  the  wing 
was  defined  as 


V  —  F  gyving  ^  y 
lift, net  lift  ^  lift 

m  2 

^23) 

Exergy  flows  from  the  propeller  and  alternator 
were  found  with  the  following  two  equations, 
respectively. 


X  —nX 

^thrust  '  i  prop^  shaft,  prop 
^ elect  ~~  ^ alt  ^  shift, alt 


(24) 

(25) ' 


Engine  performance  was  taken  from  the 
Lycoming  operator’s  handbook,  and  propeller 
efficiency  was  found  via  Hartzell’s  performance 
software.  Alternator  efficiency  was  assumed  constant 
at  80%. 


EXAMPLE  EXERGY  FLOW  DIAGRAMS 

The  three  example  exergy  flow  diagrams  are 
given  below  in  Figures  1-3  The  fuel  is  the  source  of 
exergy  for  the  aircraft.  However,  note  that  the  fuel 
itself  requires  a  lift  exergy  input  in  the  course  of  being 
kept  aloft.  The  fuel’s  exergy  is  supplied,  naturally  to 
the  engine,  which  converts  it  to  mechanical  power. 

The  engine’s  mechanical  power  is  used  to  drive 
the  propeller,  and  also  auxiliaries,  such  as  the 
alternator.  The  propeller  creates  thrust  exergy, 
supplied  to  the  wing  to  produce  lift  and  to  overcome 
the  drag  of  the  remainder  of  the  aircraft. 

The  wing’s  exergy  of  lift  is  supplied  to  itself  and 
all  other  components  of  the  airplane  Thus,  a  portion 
of  the  fuel’s  exergy  makes  a  full  circle  journey  in 
being  converted  to  lift  to  keep  itself  in  the  air. 

For  cruise  flight,  the  difference  between  inflows 
and  outflows  of  exergy  is  the  rate  of  exergy 
destruction  in  each  component,  as  the  diagrams  are  at 
steady  state.  For  climbing  flight,  there  is  a  rate  of 
change  in  the  component’s  exergy  due  to  its  increase 
in  altitude. 

Thus,  Figures  1-3  show  quantitatively  the  rates  of 
exergy  flow  through  the  aircraft  and  the  rates  of 
destruction  within  each  component. 
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CONCLUSIONS 

If  it  is  desired  to  optimize  aircraft  energ\'  systems 
using  the  methods  of  thermoeconomics,  an  expression 
for  the  exerg\'  of  lift  is  needed.  This  is  because  a 
significant  amount  of  exergy  is  required  to  hold  the 
systems  aloft.  A  suitable  expression  has  been  derived 
here,  and  employed  to  develop  exerg>^  flow  diagrams. 
These  diagrams,  when  combined  with  mone>' 
balances,  allow  costs  to  be  associated  with  the  exerg}' 
flrrvvs.  These  costs  are  key  to  thermoeconomic 
decomposition.  (See  Reference  [1].) 

The  diagrams  here  are  based  upon  a  very  simple 
breakdown  of  the  aircraft.  A  more  detailed  aircraft 
performance  model  could  further  break  down  exerg>’ 
flows  in  the  “fuselage  and  miscellaneous”  categor>'. 

One  might  ask  the  question,  “What  happens  to 
this  lift  exerg}'  supplied  to  various  aircraft 
components?”  It  is,  of  course,  ultimately  destroyed  in 
the  process  of  holding  the  system  out  of  equilibrium 
with  the  environment.  This  is  much  the  same  as  in 
other  processes,  such  as  heating  a  house.  To  keep  a 
house  at  a  higher  (or  lower)  temperature  than  its 
environment,  one  must  continuously  supply  heat 
exergy.  This  exergy  is  continuously  “used  up”. 
Furthermore,  just  as  one  can  “optimize”  a  house  and 
its  heating  systems,  such  as  by  adding  insulation  or 
with  a  more  efficient  heating  system,  one  can 
“optimize”  an  aircraft,  by  reducing  weight  and 
matching  wings  to  the  mission. 

The  advancements  in  aerospace  technolog}’ 
suggest  future  work  along  lines  similar  to  the 
foregoing.  For  example,  at  high  speeds  the  proper 
“dead  state”  temperature  might  not  be  the  temperature 
of  the  air  surrounding  the  aircraft.  Constraints  to  heat 
transfer  with  that  air  might  require  the  dead  state 
temperature  to  be  taken  as  the  stagnation  temperature 
of  the  free  stream.  Additionally,  the  possible 


employment  of  plasma  technologies  in  fiiture 
h}personic  aircraft  suggests  a  need  for  appropriate 
property  relations  of  the  exergy  of  plasma. 
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Figure  1 :  Aircraft  Exergy  Flow,  75%  Power,  8000  Feet 


Exergy  destructions  are  shown  inside  of  boxes. 


Figure  2:  Aircraft  Exergy  Flow,  65%  Power,  8000  Feet 


,:5r.,v.  Electrical - ► 

Exergy  destructions  are  shown  inside  of  boxes. 


Figure  3:  Aircraft  Exergy  Flows,  Maximum  Rate  of  Climb  (1220  FPM),  Sea  Level 
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Rational  Design  of  Vehicles 
I.  Multi-objective  Optimization 

David  M.  Paulus,  Jr.  and  Richard  A.  Gaggioli 
Marquette  University 


ABSTRACT 

The  customer  for  a  vehicle  typically  has  several 
desiderata,  such  as  top  speed,  fuel  economy,  range, 
acceleration,  ....  Generally,  these  desiderata  are 
conflicting.  So,  in  order  to  deduce  a  single  objective 
function,  a  means  is  needed  for  weighting  (implicitly 
if  not  explicitly)  the  relative  importance  of  these 
desiderata.  That  is,  for  weighting  these  “multiple 
objectives.”  This  paper  presents  a  rational 
methodology  for  developing  a  single-objective 
function  to  be  optimized  during  the  design  of  a 
vehicle.  The  methodology  does  require  answers  from 
the  customer(s)  to  a  straightforward  set  of  questions, 
referring  to  the  desiderata.  Based  on  the  answers,  the 
objective  function  follows,  mathematically,  in  a 
straightforward  manner.  An  application  to  a  light, 
personal  aircraft  serves  as  a  case  study. 

INTRODUCTION 

This  scries  of  papers  is  an  expansion  and 
clarification  of  material  CHiginally  published  in  [1]. 

The  customer  for  a  vehicle  typically  has  several 
desiderata,  such  as  top  speed,  fuel  economy,  range, 
acceleration,  .  ...  Generally,  these  desiderata  are 
conflicting.  Consider  this  example:  the  customer  of  a 
light,  personal  aircraft  might  want  a  high  cruising 
speed,  good  rate  of  climb,  long  range  and  a  large 
payload.  Clearly,  these  desiderata  may  conflict  with 
each  other.  They  all  conflict  with  yet  another 
desideratum  -  the  minimization  of  expenditure. 

If  an  optimal  design  for  a  vehicle  is  desired,  it  is 
desirable  to  define  a  single,  “overall”  objective  from 
these  multiple,  “subsidiary”  objectives,  referred  to 
here  as  desiderata^ 

The  only  means  by  which  this  can  be  achieved  is 
to  determine  (implicitly  if  not  explicitly)  the  relative 
importance  of  the  individual  desiderata.  This  requires 
that  a  common  measure  be  found  for  the  overall 
objective  and  the  importance  of  each  desideratum,  and 
that  a  function  be  determined  that  relates  the  value  of 
the  design  to  the  desiderata. 

The  common  measure  may  be  one  of  the 
desiderata,  preferable  one  with  which  all  of  the  other 
desiderata  conflict,  typically  cost.  The  cost  may  be 
monetary,  as  would  make  sense  for  the  light  aircraft 
example,  or  it  may  be  the  number  of  total  aircraft  to  be 
produced,  a  measure  which  would  be  perhaps  of  more 
utility  for  a  military  aircraft. 


^  Some  will  argue  that  there  are  better  alternatives  tlian 
defining  a  single  objective.  We  will  address  that  issue  in  our 
conclusions. 


This  paper  introduces  an  overall  objective 
function  and  a  means  of  determining  its  component 
functions. 

THE  FUNCTION 

The  function  that  relates  the  desiderata  to  the 
value  of  the  design  will  be  referred  to  here  as 

ViPi ,  D2 , , . ,  where  the  ’s  measure  the  various 

desiderata.  The  function  to  be  optimized  takes  the 
form  of  a  “profit”  to  be  maximized, 

n  =  V(D,)-Z-c,F,-c^F^-...  (1) 

where  Z  is  the  capital  invested  (e.g.,  in  monetary'  units 
or  sacrificed  production)  the  c’s  are  the  unit  costs  of 
the  feeds,  the  F!s  in  tbe  previous  formula.  The 
obvious  feed  is  fuel,  but  other  feeds  may  be  included, 
such  as  maintenance  costs  over  the  lifecycle  of  a 
vehicle. 

DETERMINING  THE  FUNTION  V 

Clearly,  before  one  can  apply  equation  (1),  one 

must  find  the  function  F(Z),  )  .  One  would  expect 

the  function  to  return  zero  for  a  very  low  level  of 
performance  for  a  given  desideratum.  That  is,  the 
level  of  performance  is  unacceptable. 

As  an  example  of  this,  consider  an  automobile. 
A  person  residing  in  the  United  States  (with  typical 
sp^  limits  in  the  100-125  km/hr  range)  might  be 
unwilling  to  consider  purchasing  a  vehicle  incapable 
of  at  least  125  km/hr. 

However,  this  person  would  likely  be  willing  to 
pay  more  for  an  automobile  that  could  achieve  a 
higher  speed  than  this  minimal  amount.  Speed  limits 
(or  more  correctly,  the  enforcement  thereof)  cause  a 
reduction  in  the  marginal  value  of  a  1  km/hr  increase 
in  speed  as  the  top  speed  of  the  car  becomes  higher 
and  higher.  This  hypothetical  person,  therefore,  might 
not  be  willing  to  pay  any  more  money  for  an 
automobile  capable  of  reaching  240  km/hr  than  an 
automobile  capable  of  reaching  200  km/hr. 
(Acceleration  would  be  a  separate  desideratum.)  It  is 
not  rational  to  invest  further  resources  in  performance 
beyond  the  200  km/hr  “maximum”  value.” 

Conversely,  a  particular  automobile  customer 
might  state  that  no  vehicle  with  a  top  speed  below  a 
“threshold  value”  of  140  km/hr  would  be  acceptable. 

Recognizing  that  the  fimctions  have  a  threshold 
value  and  some  maximum  value,  one  would  expect  the 

flmction  to  return  a  value  of  zero  up  to  a 

certain  value  of  A,  and  become  flat  after  a  certain 
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value  for  that  desideratum  is  reached.  In  between 
these  points  (the  minimum  acceptable  and  the 
maximum  useful)  lies  some  continuous  function 

A  method  is  laid  out  below  to  estimate  the 
function.  It  consists  of  four  steps.  It  should  be  noted 
that  in  order,  ultimately,  to  optimize  a  subsystem,  it  is 
absolutely  necessary  to  develop  this  information,  that 

is,  the  function  F  (Z)^- )  ,  in  some  form,  even  if  it  is 

not  done  with  the  method  used  in  this  paper.  The 
approach  taken  in  this  paper  is,  perhaps,  the  simplest 
conceivable.  If  a  designer  determines  that  the  linear 
function  is  inappropriate  for  a  given  application,  the 
method  here  would  need  to  be  extended. 

First,  it  is  assumed  that  F(Di,Z)2,...)  can  be 

reasonably  well  represented  by  a  sum  of  functions, 
each  of  one  variable: 

V{D„D^,..)  =  V,{D,)+V^{D2)+-  (2) 

Then,  the  following  four-step  procedure  is 
proposed  for  finding  each  function  F  • 

Step  1:  Determinine  Median  Performance 

The  foregoing  algebraic  “tradeoff  functions” 

F(Z),  )  for  representing  the  value  of  performance 

may  have  an  arbitrary  shape  between  the  points  of 
minimum  acceptable  and  maximum  desirable 
performance.  However,  if  a  limited  range  of 
performance  is  considered,  the  assumption  of  a  linear 
relationship  is  reasonable.  One  way  in  which  this 
linear  function  may  be  constructed  is  around  a  median 
point.  If  a  linear  function  is  not  satisfactory,  this 
information  (minimum  acceptable,  maximum 
desirable  and  standard  performance)  will  still  be  of  use 
in  the  construction  of  the  function. 

For  a  military  combat  aircraft,  one  way  to  find  a 
median  point  is  by  considering  the  performance  of  the 
aircraft’s  adversaries,  both  current  and  projected.  For 
each  desideratum,  at  least  one  of  the  aircraft  in  the 
adversarial  group  has  a  best  value.  That  value  could 
be  selected  as  the  median  of  acceptable  performance 
levels.  The  set  of  medians  would  form  a  “standard”  of 
comparison  for  ftirther  investigations. 

The  same  basic  idea  is  applicable  to  civilian 
vehicles.  Alternatively,  the  median  values  of  a  market 
segment  might  instead  be  chosen  to  set  the  “standard 
values.” 

Step  2;  Protected  Units  Costs  and  Projected 
Production 

A  realistic  estimate  of  both  the  cost  per  vehicle 
and  the  total  production  quantity  of  the  vehicle  should 
be  made.  This  step,  and  the  following,  are  particularly 
necessary  if  using  production  figures  as  a  common 
currency. 


Step  3;  Projected  Research.  Design  and 
Development  costs 

The  projected  research,  design  and  development 
costs  should  be  listed.  When  employed  with  the 
information  attained  above,  the  total  project  cost  ma>’ 
be  estimated. 

Step  4;Aigebraic  TradeofF  Functions 

The  algebraic  tradeoff  functions,  the  F  ^or 
each  desideratum  must  be  determined  through 
questioning  of  the  customers  or  end-users.  As  a  bare 
minimum,  the  following  three  questions  should  be 
asked  (in  some  form), 

a.  What  is  the  minimum  acceptable  value  for  each 
of  the  desiderata? 

b.  Is  there  a  point  beyond  which  further 
improvement  is  not  necessary? 

c.  How  much  would  a  given  improvement,  over 
and  above  the  standard  value  determined  in 
Step  1,  be  valued? 

With  the  answers  to  these  three  questions 
determined,  the  simplest  tradeoff  function,  linear,  may 
be  determined.  Once  again  visiting  the  automobile 
example,  let  us  imagine  that  the  marketing  department 
has  asked  these  three  questions  to  potential  customers 
regarding  the  top  speed  of  an  automobile  of  a  certain 
class.  Suppose  the  average  (or  weighted  average) 
answers  were,  S’- 125  km/hr  as  a  minimum  acceptable 
top  speed,  5=200  km/hr  as  a  ceiling  beyond  which 
improvement  has  little  or  no  value,  and  a  willingness 
to  pay  800  dollars  for  an  improvement  AS-10  km/hr 
over  a  standard  162  km/hr  top  speed.  One  can 
imagine  a  function,  V(S),  which  would  have  a  value  of 
zero  dollars  up  to  125  km/hr  and  rise  with  a  slope  of 
80  dollars/(kni/hr)  to  a  maximum  of  6000  dollars. 

However,  for  many  cas  the  questions  may  not 
be  best  asked  in  terms  of  dollars.  For  a  military 
vehicle,  say  an  air  superiority  fighter,  dollars  would  be 
a  poor  choice  of  units.  In  this  example,  the  end-user 
(the  Department  of  the  Air  Force  or  Navy)  is  not  the 
same  as  the  purchaser  (Congress).  The  purchase  costs 
have  reached  such  high  amounts  that  it  is  difficult  for 
the  average  person  to  comprehend  the  sums  in  rational 
terms.  Furthermore,  neither  body  (Congress  or  the 
end-user)  is  spending  their  own  money. 

In  such  a  case  the  questions  may  be  rephrased  in 
terms  of  production  sacrifices;  Question  Number  3 
could  be  changed  to:  “What  raluction,  in  number  of 
aircraft  delivered  to  you,  would  you  accept  in  order  to 
obtain  a  specified  improvement  from  the  median  value 
determined  in  Step  1?”  A  military  leader  should  have 
a  good  grasp  of  tradeoffs  between  quantity  and 
quality.  The  information  from  Steps  2  and  3  allows 
production  tradeoffs  to  be  converted  to  a  dollar 
amount  (or  a  cost  to  a  production  adjustment). 

If  a  computer  simulation  were  to  be  available  that 
would  predict  aircraft  survivability  as  a  function  of 
measured  desiderata,  it  could  be  used  to  develop,  or 
help  develop,  these  trade-off  functions. 
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The  trade-ofF  functions  developed  in  this  manner 
may  be  used  directly  in  equation  (1).  If  they  are, 
however,  equation  (1)  will  normally  result  in  a 
negative  value.  These  trade-off  curves  do  not  take 
into  account  the  base  value  of  a  vehicle  performing  at 
the  standard  levels.  Nonetheless,  the  resulting 
fiinction  will  be  adequate  for  the  purpose  of 
comparing  competing  designs,  as  the  standard  values 
cancel  in  the  comparison. 

CASE  STUDY:  THE  FUNCTION  V  FOR  A 
LIGHT  AIRCRAFT 

The  algebraic  tradeoff  functions  were  constructed 
for  a  personal,  light  aircraft  as  part  of  a  case  study. 
As  this  aircraft  was  to  be  an  amateur-built  aircraft, 
Steps  i-3  lack  the  importance  they  would  have  for  a 
production  run  of  vehicles,  and  are  not  considered 
here. 

A  User  Questionnaire 

In  order  to  construct  the  objective  function  it  was 
necessary  to  question  the  user  about  the  value  of 
different  desiderata,  in  particular  three. 

It  was  decided,  based  on  the  light  aircraft’s 
typical  “mission”  profile,  to  consider  speed  at  full 
power  cruise  at  8000  ft.  (This  works  out  to 
approximately  75%  engine  power  for  a  normally 
aspirated  piston  engine.)  as  one  desideratum,  range  at 
economy  cruise  (65%  rated  engine  power)  at  8000  ft 
as  another,  and  finally  and  climb  rate  at  sea  level. 

The  following  questions  were  posed  to  the  end- 
user,  and  the  following  answers  received. 


The  Overall  Objective  Function 

With  this  information  linear  functions  for  the 
three  desiderata  were  created.  The  linear  functions  for 
each  desideratum  are; 


=  $0  when  <  120mph 

T.  /  ^  70$ 

K,  =(^'.«»-120mph)— 

tnph 

when  120  mph  <  <  200mph 

=  $5600  when  >  200  mph 


(3) 


-  SO  dol  when  RC<700  fpm 

F,,=(i?C-700mph)^ 

fpm  (4) 

when  700  fpm  <  RC  ^  1800  fpm 

=  $2750  when  RC>1800  fpm 

=  $0  when  R<500  miles 
$6 

=:(/?- 500  mph) - 

mile  (5) 

when  500  miles  <  R  <  1000  miles 
F^  =  $3000  when/?  >  iOOOmiles 


1 .  When  you  were  selecting  an  aircraft  kit,  what 
was  the  range  of  advertised  cruise  speeds  of  the 
aircraft  (full  throttle,  8000  ft.)  in  which  you  were 
interested?  Answer  120  to  200  mph 

2.  Approximately  how  much  additional  money 
would  you  be  willing  to  spend  for  the  aircraft  in 
order  to  increase  the  cruise  speed  of  your  airplane 
by  5  mph  above  the  average  of  the  figures  in 
question  1  ?  Answer:  $350 

3.  When  you  were  selecting  an  aircraft  kit,  what 
was  the  range  of  advertised  ranges  of  the  aircraft 
in  which  you  were  interested  (at  65%  rated 
power,  8000  ft.)?  Answer:  500  to  1000  miles 

4.  Approximately  how  much  money  would  you  be 
willing  to  spend  in  order  to  increase  the  range  of 
your  airplane  by  50  miles  over  the  average  of  the 
figures  in  question  3?  Answer:  $300 

5.  When  you  were  selecting  an  aircraft  kit,  what 
was  the  range  of  advertised  climb  rates  of  the 
aircraft  in  which  you  were  interested?  Answer: 
700  to  1800 Jpm 

6.  Approximately  how  much  money  would  you  be 
willing  to  spend  in  order  to  increase  the  climb 
rate  of  your  airplane  by  100  fpm  above  the 
average  of  the  figures  in  question  1?  Answer  : 
$250 


Here  Vc  represents  cruise  speed  at  full  throttle  and 
8000  ft.,  RC  the  rate  of  climb  at  seal  level  and  R  the 
range  at  economy  cruise  and  8000  ft. 

Now  the  objective  function,  to  be  maximized, 
may  be  written  as 

^RC  +  ^  ^Capital 

It  is  desired  to  maximize  the  value  of  the  aircraft’s 
performance  minus  the  cost  of  the  aircraft. 
Maximization  of  J  is  equivalent  to  maximization  of 
the  profit  P  shown  in  Equation  (1).  They  differ  only 
by  a  constant,  the  hypothetical  profit  of  the  standard  of 

comparison:  J^^n-Ustandard- 

Application 

The  foregoing  objective  function  was  employed 
for  two  simple  applications,  relevant  to  a  hght 
personal  aircraft.  A  kit  had  been  purchased  for 
construction  of  the  aircraft,  a  Glastar'^'^.  The  kit  is  not 
complete.  The  owner  must  complete  the  design  of  the 
airplane  by  selecting  (or,  in  theor>%  designing)  several 
component,  such  as  an  engine,  an  alterator,  avionics... 

The  Glastar'^  homebuilt  aircraft,  developed  by 
Arlington  Aircraft  Development,  Inc,  is  typically 
powered  by  a  Lycoming  0-320  or  0-360  engine 
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prcxlucing  between  150  and  200  hp,  turning  either  a 
constant-speed  or  fixed  pitch  propeller.  The  aircraft 
has  a  gross  weight  of  1960  pounds.  The  Glastar  has  a 
NASA  GAW-2  airfoil,  a  wingspan  of  35  feet,  and  an 
aspect  ration  of  9.6.  For  the  purpose  of  this  article,  the 
“standard”  aircraft  will  be  assumed  to  be  powered  by  a 
160  hp  0-320  engine  with  a  constant-speed  Hartzell 
propeller  and  to  have  an  empty  weight  of  1200 
pounds. 

For  all  calculations,  it  was  assumed  that  the 
aircraft  would  have  a  lifespan  of  20  years,  be  flown 
200  hours  per  year,  with  50%  of  those  hours  at 
economy  cruise  and  the  remaining  50%  at  high-speed 
cruise.  Hours  spent  in  climb  are  considered 
negligible.  An  interest  rate  of  8%  was  used. 

Alternator  Selection 

First  it  was  decided  whether  it  is  better  to 
purchase  a  “standard”  alternator  or  a  lightweight 
model  for  the  light  aircraft.  As  no  efficiency  data  is 
available,  both  alternators  were  assumed  equally 
efficient.  The  standard  alternator  has  an  initial  cost  of 
$294  and  a  mass  of  5.9-kg  (13-lbm).  The  alternative 
costs  $450,  but  has  a  mass  of  only  2.7-kg  (6-lbm). 

Performance  of  the  aircraft  was  simulated  with 
both  alternators  using  a  simple  flight  model. ^  The 
developed  objective  function  predicts  a  decrease  in 
present  value  profit  of  $118  by  choosing  the 
lightweight  alternator. 

Engine  Selection 

Now,  the  objective  function  is  used  to  choose 
between  two  possible  engines.  The  160-hp  engine  is 
the  “standard”  engine,  chosen  most  widely  by 
builders.  The  180-hp  engine  will  produce  gains  in  fiill 
power  cruise  and  climb  rate,  but  at  the  expense  of 
indexed  weight,  decreased  range,  increased  fuel 
consumption  and  increased  capital  cost.  Table  1 
shows  predicted  performance  with  both  engines, 
estimated  with  the  previous  flight  model. 


Engine  (hp) 

Cost  (k 
dollars) 

Mass  (kg) 

Full  Power 
Cruise 
(km/hr)/ 

Fuel  Bum 
(kg/hr) 

Fuel  Bum  @ 
244  km/hr 
(kgflir) 

Rate  of  Climb 
(m/min) 

160 

22.3 

116 

259/22.9 

18.8 

372 

180 

26.2 

122 

271/29.9 

21.8 

442 

Table  1:  Performance  with  different  engines 


The  objective  function  shows  a  $10,800  decrease 
in  present-value  profit  with  the  more  powerful  engine. 

CONCLUSION 

A  rational  methodology  for  developing  a  single¬ 
objective  function  for  a  vehicle  has  been  presented 


here.  It  uses  information  firom  the  customer,  who  is 
the  ultimate  judge  of  trade-offs  in  cost  and 
performance,  to  develop  this  function.  In  Part  2,  this 
information  will  be  us^,  along  with  thermoeconomic 
methods,  to  optimally  select  vehicular  energy* 
conversion  equipment. 

It  should  be  mentioned  that  two  aspects  of  the 
foregoing  development  were,  incidentalh,  simplistic 
and  could  be  generalized.  One  was  that  the  “user 
questionnaire”  was  addressed  to  only  a  single 
customer;  this  could  easily  be  broadened  b\’  using 
modem,  statistical  market  analysis  methods.  Also,  the 
function  F(Z)/,  D?,  .  .  .)  was  assumed  to  be  a  sum  of 
single-variable  fimetions  of  each  /),,  and  these 
functions  were  represented  linearly.  With  a  more 
elaborate  battery  of  questions,  higher  order 
representations  could  be  emplo>^ed,  including 
“coupling”  terras  involving  more  than  one  of  the  D,. 

Another  point  needs  to  be  made,  with  reference 
to  Footnote  1  presented  earlier.  Some  might  argue 
that  it  is  not  necessary  to  define  a  single  overall 
objective.  (Although  we  did  not  say  that  it  was 
necessary,  but  “desirable.”)  Our  value  function  V{Pi, 

.  .  0  is  a  means  for  weighting  the  relative 
importance  of  the  different  desiderata.  In  order  to 
finally  determine  an  optimum,  the  so-called  “multi¬ 
objective  optimization”  methodologies  (e.g..  Fan  and 
Shieh,  1983)  must  and  do,  ultimately,  invoke 
weighting  factors.  Thereby,  they  establish  the 
equivalent  of  a  single  objective.  Concerns  that  use  of 
a  single-objective  approach  like  ours  could  lead  to 
finding  a  local  optimum  and  not  the  global  optimum, 
would  be  tantamount  to  suggesting  that  single 
objective  functions  should  not  be  used  even  when 
there  is  only  one  objective.  Moreover,  in  general  it  is 
easier  to  cany  out  an  appropriate  global  search  than  to 
employ  such  techniques  as  the  generation  of  tradeoff 
(Pareto)  and/or  indifference  “cxirves,”  especially  in 
those  instances  (usual)  when  there  are  more  than  two 
objectives  -  when  the  curves  become  surfaces  and 
hypersurfaces.  Furthermore,  the  value  function  V  is 
determined  by  a  rational  procedure  before  the 
optimization  process  is  begun,  avoiding  subjectivity 
that  could  enter  when  making  “tradeoffs”  toward  the 
end  of  optimization  procedures.  Finally,  we  believe 
that  with  the  approach  advocated  it  is  easier  to  keep 
the  distinction  clear  between  “desiderata”  (desirable 
features  of  the  product)  and  “objectives”  (to  be 
extremized). 
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ABSTRACT 

Overall  vehicle  design  entails  the  subsidiary 
design  -  concurrent  engineering  -  of  the  numerous 
subsystems  that  make  up  the  vehicle.  The  engineer(s) 
in  charge  of  the  overall  design  must  coordinate  the 
work  and  integrate  the  designs  of  several 
‘‘independent'’  subsystem  design  teams.  That 
coordination  and  integration  is  achieved  by 
specification  of  a  particular  set  of  constraints,  imposed 
upon  each  design  team.  As  the  overall  design 
proceeds,  as  a  result  of  feedback  from  and  review  of 
the  work  of  the  individual  teams,  these  constraints 
may  be  adjusted.  Then  many  of  the  subsystems  need 
to  be  re-designed.  This  procedure  is  carried  out, 
iteratively,  until  “convergence”  is  achieved.  That  is, 
until  the  overall  design  is  “suitable  and  satisfactory”  in 
the  judgment  of  the  engineer  in  charge  of  the  overall 
design,  in  consultation  -  directly  or  indirectly  -  with 
the  “customer(s).”  Better  yet,  die  iterative  procedure 
is  continued  until  a  rational  objective  function  for  the 
overall  vehicle  is  optimized.  In  order  to  facilitate  this 
concurrent  engineering  process,  it  is  desirable  that  the 
design  decisions,  being  made  by  the  subsystem  design 
teams,  explicitly  take  into  account  the  customers’ 
weighted  desiderata  for  the  overall  system.  In  this 
paper,  it  is  shown  how  the  methods  of 
thermoeconomic  decomposition  can  be  applied  (in 
conjunction  with  the  rational  objective  function  for  the 
overall  design  presented  in  Part  I)  to  deduce  objective 
functions  for  optimizing  the  design  of  energy- 
conversion  subsystems.  An  illustrative  case  study  is 
presented,  for  the  design-selection  of  an  engine  and  a 
generator  for  a  light  aircraft. 

INTRODUCTION 

The  design  method  to  be  presented  here  depends 
upon  the  initial  development  of  a  preliminary 
conceptual  design. 

Once  a  conceptual  design  has  been  accepted,  the 
detailed  design  should  proceed  in  an  efficient  fashion. 
It  is  not  possible  to  optimize  and  design  a  vehicle,  or 
even  an  energy  subsystem,  as  an  entirety.  Vehicular 
energy  systems  are  too  complex  to  be  designed  by  a 
single  individual  or  team,  as  the  number  of  decision 
variables  becomes  unmanageable.  Therefore  the 
design  or  selection  of  individual  devices  and/or 
subsystems  is  delegated  to  subordinate  teams  or 
individuals.  It  is  desirable  for  these  designers  to  have 
a  methodology  and  information  that  allows  them  to 
make  decisions  in  accordance  with  the  overall  goals  of 
the  vehicle.  ^ 


*  Naturally,  it  must  be  determined  which  decision  variables 
are  local  and  which  are  global.  That  is,  which  variables  a 


A  vehicle  typically  relies  on  one  fuel  to  achieve 
its  goals.  Because  the  subs3'stem  or  de\’ice  design 
team  is  not  optimizing  a  whole  vehicle,  but  something 
that  may  produce  or  consume  commodities  not 
normally  considered  when  looking  at  a  whole  vehicle, 
its  objective  function  will  vary'  from  that  of  an  overall 
vehicle.  For  example,  the  alternator  on  a  light  general 
aviation  aircraft  does  not  directly  consume  fuel,  nor 
does  it  directly  influence  performance.  Nonetheless, 
an  aircraft  wiA  an  alternator  that  is  both  lighter  and 
more  efficient,  will,  with  all  else  remaining  equal, 
perform  better  (and  cost  more).  The  person  designing 
(or  selecting)  an  alternator  should  not  be  burdened 
with,  nor  at  this  stage  of  the  design  be  necessarily 
capable  of,  directly  calculating  the  impact  upon 
aircraft  performance 

The  alternator  mentioned  above  still  has  only  one 
product  (electrical  energy).  An  onboard  energy^  device 
or  subsystem  may  use  several  “fuels”  and/or  supply 
several  “products”.  These  fuels  and  products,  besides 
electricity  and  shaft  power,  include,  but  are  not  limited 
to,  compressed  air,  hydraulic  power  and  heat  (or 
cooling).  Additionally,  for  aircraft,  one  must  account 
for  the  lift  required  to  hold  the  device  in  the  air  (or  to 
make  it  go  up). 

ASSUMPTIONS 

The  techniques  presented  here  are  for 
employment  after  a  conceptual  design  has  been 
penned.  The  (conceptual  design)  engineer(s)  has 
estimated  the  weight,  power  and  exergy  consumption 
of  the  energy  conversion  systems.  Therefore,  only 
changes  from  the  conceptual  design  are  analyzed,  as 
the  “black  boxes”  are  “filled  in”  in  detail. 

In  this  paper,  when  monetary  values  needed  to  be 
converted  from  present  to  time  vdues,  or  vice-versa,  it 
was  assumed  that  the  aircraft  would  have  a  lifespan  of 
20  years,  be  flown  200  hours  per  year,  with  50%  of 
those  hours  at  economy  cruise  and  the  remaining  50% 
at  high-speed  cruise.  Hours  spent  in  climb  are 
considered  negligible.  An  interest  rate  of  8%  was 
used. 

THE  DECOMPOSED  OBJECTIVE 
FUNCTION 

In  general,  a  vehicular  energy  subsystem  will  fall 
into  one  of  two  categories.  It  may  be  an  auxiliary.  In 
the  context  of  this  paper,  such  a  system  is  necessary 


subsystem  design  team  is  free  to  vary'  in  optimizing  its 
individual  subsystem,  and  which  are  constraints  to  the  design 
team  so  as  not  to  significantly  affect  the  designs  being  carried 
out,  simultaneously,  by  other  design  teams.  These 
determinations  may  at  times  be  made  through  “common 
sense”,  but  at  other  times  may  require  a  sensitivity  analysis. 
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(or  desirable)  for  the  operation  of  a  vehicle,  but  the 
system  does  not  contribute  to  any  of  the  various 
desiderata,  it  can  only  detract,  by  adding  weight  to  the 
vehicle  or  consuming  exerg>\  Examples  include 
components  of  the  charging  system  and  HVAC 
systems.  Here  such  systems  will  be  referred  to  as 
non-propulsive  subsystems.  Other  on-board  energy 
systems,  notably  components  of  the  propulsion 
system,  directly  contribute  to  achieving  the  desiderata. 
The  form  of  the  objective  function  for  a  component 
will  thus  take  one  of  two  forms,  dependant  upon 
which  function  the  component  has. 


The  Objective  Function  for  an  On«board  Non* 
proDuisive  Subsystem 


An  on-board  subsystem  is  required,  in  general,  to 
deliver  a  certain  output  as  a  design  constraint.  The 
goal  in  the  design  of  such  a  subsystem  is  then  to 
minimize  its  contribution  to  the  total  cost  of  the 
vehicle  (which  will  maximize  the  “profit”  of  the 
vehicle;  see  Part  I).  The  subsystem  objective  function, 
to  be  minimized,  takes  the  form 


j  =  z+YJVi 


w + y  a-  ,Xi  ,■ 


(1) 


'■  V  J  J 

Here  i  represents  each  desideratum.  Wi  measures 
the  relative  importance  of  each  desideratum.  Average 
unit  costs  (C ’s)  of  either  exergy  or  mass  (subscript 
m)  are  evaluated  with  the  vehicle  operating  in  a  mode 
where  the  individual  desideratum  is  measured  (e.g., 
full  throttle  at  a  specified  altitude).  The  subscript  j 
represents  an  exergy  feed.  For  reasons  described 
below,  it  is  more  convenient  to  associate  a  unit  cost 
with  mass  than  with  lift  exergy.  (This  is  of  course 
only  relevant  to  aircraft.)  Therefore  the  costs 
associated  with  lift  exergy  are  accounted  for  with  the 


mass  cost,  . 


From  (1),  an  approximate  expression  for  the 
change  in  the  total  cost  of  a  vehicle  due  to  a  change  in 
subsystem  design  may  be  written  as 


,Am  +  yc.^.AZ.^ 


(2) 


This  assumes  relatively  small  changes  in  Z,  m  and  X . 
The  c’s  now  represent  marginal  exergy  costs. 


The  Objective  Function  for  a  Propulsion 
Subsystem  Component 

While  a  non-propulsive  subsystem  typically  has 
its  output  constrained,  the  propulsion  system  itself  is  a 
balance  between  the  performance  it  supplies  to  the 
vehicle  (e.g.,  as  thrust  to  an  aircraft)  and  its  capital 
cost,  weight  and  exergy  consumption.  Therefore, 
instead  of  only  minimizing  costs  to  achieve  an  optimal 
design,  it  is  necessary  to  maximize  its  profit  -  the 
difference  between  performance  delivered  and  the  cost 
of  providing  it.  The  change  in  the  overall  vehicle’s 
profit  may  be  approximated  with 


(  \ 

'  V  J 


(3) 


i  k 

The  j  represents  an  exergy'  output,  k  an  exergy’  input. 
Two  new  symbols  appear  in  this  function.  The 
marginal  unit  revenue  of  an  exergy*  output  is 
represented  with  r.  The  subscript  c  distinguishes  this 
components  “profit”  fi-om  the  profit  of  the  vehicle  as  a 
whole. 

A  time  factor,  J,  is  also  now  necessary^  It 
represents  the  fraction  of  a  vehicle’s  time  spent 
operating  in  the  regime  where  an  individual 
desideratum  is  measured.  The  propulsion  system 
component  is  only  charged  for  its  projected  inputs.  As 
an  example,  the  top  speed  of  an  automobile  may  be 
quite  important  to  the  purchaser  of  a  sports  car,  but  the 
vehicle  will  spend  little  of  its  lifetime  at  full  throttle, 
say  0.1%.  This  automobile  will  only  see  the  high  fuel 
flow  associated  with  this  operation  0.1%  of  the  time. 

The  methods  for  finding  the  above-mentioned 
weighting  factors  and  unit  costs  will  now  be 
explained.  Creation  of  exergy  flow  diagrams  for  light 
aircraft  is  explained  in  detail  by  Paulus  and  Gaggioli 
[1]. 


DETERMINING  THE  WEIGHTING 
FACTORS 

In  Part  I.  of  this  series,  an  overall  objective 
function  was  presented.  For  generic,  global 
optimization  of  a  vehicle,  this  is  sufficient.  However, 
for  concurrent  detailed  design  to  take  place,  and 
equations  (1)  and  (3)  to  be  employed,  weighting 
factors  are  needed. 

The  relative  weighting  factor  for  a  desideratum 
,  as  employed  in  this  papr.t ,  is  deduced  from  the 

algebraic  tradeoff  functions  presented  in  Part  I.  It  is 
given  by 


Where  F)  represents  the  value  returned  by  the 
tradeoff  function  VfDJ,  Z),  a  performance 
desideratum,  and  AF;  the  increase  in  price  an  end  user 
would  pay  for  a  AD/  increase  in  performance.  The 
numerator,  then,  represents  the  percent  of  expenditure 
increase  the  customer  is  willing  to  make,  per  percent 
increase  in  performance  desideratum  D,.  The 
denominator  is  the  sum  over  all  performance 
desiderata;  so  represents  the  relative  importance  of 
D/.  F  cancels,  leaving 
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/AD, 


/D, 


(5) 


DETERMINING  THE  UNIT  COSTS 


The  proper  unit  exergy  costs 

As  mentioned  above,  it  is  assumed  that  before 
detailed  component  design  takes  place,  a  “generic” 
design  of  the  overall  vehicle  is  present.  Exergy  and 
other  unit  costs  to  be  used  in  Equations  (1)  and  (2)  for 
conceptual  design  are  derive^  at  first,  from  the 
conceptual  design.  During  detailed  design  of  a 
component,  the  “filling  in”  of  “black  boxes”,  the 
component  design  engineer  has  no  immediate  effect 
on  equipment  specifications  outside  of  his  “box”. 
(Obviously,  the  engineer’s  decisions  will  affect,  to  at 
least  some  extent,  other  equipment’s  performance.  If 
need  be,  one  can  account  for  those  effects  in  a  second 
iteration  of  the  overall  component  designs.  And  so  on 
and  so  forth.)  Therefore,  the  proper  exergy  costs  (and 
values)  to  employ  are  marginal  costs.  In  some  cases 
these  may  be  estimated  by  computing  average  unit 
exergy  costs  without  consideration  of  capital  costs.  In 
other  instances  (e.g.,  where  a  jump  to  afterburning  in  a 
jet  engine  is  necessary  for  increased  thrust,  or  when 
neither  a  system’s  input  nor  its  output  is  fixed), 
however,  such  methodology  will  result  in  error  due  to 
the  average  costs  not  approximating  marginal  costs 
[2]. 

Because  the  lift  exerg>'  requirements  for  an 
aircraft  go  up  with  the  square  of  its  weight,  unit  costs 
for  lift  exergy  are  an  area  where  the  use  of  average 
unit  costs  without  consideration  of  capital  costs  will 
result  in  error.  For  this  reason,  marginal  exergy  costs 
for  lift  are  instead  accounted  for  in  (1)  and  (3)  as  unit 
mass  costs,  derived  below. 

Fixing  Marginal  Costs  and  Revenues 

If  marginal  exergy  costs  and/or  revenues  may  be 
estimated  by  calculating  average  costs  while 
neglecting  capital  costs,  they  are  found  by  applying 
money  balances  to  each  component  of  the  vehicle. 
The  end  goal  of  these  money  balances  is  to  find  a 
marginal  unit  cost  and/or  revenue  of  each  exergy 
stream.  The  “source”  of  these  costs  is  generally  ^e 
vehicle’s  fuel,  while  the  revenues  result  from 
vehicular  performance. 

While  applying  Equation  (1)  to  non-propuisive 
subsystems,  necessary  auxiliary  relations  fixing  a  cost 
of  an  exergy  stream  the  exergy  costing  for  each 
desideratum  Dj  may  be  accomplished  in  one  of  two 
ways,  depending  on  the  nature  of  the  desideratum. 

Certain  desiderata  are  measured  at  full  engine 
power.  Changes  in  subsystem  exergy  draw  will  not 
affect  the  overall  fuel  flow  rate.  It  will,  however, 


certainly  affect  the  vehicle’s  performance.  Therefore, 
for  the  case  of  desideratimi  ,  unit  exerg>'  cost  C^j 


of  a  feed  X^j  should  be  calculated  by  determining  the 

marginal  cost  of  the  exerg\  flow  which  has  the 
greatest  direct  impact  on  the  desideratum.  For 
example,  this  force  is  thrust  exerg>'  for  an  aircraft. 
This  cost  is  foimd  with 


%■ 


dU  dP. 
dD,  dX,j 


(6) 


Here  represents  the  performance  desideratum  and 


n  the  overall  vehicle  profit,  (See  Part  I  of  this  paper.) 

Other  desiderata  are  measured  at  part  throttle. 
Here,  performance  may  be  kept  constant  by  varying 
engine  power  output  while  varying  subsystem  exergy^ 
draw.  Then,  the  consequential  exergy'  cost  of  the  fuel 
may  be  assigned. 

Changes  in  subsystem  exergy’  draw  effect  not 
only  speed,  etc.  but  also  range.  Therefore,  when  range 
is  the  considered  desideratum  (A)  not  only  must  the 
direct  cost  of  fuel  be  considered,  but  also  the  indirect 
cost  of  the  increased  fuel  bum  on  the  range  R  of  the 
vehicle.  This  can  be  accomplished  by  adding  the 
following  term  to  the  base  cost  of  fuel. 

an  dR 


^Juel^range 


dR  dm 


(7) 


■fuel 


Although  this  marginal  cost  is  shown  per  unit 
mass,  it  could  also  be  calculated  on  a  per  unit  exergy 
basis. 

Components  in  the  propulsion  or  drivetrain 
system  of  a  vehicle  both  consume  exergy%  and  supply 
exergy  that  drive  vehicle  performance.  Hence  both 
terms  for  unit  exergy'  cost  and  revenue  appear  in 
equation  (3).  In  this  equation,  marginal  costs  (c’s), 
may  be  determined  from  the  known  unit  price  of  fuel 
(adding  the  cost  to  range  if  necessary)  while  revenues 
(r’s)  should  be  determined  in  the  same  manner  as  with 
equation  (6),  with 


50  5D,. 
aD,  dX, 


(8) 


The  Marginal  Cost  of  Mass 


The  marginal  cost  of  mass,  appearing  in  both 
Equations  (1)  and  (2),  may  be  found  with 


(9) 


where  is  the  exergy  of  lift  and  its  unit  cost. 

With  the  exergy  of  lift  devoted  to  a  component  of 
mass  m^  given  by  [1] 


^lift  = 


napAV^^ 


(10) 
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equation  (8)  becomes,  when  applied  to  the  whole 


aircraft, 


V 


KapAV^  ^ 


^lift 


(11) 


CASE  STUDY:  A  LIGHT  PERSONAL 
Aircraft 

These  methods  were  apphed  to  the  Giastar 
aircraft  described  in  Part  1.  to  select  optimally  an  on¬ 
board  non-propulsive  subsystem,  and  to  optimally 
select  an  engine. 

The  Wei2htine  and  Time  Factors 

First,  weighting  factors  are  calculated  for  the 
three  desiderata,  described  in  Part  I.:  maximum  cruise 
speed  at  8000  ft.,  range  at  65%  rated  engine  power  and 
8000  ft.  and  maximum  sea-level  rate  of  climb.  The 
tradeoff  functions  were  given  in  Figures  1-3  of  Part  I. 
Using  equation  (4),  the  weighting  factors  were 
calculated  to  have  the  values  in  Table  1. 

As  stated  in  Part  1,  the  aircraft  was  assumed  to 
spend  50%  of  its  flight  time  in  maximum-speed  cruise 
(full  throttle,  i.e.,  75%  rated  power,  at  8000  ft.)  and 
50%  of  its  flight  time  in  economy  cruise  (65%  rated 
power  at  8000  ft.).  Time  spent  in  climb  was 
considered  negligible.  These  assumptions  \ield  time 
factors  (7))  of  0.5, 0.5  and  0.0. 


fuel  .total  ^ fuel  ^ 


(12) 


Engine: 

^ fuel. total  ^ fuel. ch  ^ lift  ^ lift  .engine 


-C, 


shcft 


^^shaft.^ 


4.  Y 

prop  ^  shaft. alt 

(13) 


)=o 


Propeller: 

^ shaft^ shaft  .prop  "1"  ^ lift  ^ lift,  prep 

— F  Y  —  0 

^thrust^  thrust  ““  ^ 

Alternator: 

^shaft^ shaft. alt  ^lift^ lift. alt 

^^elect^ elect  ^ 

Wing: 

^thrust ^thrust. wing  ^ lift ^ lift, net  “  ^ 

Fuselage  and  Miscellaneous: 

^ thrust^ thrust, misc  ^ lift ^ lift. misc  ~~  ^misc 


(14) 


(15) 


(16) 


(17) 


Desideratum 

Weighting  Factor 

Cruise  Speed 

0.565 

Range 

0.283 

Climb  Rate 

0.153 

Table  1:  Weighting  Factors 


Txergv  Costing 

For  the  case  of  a  light  aircraft,  it  was  assumed 
that  marginal  exergy  costs  could  be  approximated  with 
average  exergy  costs,  which  were  calculated  from 
money  balances  while  ignoring  capital  costs.  When 
finding  these  costs,  all  lift  exergy  and  all  shaft  power 
was  assigned  the  same  cost.  In  climb,  stored  exergy 
(due  to  the  available  energy  of  the  aircraft  and 
gravitational  field)  was  assigned  a  value  of  zero.  Fuel 
was  assumed  to  have  a  cost  of  two  dollars  per  gallon,  a 
density  of  six  pounds  per  gallon  and  have  a  chemical 
exergy  content  of  46,000  kJ/kg. 

Balances  The  following  seven  money  balances  were 
used,  one  for  each  major  airplane  component  (shown 
on  the  exergy  flow  diagrams  in  a  companion  paper, 
“Exergy  of  Lift,  and  Aircraft  Exergy  Flow  Diagrams”, 
Paulus  and  Gaggioli,  this  volume). 

Fuel: 


Cargo: 

^ lift  ^ lift  .cargo  “  ^ cargo  (^^) 


The  C’s  are  the  rate  at  which  money  is  flowing 
to  a  component  to  keep  it  aloft  with  the  rest  of  the 
aircraft.  That  is,  for  the  cargo,  the  cost  of  transporting 
it  and  for  the  fuselage,  the  cost  of  having  a  place  to  put 
the  cargo,  etc. 

For  each  of  the  modes  of  operation  (full  power 
cruise,  economy  cruise  and  climb)  the  exergy  flow 
rates  were  taken  from  the  exergy  flow  diagrams  in  the 
companion  paper.  Then  (for  each  mode)  there  are 

eight  quantities,  C ’s  and  C ’s,  to  be  determined  with 
the  help  of  these  seven  equations,  (1 1)-(17).  One  of 
the  eight  needs  to  be  and  can  be  assigned  rationally, 
for  each  mode. 

Assigned  Unit  Costs  During  maximum  speed  cruise 
and  sea  level  climb,  the  engine  is  operating  at  full 
throttle.  Therefore,  a  cost  of  thrust  exergy  must  be 
found  with  equation  (6),  using  the  aircraft  flight  model 
(see  [1])  and  the  algebraic  trade-off  functions  shown 
in  Figures  1  and  2  in  Part  I.  The  partial  derivatives  of 
cruise  speed  and  climb  rate  with  respect  to  thrust 
exergy  were  found  to  be  0.381  m/s-kW  and  0. 1 15  m/s- 
kW  respectively.  This,  combined  with  marginal  profit 
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$/ 

for  speed  of  2,22  10’^  from  Figure  1  of  Part  I  and 

m/ 

/s 


a  marginal  profit  of  climb  rate  of  5,57  10"^ 


—  from 
m/ 

/s 


Figure  2  of  Part  I  yielded  unit  exergy  costs  of  thrust  of 
8.43  10“^  $/kJ  for  cruise  and  6.39  lO"^  $/kJ  for  climb. 

Range  is  evaluated  at  65%  rated  power  at  8000 
feet,  a  part  throttle  condition.  'P'e  range  of  a  propeller 
driven  aircraft  is  well  approximated  with 
(McCormick,  1995) 


Applying  equation  (7)  to  equation  (19) ,  using  Figure 

an 

3  from  Part  I  to  evaluate - yields  an  added  cost  to 

dR 

fuel  (for  decreasing  the  range  of  the  aircraft)  of  0.414 
$/gal. 


Exergy  Costs  Equations  (12)  through  (18)  were 
solved  simultaneously  using  the  unit  thrust  exergy  or 
fuel  costs  given  above  to  approximate  the  marginal 
unit  exergy  costs  for  full  throttle  cruise,  climb  and 
economy  cruise. 

Solving  these  equations  simultaneously  does  not 
result  in  a  good  approximation  for  the  marginal  cost  of 
fuel  in  a  full-throttle  condition.  For  the  full-throttle 
cases,  the  marginal  cost  of  fuel  (exergy  basis)  is  found 
from  the  price  of  fuel. 

The  approximate  marginal  exergy  costs  are  given 
in  Table  2. 


FuII- 

Throttle 

Cruise 

Climb 

Economy 

Cruise 

Lift 

4.26E-5 

2.01E-6 

3,39E-4 

Thrust 

8.43E‘-6 

6.94E-6 

s.oeE-s 

Shaft  Power 

6.80E-6 

4.71E-6 

6.51E-5 

Fuel 

7.24E-6 

7.24E-6 

1.98E-5 

Electric 

1.09E-4 

7.86E-6 

2.14E~4 

Table  2;  Approximate  Marginal  Exergy  Costs  for 
the  Light  Aircraft,  $/kJ 


Alternator  Selection 

Here  it  will  be  decided  whether  it  is  better  to 
purchase  a  “standard”  alternator  or  a  lightweight 
model  for  the  light  aircraft.  As  no  efficiency  data  is 
available,  both  will  be  assumed  equally  efficient.  The 
standard  alternator  has  an  initial  cost  of  $294  and  a 
mass  of  5.9-kg  (13-lbm).  The  alternative  costs  $450, 
but  has  a  mass  of  only  2.7-kg  (6-lbm). 

To  make  this  decision.  Equation  (1),  for  a  non- 
propulsive  subsystem,  is  used.  Terms  for  shaft  power 
are  unnecessary,  as  the  power  inputs  to  both  units  are 


the  same.  The  cost  of  switching  from  the  standard  to 
the  lightweight  alternator  will  be  given  by 

~  ^lightweight  ^standard 

{^lightweight  ""  ^standard) 
i 

The  individual  terms  have  been  given  in  tables  1 
and  2.  Equation  (20)  predicts  that  the  lightweight 
alternator  will  have  a  present  value  cost  of  $1 18  more 
than  the  standard.  The  standard  alternator  is  the  better 
choice. 

This  decision  could  be  made  by  alternate  means, 
directly  employing  the  aircraft  flight  model  and  using 
equation  2,  as  was  done  in  Part  I.  This  method 
predicted  an  added  cost  of  $121.  The  $118  value 
computed  using  the  decomposed  optimization  of  the 
alternator  compares  favorably  to  this  $121  value. 


Engine  Selection 


A  160-hp  engine  is  the  “standard”  engine,  chosen 
most  widely  by  builders.  A  180-hp  engine  will 
produce  gains  in  full  power  cruise  and  climb  rate,  but 
at  the  expense  of  increased  weight,  decreased  range, 
increased  fuel  consumption  and  increased  capital  cost. 
(See  Part  I.). 

Equation  (3)  for  a  propulsive  subsystem  may  be 
employed  in  the  form 


-(2,„ 


Zv  (x  -  X 


(21) 


The  unit  values  of  power  (i.e.,  exergy)  may  be 
taken  from  the  exergy’  cost  table,  Table  2  for  full- 
throttle  cruise  and  the  rate  of  climb,  as  the  calculation 
methods  are  identical.  For  the  part  throttle  flight  at 
economy  cruise,  additional  power  capacity  has  no 
value,  as  speed,  and  therefore  power  output,  is  held 
constant.  The  weighting  factors  come  from  Table  1. 
The  time  factors,  as  stated  previously,  are  0.5  for  full 
power  cruise,  0.5  for  economy  cruise  and  0  for  climb. 
Fuel  costs  are  calculated  on  a  per  unit  mass  basis  here, 
as  fiiel  is  the  raw  energy  source.  The  value  used  is 
0,735  dollars/kg.  (Exergy  costing  with  the  fuel  exergy 
costs  assigned  would  yield  the  same  values.)  The  cost 
of  fiiel  bum  to  range  must  be  added  to  the  cost  of  fuel 
for  evaluation  of  range  at  economy  cruise,  as  per 
equation  (7),  increasing  the  per  kilogram  cost  of  fuel 
to  0.912  dollars. 

Evaluation  of  equation  (21)  shows  a  decrease  of 
present  value  profit  of  11,400  dollars  for  the  180-hp 
engine. 

Alternatively,  the  method  used  in  Part  I  directly 
employ  the  overall  objective  function.  This  results  in 
decrease  of  profit  (present  value)  of  10,800  dollars. 

Although  this  differs  from  the  11,400-doUar 
decomposed  value  by  5.5%,  the  same  decision  would 
be  reached  in  either  case.  The  error  likely  results  from 
the  large  shift  from  the  preliminary  design  values, 


shifts  in  engine  power  of  12%,  (engine)  weight  of  6% 
and  fuel  consumption  of  14%.  In  such  cases  iteration 
in  the  exergy  costing  step  would  improve  results  if 
necessary. 

CONCLUSION 

The  combination  of  a  rational  objective  fimction 
for  a  vehicle  and  the  methods  of  thermoeconomics 
were  used  to  decompose  a  vehicle  into  systems  and 
subsystems  in  orders  to  make  design  decisions.  These 
decisions  were  made  with  the  pnie^ted  objective 
ftinctions  for  a  non-propulsive  .  or  propulsive 
subsystem. 

The  results  from  these  objective  functions 
showed  excellent  agreement  with  the  results  of  whole- 
vehicle  simulation  in  combination  with  the  overall 
objective  function  for  a  vehicle.  The  main  conclusion, 
then,  is  that  decomposition  of  a  vehicle,  for  the 
purpose  of  effective  concurrent  engineering  of 
subsystems,  can  evidently  be  achieved.  In  addition  to 
exergetic  costing,  the  key  is  the  utilization  of  a  rational 
objective  function  for  the  overall  vehicle,  to  get  the 
weighting  factors  needed  for  subsystem  objective 
functions. 

The  next  step  in  testing  the  methodology 
presented  here  is  to  apply  it  to  subsystem  design, 
instead  of  only  subsystem  selection.  Then,  the 
application  to  more  complex  vehicles  could  be 
pursued. 
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